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Summary 



Solar coronal mass ejections (CMEs) are large-scale eruptions of plasma 
and magnetic field from the Sun into the corona and interplanetary 
space. They are the most significant drivers of adverse space weather 
at Earth and other locations in the heliosphere, so it is important 
to understand the physics governing their eruption and propagation. 
However the diffuse morphology and transient nature of CMEs makes 
them difficult to identify and track using traditional image process- 
ing techniques. Furthermore, the true three-dimensional geometry 
of CMEs has remained elusive due to the limitations of coronagraph 
plane-of-sky images with restricted fields-of-view. For these reasons 
the Solar Terrestrial Relations Observatory (STEREO) was launched 
as a twin-spacecraft mission to fly in orbits ahead and behind the 
Earth in order to triangulate independent observations of CME struc- 
ture. It is the first time CMEs have been observed from vantage points 
off the Sun-Earth line and each spacecraft carries an instrument suite 
designed to image from the low solar corona out to the orbit of Earth 
in order to observe and study CME propagation towards Earth. 

In this thesis the implementation of multiscale image processing tech- 
niques to identify and track the CME front through coronagraph im- 
ages is detailed. An ellipse characterisation of the CME front is used 
to determine the CME kinematics and morphology with increased 
precision as compared to techniques used in current CME catalogues, 
and efforts are underway to automate this procedure for applying to 
a large number of CME observations for future analysis. It was found 
that CMEs do not simply undergo constant acceleration, but rather 
tend to show a higher acceleration early in their propagation. The 



angular width of CMEs was also found to change as they propagate, 
normally increasing with height from the Sun. However these results 
were derived from plane-of-sky measurements with no correction for 
how the true CME geometry and direction affect the kinematics and 
morphology observed. 

With the advent of the unique dual perspectives of the STEREO 
spacecraft, the multiscale methods were extended to an elliptical tie- 
pointing technique in order reconstruct the front of a CME in three- 
dimensions. Applying this technique to the Earth-directed CME of 
12 December 2008 allowed an accurate determination of its true kine- 
matics and morphology, and the CME was found to undergo early 
acceleration, non-radial motion, angular width expansion, and aero- 
dynamic drag in the solar wind as it propagated towards Earth. This 
study and its conclusions are of vital importance to the fields of space 
weather monitoring and forecasting. 



*This is the online version of this thesis which has reduced quality 
images due to file size concerns. 
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Berghmans (2004) E7J 



xvi 



LIST OF FIGURES 



3.3 Example of the SEEDS CME detection and height determination, 



reproduced from Olmedo et al. (2008). (a) shows the running- 



difference image unwrapped into Sun-centre polar coordinates, show- 
ing a CME observed by LASCO/C2 on 12 September 2002. The 
black line distribution across the image represents the positive 
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3.4 An example of how the synoptic maps are generated for the ARTEMIS 



catalogue, reproduced from Boursier et al. (2009b). At a chosen 
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from the high-pass filtering at one scale of the multiscale decom- 
position (called the rows and columns respectively). Bottom left, 
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the angle information (0 - 360°) taken from the gradient space, for 



a CME observed in LASCO/C2 on 1 April 2004 ( |Byrne et ggjQ9j ). M 
3.6 The vectors plotted represent the magnitude and angle determined 
from the gradient space of the high-pass filtering at a particular 
scale. The CME of 2004 April 1 shown here is highlighted very 
effectively by this method (Byrne et al, 2009[ ) 



3.7 A normalised, background subtracted, LASCO/C2 image (left) of a 
CME on 1 April 2004, and the resulting NRGF image (right). The 
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3.8 A chosen scale of the decomposed NRGF image provides a mag- 
nitude image of the edge strengths displayed on the left, which 
is thresholded at one standard deviation from the mean intensity 
to obtain contoured regions of interest that could contain a CME 
(sample contours indicated in red). As shown, the streamers have 
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olded regions of strongest magnitude and angular distribution at 
four scales of the decomposition. The 3D representation on the 
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3.11 The NRGF (left) and resulting detection masks (middle and right) 
for different frames of a CME observed by LASCO/C2 at times 
of 00:00 UT (top), 00:40 UT (middle) and 01:20 UT (bottom) 
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3.12 Example of the difficulty in detecting the faint CME observed by 
LASCO on 23 April 2001. If the CME is far-sided and/or of low 
intensity it becomes difficult to threshold its edges in the image 
compared to the coronal streamers, (a) is the pre-processed CME 
image, (b) is the NRGF image, (c) and (d) show the magnitude 
and angular information from the multiscale decomposition with 
intensity contours overlaid in red. (e) and (f) show the resulting 
CME detection mask in 2D and 3D respectively 11081 
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tion of a CME observed by LASCO/C2 on 24 January 2007. The 
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cor2 cm 

4.3 Kinematic and morphological profiles for the ellipse fit to the mul- 
tiscale edge detection of the 2 January 2000 CME observed by 
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xix 



LIST OF FIGURES 
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scale edge detection of the 18 April 2000 CME observed by LASCO/C2 
and C3. The plots from top to bottom are height, velocity, angular 
width, and ellipse tilt. The CDAW heights are over-plotted with 

a dashed line. The height and velocity fits are based upon the 
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width, and ellipse tilt. The CDAW heights are over-plotted with 

a dashed line. The height and velocity fits are based upon the 
constant acceleration model 11291 

4.7 Kinematic and morphological profiles for the ellipse fit to the multi- 
scale edge detection of the 21 April 2002 CME observed by LASCO/C2 
and C3. The plots from top to bottom are height, velocity, angular 
width, and ellipse tilt. The CDAW heights are over-plotted with 

a dashed line. The height and velocity fits are based upon the 
constant acceleration model 11311 

4.8 Kinematic and morphological profiles for the ellipse fit to the multi- 
scale edge detection of the 1 April 2004 CME observed by LASCO/C2 
and C3. The plots from top to bottom are height, velocity, angular 
width, and ellipse tilt. The CDAW heights are over-plotted with 

a dashed line. The height and velocity fits are based upon the 
constant acceleration model 11331 
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4.9 Kinematic and morphological profiles for the ellipse fit to the mul- 
tiscale edge detection of the 8 October 2007 CME observed by 
SECCHI/COR1 and COR2 on STEREO-B. The plots from top to 
bottom are height, velocity, angular width, and ellipse tilt. The 
height and velocity fits are based upon the constant acceleration 
model ECU 

4.10 Kinematic and morphological profiles for the ellipse fit to the mul- 
tiscale edge detection of the 16 November 2007 CME observed by 
SECCHI/COR1 and COR2 on STEREO-A. The plots from top to 
bottom are height, velocity, angular width, and ellipse tilt. The 
height and velocity fits are based upon the constant acceleration 
model 11361 

5.1 Schematic of the epipolar geometry used to relate the observa- 



tions from the two STEREO spacecraft, reproduced from Inhester 



(2006). This geometry enables us to localise features in 3D space 

by the triangulating lines-of-sight across epipolar planes 11431 

5.2 Schematic of how tie-pointing a curved surface within an epipolar 
geometry is limited in its ability to resolve the true feature, since 
lines-of-sight will be tangent to different edges of the surface and 



not necessarily intersect upon it. Reproduced from Inhester ( 2006 ) . [T441 



5.3 Schematic of the Graduated Cylindrical Shell (GCS) model of a 



flux rope CME, reproduced from Thernisien et al. (2009). Indi- 
cated are the model parameters of front height hf ront , leg height 
h, angle between the legs 2a, cross-sectional radius a, and distance 
from Sun centre O to a point on the edge of the shell r. Two views 
of the GCS model are shown; (a) 'face-on', and (b) 'end-on'. The 
positional parameters of longitude 0, latitude 9, and orientation 
angle 7 are illustrated in (c) 11461 
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5.4 The elliptical tie-pointing technique developed to reconstruct the 
3D CME front, shown here for the 12 December 2008 event. One of 
any number of epipolar planes will intersect the ellipse characteri- 
sation of the CME at two points in each image from STEREO-A 
and B. (a) illustrates how the resulting four sight-lines intersect in 
3D space to define a quadrilateral that constrains the CME front 
in that plane. Inscribing an ellipse within the quadrilateral such 
that it is tangent to each sight-line provides a slice through the 
CME that matches the observations from each spacecraft, (b) il- 
lustrates how a full reconstruction is achieved by stacking multiple 
ellipses from the epipolar slices to create a model CME front that 
is an optimum reconstruction of the true CME front, (c) illustrates 
how this is repeated for every frame of the eruption to build the 
reconstruction as a function of time and view the changes to the 
CME front as it propagates in 3D. While the ellipse characterisa- 
tion applies to both the leading edges and, when observable, the 
flanks of the CME, only the outermost part of the reconstructed 
front is shown here for clarity 11481 

5.5 An ellipse inscribed within a convex quadrilateral. An isometry of 
the plane is applied such that the quadrilateral has vertices (0, 0), 
(A,B), (0, C), (s, t). The ellipse has centre (h,k), semimajor axis 
a, semiminor axis b, tilt angle 5, and is tangent to each side of the 
quadrilateral 11491 

5.6 The back-projection of a STEREO 3D CME front reconstruction 
onto the SOHO/LASCO plane-of-sky, from the observations by 
STEREO-A (red) and STEREO-B (green) at 16:22 UT on 26 April 
2008 fT52l 
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5.7 The back-projection of a STEREO 3D CME front reconstruction 
onto the SOHO/LASCO plane-of-sky, at 16:22 UT on 26 April 
2008. The left panel shows the multiscale edge detection and ellipse 
characterisation of the CME front observed off the east limb from 
STEREO-A. The right panel shows the running-difference point- 
and-click ellipse characterisation of the CME front observed as 
a halo event from STEREO-B. . The middle panel shows the 
LASCO/C2 observation of the CME at 16:30 UT (8 minutes later 
than the STEREO observations) with the back-projected 3D CME 
front reconstruction overplotted for comparison 11541 

5.8 Composite of STEREO-Ahead and Behind images from EUVI, 



COR1, COR2, and HI1 (Byrne et al, 2010). (a) indicates the 



STEREO spacecraft locations, separated by an angle of 86.7° at 
the time of the event, (b) shows the prominence eruption observed 
in EUVI-B off the NW limb from approximately 03:00 UT which 
is considered to be the inner material of the CME. The multi- 
scale edge detection and corresponding ellipse characterisation are 
overplotted in COR1. (c) shows that the CME is Earth-directed, 
being observed off the east limb in STEREO-A and the west limb 
in STEREO-B USSl 

5.9 STEREO-Ahead and Behind COR1 images at 07:35 UT on the 12 
December 2008. Overplotted are: the multiscale edge detections of 
the CME front (red); the ellipse characterisations (blue); and the 
resulting 3D reconstructions back-projected onto the plane-of-sky 

(white) nsa 

5.10 STEREO-Ahead and Behind COR2 images at 14:52 UT on the 12 
December 2008. Overplotted are: the multiscale edge detections of 
the CME front (red); the ellipse characterisations (blue); and the 
resulting 3D reconstructions back-projected onto the plane-of-sky 
(white) [1571 
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5.11 STEREO- Ahead and Behind HI1 images at 01:29 UT on the 13 
December 2008. Overplotted are: the running difference edge 
detections of the CME front (red); the ellipse characterisations 
(blue); and the resulting 3D reconstructions back-projected onto 

the plane-of-sky (white) 11581 

5.12 The 3D CME front reconstruction as it propagates along the Sun- 
Earth line through interplanetary space. This particular frame of 
the reconstruction is from the observations of COR1 at 07:35 UT 

on 12 December 2008 11551 

5.13 The 3D CME front reconstruction as it propagates along the Sun- 
Earth line through interplanetary space. This particular frame of 
the reconstruction is from the observations of COR2 at 14:52 UT 

on 12 December 2008 fTHOl 

5.14 The 3D CME front reconstruction as it propagates along the Sun- 
Earth line through interplanetary space. This particular frame of 
the reconstruction is from the observations of HI1 at 01:29 UT on 

13 December 2008 [T6D 

5.15 The graphical user interface for the 3D visualisation of the CME 
front reconstruction, which may be used to change the orientation 
of the observer, or the location of the definable parameters such as 
time-stamps and distance measures (as displayed on the right). . . 11621 

5.16 The 3D CME front reconstruction from the COR2 Ahead and Be- 
hind frames at 14:52 UT on 12 December 2008. The lines drawn 
from Sun-centre indicate the 'Midpoint of Front' (solid blue), the 
'Northern / Southern Flanks' (solid red /brown), and the 'Mid- 
top / Midbottom of Front' at angles in between (dashed red / brown). 
By taking these measurements across all frames we may determine 
the kinematics and morphology of the CME as plotted in Fig- 
ures KM and KM [T641 
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5.17 The error trapezoid on the tie-pointing of two lines-of-sight in 3D 



space, reproduced from Inhester (2006). The error w in localising 
a point on each plane-of-sky results in a trapezoid with a diag- 
onal measuring wj sin(a/2) as shown, where a is the spacecraft 
separation angle 11651 

5.18 The kinematics of the prominence and 3D reconstructed CME 
front of 12 December 2008. The prominence is observed as the in- 
ner material of the CME, with both undergoing acceleration from 
~ 06:00 -07:00 UT, peaking at ~40 m s" 2 and ~94 m s" 2 re- 
spectively, before reducing to scatter about zero. Measurement 
uncertainties are indicated by one standard deviation error bars. . 11681 

5.19 Kinematic and morphological properties of the 3D reconstruction 
of the 12 December 2008 CME front. The top panel shows the 
velocity of the middle of the CME front with corresponding drag 
model and, inset, the early acceleration peak. Measurement uncer- 
tainties are indicated by one standard deviation error-bars. The 
middle panel shows the declinations from the ecliptic (0°) of an 
angular spread across the front between the CME flanks with a 
power-law fit indicative of non-radial propagation. The bottom 
panel shows the angular width of the CME with a power-law ex- 
pansion 11691 

5.20 The 3D CME front parameters are used as initial conditions for an 



ENLIL with Cone Model MHD simulation flXie et al\ |2004j) and 
the output density (top) and velocity (bottom) profiles of the inner 
heliosphere are illustrated here for the time-stamp of 06:00 UT 
on 14 December 2008. Beyond distances of ~ 50 R the CME is 
slowed by its interaction with the upstream, slow-speed, solar wind 
flow along its trajectory towards Earth, and this accounts for its 
arrival time as detected in-situ by the ACE and WIND spacecraft 
at the LI point near Earth 11741 
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5.21 The in-situ solar wind plasma and magnetic field data observed 
by the WIND spacecraft. From top to bottom the panels show 
proton density, bulk flow speed, proton temperature, and magnetic 
field strength and components. The red dashed lines indicate the 
arrival time of the density enhancement predicted from our ENLIL 
with Cone Model run providing 08:09 UT on 16 December 2008, 
with a potential offset error between our reconstruction and the 
derived model height profiles accounting for an arrival time up to 
13:20 UT. We observe a magnetic cloud signature behind the front, 

as highlighted with blue dash-dotted lines 11751 

5.22 Top: Ellipses characterising the 3D CME front reconstruction in 
the out-of-ecliptic plane along the Sun-Earth line. Bottom: The 
ellipse tilt angle is indicative of the initial effects of 'pancaking', as 
the curvature of the CME front decreases with increasing height 

due to its changing morphology in the solar wind 11781 
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6.1 An example interface of the automated multiscale detection algo- 
rithm based upon thresholding the magnitude and angular infor- 



mation, discussed in Section 3.3 Left: The algorithm can be made 
to cycle through different scales, and threshold any number of con- 
toured regions of magnitude (edge strength). Right: Contours that 
contain a wide distribution of angles signal a CME detection. The 
angular information is normalised to 1 and folded into a 0-180° 
range due to symmetry of the edge normals. A threshold on the 
normalised angular distribution is specified to flag regions as CMEs 
or otherwise (e.g., > 20%). A detection mask is then built through 
multiple scales. The limitations to be overcome for a robust au- 
tomation of this technique include developing dynamic thresholds 
such that multiple contours of CME edges are not fragmented, and 
increased angular distributions due to streamer deflections are not 
mistakenly labelled as CMEs (a non-trivial task). Moreover, an 
automated pixel chaining of the CME edges must be included in 
order to produce an ellipse characterisation of the CME front in 
the image 11891 

6.2 An example interface which could be developed for the potentially 
automated multiscale filtering and ellipse tracking of CMEs ob- 
served by STEREO-A and B. Images from the respective instru- 
ments appear on the left, and the resulting parameters from the 
characterisation of the CME front appear on the right. Any man- 
ner of information may be chosen for display, and clocked through 
time as the CME progresses 11901 
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6.3 The magnitude (edge strength) from the multiscale filtering of a 
CME observation, unwrapped into polar coordinates (r, 9) with 
axis units of pixels and degrees, respectively. A scanning-line runs 
over the image to produce the scatter of edge normals at the bot- 
tom of the image, where each normal is plotted with a vector mag- 
nitude and associated angular information from the multiscale fil- 



ter (cf. the edge normals of Figure 3.6). An end-projection of 
the normals along the scanning-line is plotted in the bottom right 
of the image. The scanning-line is located at angle 328°, passing 
over a CME in the image. The resulting edge normals show a slice 
through the angular distribution of the CME, and as the scanning- 
line moves along the image, the end-projection shows a continuous 
rotation of angles due to the curvature of the CME front. Detect- 
ing this continuous rotation of angles may be used to distinguish 
CMEs from streamers which show an abrupt, stepwise change in 
angle as the scanning-line crosses over them. This may alleviate a 
need for the intensity thresholding and contouring which can fail 
to detect faint CMEs or portions thereof 11921 



6.4 An example ridgelet, reproduced from Gallagher et al. (2010). The 
first graph shows a typical ridgelet, and the second to fourth graphs 
are obtained from simple geometrical manipulations, namely rota- 
tion, rescaling and shifting 11941 

6.5 The curvelet filtered image of a CME observed by LASCO/C2 on 



18 April 2000, reproduced from Gallagher et al. (2010). The detail 
along the curved CME front is enhanced as a result of the curvelet 
technique following the removal of certain coefficients probably due 
to noise 11951 
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6.6 A schematic of the possible 3D magnetic field topology of a CME, 
as inferred from typical observations of a leading edge, dark void 
and bright core. There may also be an associated on-disk post-flare 
arcade at the base of the postulated current sheet. 3D reconstruc- 
tions of CMEs are important for testing the validity of theoretical 
models. 

Image credit: http://www.nrl.navy.mil 11971 

6.7 'Top' view of the polarimetric reconstruction method applied to 
both STEREO-A and B observations of a CME at 01:30 UT on the 
31 December 2007 (rotated to the Sun-Earth frame), reproduced 
from Moran et al. (2010). The red points are the reconstruction 
from STEREO-A, and the blue points from STEREO-B, with the 
green points showing the regions of overlap 11991 

6.8 A theoretical model for a CME with constant acceleration 2 m s~ 2 
and initial velocity 300 km s -1 , and two simulations of how the 
resulting profiles for a noisy sample of datapoints behave using 
3-point Lagrangian interpolation 12011 

6.9 The resulting velocity profile for the 3D reconstructed CME front 



of the 12 December 2008 using the inversion technique of Kontar 



fc MacKinnon (2005) [2021 
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Chapter 1 
Introduction 



The Sun, as provider of light and heat to all life on Earth, has been a con- 
stant source of mystery and wonder to humankind. History recounts numerous 
tales inspired by our connection with the Sun: from its worship as a deity in 
the earliest civilisations, to the appreciation of its seasonal influence marked by 
structures like Newgrange, and the eventual observance of its complex behaviour 
with the development of telescopes and scientific intrigue. As our nearest star, 
astronomers have increasingly taken interest in the complexities of the Sun, and 
now in the modern age of space exploration numerous observatories have been 
built specifically to monitor solar activity and further our understanding of its 
dynamic behaviour. 
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1.1 The Solar Interior 



The Sun is a G2V main sequence star of luminosity L = 3.85 x 10 W, mass 



M = 1.99 x 10 30 kg and radius R Q = 6.96 x 10 8 m (Prialnik, 2000). It was born 



from the gravitational collapse of a molecular cloud approximately 4.6xl0 9 years 
ago, is currently in a state of hydrostatic equilibrium (VP = —pg), and predicted 
to enter a red giant phase in another ~ 5 billion years before ending its life as 



a white dwarf (Phillips, 1999). Since we cannot directly observe the interior 



of the Sun, its structure and evolution are fundamentally realised with the use 



of the 'standard solar model' (SSM; Bahcall, 1989), which is a mathematical 



treatment of stellar structure described by several differential equations derived 
from basic physical principles. The SSM is constrained by the well-determined 
boundary conditions of the Sun's luminosity, radius, age and composition, and 
thus provides a basis for understanding the mechanisms of energy transport in 
the solar interior. It assumes hydrostatic equilibrium, with energy generated by 
nuclear fusion, although small effects of contraction or expansion are included, 
and any abundance changes are caused solely by the nuclear reactions. The 
SSM is the end product of an iterative process that converges on an optimum 
description of the internal energy generation and transport, and overall evolution 
of the Sun. 
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Figure 1.1: Illustration of the structure of the Sun. The core is the source of 
energy, where fusion heats the plasma to ~ 15 MK. Energy is transported from the 
core by radiative processes in the radiation zone. The convection zone is heated 
from the base at the tachocline, allowing convective currents to flow to the photo- 
sphere. Locations of strong magnetic fields inhibit convection and appear as dark 
sunspots on the photosphere. These strong magnetic fields extend into the upper 
atmosphere of the Sun, responsible for coronal loops, prominences and streamers. 
Image credit: eu.spaceref.com. 
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The fundamental energy process driving the Sun is nuclear fusion in the core, 
through the proton-proton chain at temperatures of ~ 15 MK: 

JH + JH -> lR + e + + v e (1.1) 
lE + \E -> pie + 7 (1.2) 
pie + pie -> pie + 2 }H (1.3) 

where JH is a proton, is the deuteron isotope of hydrogen, |He and gHe are 
helium isotopes with 1 and 2 neutrons respectively, e + a positron, v e an electron- 
neutrino and 7 a gamma ray. The resulting energy release for one complete 



reaction chain is approximately 4.3 x 10 J (Phillips, 1995). The core extends 



from the centre out to ~ 0.25 R , followed by the radiation zone out to ~ 0.75 R, 



0) 



then the convection zone out to the solar surface at 1 R (Figure 1.1). The 
temperature across the radiation zone drops to ~ 5 MK with radiation being 
the most efficient method of energy transport. This radiation field is closely 
approximated by a black body, for which the spectral radiance is described by 
the Planck equation: 

2/ic 2 /i 2 

5a(T) = A 5 [exp (hc/XkT) - 1] ( 1A > 

where B\(T) is the intensity of radiation per unit wavelength interval (at tem- 
perature T), h is the Plank constant, c is the speed of light, fi is the refrac- 
tive index of the medium, and k is the Boltzmann constant. By Wien's law 
XmaxT = 2.8979 x 10~ 3 m K we determine that the radiation is in the form of 
X-rays, and these high-energy photons undergo random walks in the plasma, es- 
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caping into the convection zone on time-scales of 10 6 years. The optically thick 
convection zone then transports energy by fluid motion across the temperature 
gradient between its base (~ 1-2 MK) and the solar surface (~ 5,800 K). Plasma 
elements move sufficiently rapidly for the energy interchange with their surround- 
ings to be negligible, i.e., they change adiabatically. A useful measure of when 
convection is likely to occur is given by the Schwarzschild criterion: 



dlogT 



dlogP 



> izi (1.5) 

star 



where 7 = Cp/Cy is the ratio of specific heats, equal to 5/3 for a perfect 
monatomic gas. Essentially convection occurs once the absolute magnitude of 
the radiative gradient becomes larger than the absolute magnitude of the adi- 
abatic gradient, so that rising elements of plasma remain buoyant and move 
towards the surface before they can lose heat to their surroundings. The rising 
and falling parcels of plasma create the granulation effects observed on the sur- 
face, with granules ranging in size from hundreds to thousands of kilometres and 
dissipating over tens of minutes. (Details of the above radiative and convective 
processes are found in, e.g., |Kitchin| ( |l987| >; [Zirin] ( p98| )). 



Between the radiation and convection zones is a relatively thin interface called 
the tachocline, where the solid body rotation of the radiative interior meets the 
differentially rotating outer convection zone. It thus has a very large shear profile 
which could account for the formation of large scale magnetic fields in the solar 
dynamo. The magnetic field of the Sun has an overall dipolar configuration, with 
opposite polarities dominant at each pole. The differential rotation of the Sun's 
convection zone causes a large-scale winding up of the magnetic field, named the 
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f2-effect, while the effects of the coriolis force and smaller scale motions of the 



plasma can give twist and writhe to the field, named the a-effect (Figure 1.2). 
Buoyancy effects cause the magnetic field to rise up through the convection zone 
and protrude through the surface of the Sun, observed as sunspots on-disk mark- 
ing the footpoints of over-arching concentrations of magnetic flux extending up 
through the solar atmosphere. In a given hemisphere one magnetic polarity leads 
the sunspot group and the opposite follows, while in the opposite hemisphere 
the situation is reversed (Hale's law). The tilt angle between leading and trail- 
ing sunspots has an average value of 5.6° relative to the E-W line (Joy's law). 
Furthermore, sunspots are observed to migrate from high latitudes towards the 
equator over an 11 year cycle due to the continual build-up of field by the afl- 
effect (Sporer's law). These combined effects lead to an increase of oppositely 
oriented poloidal field at each of the poles, neutralising the field there and re- 
sulting in the magnetic dipole flipping every 11 years. This 22 year periodicity 
is known as the solar cycle, and gives rise to periods of increased and decreased 
solar activity manifested by the frequency of phenomena such as active regions, 



flares and transients in the solar atmosphere (Schrijver & Zwaan, 2000). 
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Figure 1.2: Illustration of the aVt effect of winding-up magnetic field due to the 



differential rotation of the Sun, reproduced from Babcock (1961). Sunspots visible 
on the disk are as a result of protruding field with positive p and negative / polarity 
as shown. 



1.2 The Solar Atmosphere 

The Sun's atmosphere is composed of all regions extending from the photosphere 
out into the heliosphere. It may be separated into distinct regimes dependent 
on the density and temperature profiles. These are plotted in Figure 1.3| for a 



ID static model of the solar atmosphere. The layers are generally stratified into 
photosphere, chromosphere, transition region, and corona; having a decreasing 
density with increasing height from the photosphere, but from the chromosphere 
up the temperature increases with a dramatic jump in the transition region giving 
rise to the so-called 'coronal heating problem'. However this stratification is a 
simplified view and the solar atmosphere in reality is an inhomogenous mix of 
photospheric, chromospheric and coronal zones due to complex dynamic processes 
such as heated upflows, cooling downflows, intermittent heating, nonthermal elec- 
tron beams, field line motions and reconnections, emission from hot plasma, ab- 
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sorption and scattering in cool plasma, acoustic waves, and shocks (Aschwanden 



2005). The interplay between the magnetic and gas pressure represents an im- 
portant determining factor in the behaviour of structures throughout the solar 
atmosphere, quantified by the plasma-/? term: 



B = = UkT (1.6) 

Pmag (i^/M 1 ) 



This is illustrated in Figure L4 for the different layers of the atmosphere. At 
the photospheric level the plasma-/? is large, and plasma motions dominate over 
the magnetic field forces. Through the chromosphere and corona the plasma-/? 
decreases to low values where the magnetic field structures are seen to suspend 
plasma in loops and filaments. Finally in the extended upper atmosphere the 
plasma-/? rises again, and the magnetic field is advected out with the solar wind 
plasma flow to ultimately form the Parker spiral. 



1.2.1 Photosphere 

The surface of the Sun is the photosphere defined as the point where the optical 
depth equals 2/3 for wavelengths of visible light, centred on 5,000 A (r 50 oo ~ 2/3 
for I/Iq = e~ T ). The spectrum of light emitted has a profile like that of a black 
body with an effective temperature of 5,800 K interspersed with the Fraunhofer 
absorption lines due to the tenuous layers above the photosphere. It has a particle 
density of ~ 10 23 m~ 3 and a thickness of less than 500 km. Cooler regions called 
sunspots have temperatures of 4,000-4,500 K and are due to intense magnetic 
field activity that acts to suppress convective plasma motion. Granulation of the 
photosphere is observed as the manifestation of plasma motion in the convection 
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Figure 1.3: A ID static model of electron density N e [cm -3 ] and temperature 



T e [K] profiles in the solar atmosphere, reproduced from Gabriel & Mason (1982). 
In the chromosphere, the plasma is only partially ionised. The plasma becomes 
fully ionized at the sharp transition from chromospheric to coronal temperatures. 



zone below, with typical cell sizes on the order of 1,000 km in diameter. They 
occur when hot plasma rises to the surface and is transported along it to the 
granule edges, which appear darker as the plasma cools and descends. The gas 
pressure dominates the magnetic pressure (/3 > 1), and the magnetic field is 
effectively coupled to the plasma motion which sweeps it into the inter-granular 
network. 



1.2.2 Chromosphere 

Above the photosphere lies the chromosphere where the temperature initially 
drops to a minimum of ~ 4,500 K before increasing to ~ 20,000 K with increasing 



height from the Sun (Figure 1.3). It is approximately 2,000 km thick and the 
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Figure 1.4: Plasma f3 as a function of height for a regime of magnetic field 
strengths between ~ 100 and ~ 2,500 G, reproduced from Gary (2001). The dotted 
lines segregate the layers of photosphere (/3 > 1), chromosphere and corona (/? < 1), 
and the solar wind (/3 > 1). 



density falls by a factor of almost a million from bottom to top, so the magnetic 
field begins to dominate the chromospheric structure (/3 < 1). The second law of 
thermodynamics does not permit heating of the chromosphere with the thermal 



energy of the cooler photosphere below. Biermann (1948), Schwarzschild (1948) 



and Schatzman (1949) put forward ideas on the acoustic wave heating of the 
chromosphere as a result of the convective plasma motions in the photosphere 
and convection zone beneath. Referred to as the BSS model, the hypothesis is 
that acoustic waves transport energy upward with little dissipation once the ve- 
locity is below the sound speed. As the density drops and the velocity reaches 
the sound speed, the waves steepen into shocks and rapidly dissipate the energy, 
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consequently heating the material (Zirin, 1998). However, acoustic wave heating 
does not apply in regions of strong magnetic field where motions, and therefore 
heating, are suppressed. This has led to work on Alfven wave heating theories, 



first introduced by Osterbrock (1961). An Alfven wave is a type of magnetohy- 
drodynamic wave that propagates in the direction of the magnetic field with the 
magnetic tension providing the restoring force and the ion mass density provid- 
ing the inertia. In Alfven wave heating theories the magnetic field itself is thus 
responsible for depositing energy from the subsurface into the chromosphere and 
above. These theories better sit with observations of vigourous heating above 
plages and emerging flux regions, since they imply the amount of heating is pro- 
portional to the rate of magnetic change. 

While the brightness of the photosphere overwhelms that of the chromosphere 
in the optical continuum, the hotter chromospheric temperatures lead to the hy- 
drogen being ionised, resulting in strong Ha emission. Filaments are observed 
as dark channels on-disk in Ha images (called prominences when seen on the 
limb). Numerous plasma columns called spicules are also observed on the limb, 
that typically reach heights of ~ 3,000 -10,000 km above the Sun's surface and 
are very short-lived (rising and falling over ~ 5 — 15 minutes). 

Between the chromosphere and corona lies the transition region where the 
temperature jumps rapidly to over 1 MK. It is only about 100 km thick and 
it marks the point where magnetic forces dominate completely over gravity, gas 
pressure and fluid motion (/3 <C 1). The extreme temperatures result in prominent 



UV and EUV emission from carbon, oxygen and silicon ions (Mariska, 1992). 
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1.2.3 Corona 



The outermost part of the solar atmosphere is the corona, with electron den- 
sities ranging from ~ (1-2) x 10 14 m~ 3 at its base height of ~ 2,500 km above 



the photosphere, to < 10 12 m 3 for heights >1 R above the photosphere (As- 



chwanden , |2005 ) . The density varies across coronal holes which can have a base 



density of ~ (0.5-1) x 10 14 m 3 , or across streamer regions with higher densi- 
ties of ~ (3-5) x 10 14 m~ 3 . Active regions that suspend and confine plasma in 
strong over-arching magnetic fields usually have the highest coronal densities of 
~2 x 10 14 -2 x 10 15 m -3 . The temperature of the corona is generally > 1 MK, 
as indicated by emission from highly ionised iron lines, for example, which again 
appears to contradict the second law of thermodynamics given the much cooler 
layers of the chromosphere and photosphere below (the 'coronal heating prob- 
lem'). Its temperature structure is far from homogeneous, revealed in images 



such as that of Figure 1.5 from the Solar Dynamics Observatory (SDO). Loop 
structures are observed at temperatures of 2-6 MK across regions of increased 
magnetic field density (such as above active regions/sunspots), and closed field 
regions are observed at temperatures of 1-2 MK across the quiet Sun, while 
open field regions of coronal holes have temperatures < 1 MK. These high tem- 
peratures lead to EUV and X-ray emission due to ionisation and recombination 
processes from the interactions between photons, electrons, atoms and ions. The 
continuum and line emission of the solar corona result from the contributions of 
bound-bound transitions (excitations and de-excitations), bound- free absorption 
(photoionisation), free- free absorption (and its inverse bremsstrahlung emission), 
and electron scattering. It is the latter process, called Thomson scattering, that 
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Figure 1.5: Composite EUV image of the low solar corona recorded by the Atmo- 
spheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) 
on 23 May 2010. 

Image credit: http://sdowww.lmsal.com/ 
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produces the white-light corona as visible during a solar eclipse or with the use of 
a coronagraph to occult the solar disk, which is six orders of magnitude brighter 
in optical wavelengths. 

The corona we observe comprises several parts: 

• The K-corona has a strongly polarised continuous emission spectrum due to 
Thomson scattering of photospheric light by the free electrons of the coronal 
gas, and it dominates within the first few R . It produces a polarised 
white-light continuum without the Fraunhofer lines which are broadened by 
Doppler shifts due to the fast electron motions at such high temperatures. 



The intensity of the K-corona gives the coronal electron density (Koutchmy 



et al. 1991). 



• The F-corona is due to scattering of sunlight by interplanetary dust parti- 
cles, and contains the Fraunhofer lines. It is roughly equal in intensity to 
the K-corona at ~4 R , and dominates at greater distances. 

• The E-corona is due to emission from highly ionized coronal atoms such as 
iron and calcium. 

• The T-corona is caused by thermal (infrared) emission of the interplanetary 
dust. It is an unpolarised continuum, insignificant in the visible part of the 
spectrum. 

In contrast to the chromosphere, solar interior, and indeed the heliosphere, the 
magnetic pressure in the corona dominates over the gas pressure and so governs 
the coronal plasma dynamics (/3 < 1). The coronal structure we observe is thus 
shaped by the magnetic fields of the Sun, resulting in extended polar regions 
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where there is mainly open magnetic field, and 'helmet-streamers' spanning the 
equatorial latitudes where, except for coronal holes, the field is mostly closed. 
Since these features are magnetically governed, the shape of the corona varies 
greatly over the solar activity cycle: it appears rounder at solar maximum, when 
multiple streamers emerge at various latitudes distributed across the Sun; and it 
appears more elliptical at solar minimum, when only a few streamers are present, 
lying closer to the equator. 



Following Chapman & Zirin (1957) the description of a static corona leads 
to an unreasonable pressure value at large distances from the Sun. This is out- 
lined below, beginning with the assumption that the corona is in hydrostatic 
equilibrium: 

dP = _ p GMe (1 7) 

dr r 2 

The plasma density is p = nm p , the pressure contribution from the protons and 
electrons is P = 2nksT, and the coronal heat flux is q = kVT with thermal 
conductivity k = k T 5j/2 . In the absence of heat sources or sinks V • q = so in 
a spherically symmetric system we can write: 

r 2 dr \ dr J 

Applying the boundary condition that the temperature tends to zero at large 
distances from the Sun, we obtain: 



T = T (^) 2/? (1.9) 
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where Tq = 2 MK is the temperature of the low corona at height r$ = 1.05 R 
from Sun centre. This would mean T w 4 x 10 5 K at Earth (1 AU~215 R©), 
close to measured values. Rewriting in terms of pressure and integrating, results 
in: 

which implies that as r — > oo the coronal pressure tends towards a finite constant 
value significantly larger than the pressure of the interstellar medium (ISM); 
P ^> Pism- This means the static coronal model is unphysical, and a dynamic 
model in which the material flows outward from the Sun must be considered, 
leading to a description of the solar wind. 



P{r) = Po ex P 



7 GM Q m p 
5 2kT Q r 



1.2 A Solar Wind 

The solar wind is the constant out-stream of charged particles of plasma from 
the Sun's atmosphere due to the persistent expansion of the solar corona. The 
wind consists mostly of electrons and protons at energies of ~ 1 keV, observed 
in two regimes of propagation: the slow solar wind with speeds of ~400 km s -1 ; 
and the fast solar wind with speeds of ~800 km s -1 , originating from regions of 
open magnetic field such as coronal holes. Thermal velocities of the particles are 
calculated at ~260 km s _1 for coronal temperatures on the order of 3xl0 6 K, 
while the escape velocity in the Sun's gravitational field in the low corona can 
be ~500 km s _1 . The additional energy to accelerate the solar wind is imparted 
by the pressure gradient Ps un Pism to attain the measured solar wind speeds 



(Parker, 1958). The Parker model assumes the outflow is steady, spherically 



symmetric and isothermal. The momentum conservation equation of the corona 
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Figure 1.6: The five classes of Parker's solar wind solution for a steady, spherically 
symmetric, isothermal outflow. 



takes the form: 



pv 



dv 
dr 



dP GM Q 
dr ^ r 2 



(1.11) 



Considering mass conservation m = 4nr 2 pv = constant, we obtain: 



d_ 
Or 



(r 2 pv) = 



1 dp ldv 2 

p dr v dr r 



(1.12) 



So for a perfect gas P = RpT Equation 1.11 can be written 



RT\ dv 2RT GM C . 



v — 



v J dr r 



+ 



(1.13) 
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A critical point occurs when d r v — > so we define: 



GMq 
2vl 



where v c = V RT 



;i.i4) 



and rewrite Equation 1.13 as: 



(v 2 ~ vl) 



1 dv 
v dr 



2^(r-r c ) 



;i.i5) 



Integrating Equation 1.15 gives Parker's 'solar wind solutions' 



-In 



v c 



Hi) [ - ) + 4- + C 

r r r 



;i.i6) 



where C is a constant of integration, leading to five potential solutions as plotted 



in Figure L6} Solutions I and II are double-valued, with II being disconnected 
from the surface. Solution III is too large (supersonic) close to the Sun. So- 
lution IV is called the 'solar breeze' as it remains subsonic. Solution V is the 
standard solar wind solution, although the assumptions of radial expansion and 
isothermality are not completely true in reality, so it is only an approximate 
characterisation of the observed solar wind. Nonetheless it is sufficient to convey 
the dynamic expansion of the corona and ultimate supersonic regime of outflow, 
often described akin to a de Laval nozzle which is used to accelerate flows from 



subsonic to supersonic speeds (as detailed in Goossens, 2003, for example). The 
presence of the solar wind was confirmed in 1959 by the Lunik I probe, and in 



1962 by the Mariner II mission en route to Venus (Neugebauer & Snyder, 1962) 



Since the gas pressure dominates over the magnetic pressure in the solar at- 
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Figure 1.7: Schematic of the Parker Spiral in the heliosphere. The streamlines of 
the solar wind act to drag out the magnetic field lines of the Sun, which become 
wound up in an Archimedean spiral as a result of the Sun's rotation. 
Image credit: Steve Suess, NASA/MSFC. 



mosphere (/3 > 1), the solar wind acts to drag out the magnetic field lines of 
the Sun which become wound up as a result of solar rotation to form the Parker 



Spiral (Figure 1.7). This is an Archimedean spiral drawn by the magnetic field 
lines as they are advected outward by the solar wind, described by the equation: 



r-r = fi(0-9 ) (1.17) 
where 9 is the polar angle, Q = 2.7 x 10~ 6 rad s _1 is the angular rotation rate 



of the Sun, r is the distance, and v is the solar wind speed (Parks, 2004 Zirin 



1998). The different speed streams can also lead to the formation of co-rotating 
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interaction regions (CIRs) where the fast wind encounters the slow wind ahead 
of it in the Parker spiral, and can form shocks in the solar wind. 

The solar wind does not extend infinitely, but eventually terminates when it 
reaches the edge of the heliosphere. The point where the solar wind slows from 
supersonic to subsonic speeds is called the termination shock, the observation of 
which is reported in a series of papers discussing the data from Voyager II as 



it began to cross the shock in August 2007 (Burlaga et al., 2008 Decker et al. 



2008 Richardson et al, 2008). Beyond the termination shock, the wind comes 



into pressure balance with the ISM to form the heliosheath, whose outer boundary 



is called the heliopause (Figure 1.7). In the heliosheath the continually slowing 
wind is compressed and becomes turbulent through its interaction with the ISM 



(Opher et al. 2009). As the heliosphere moves through interstellar space, a bow 
shock is thought to form ahead of the heliopause as it encounters the ISM. 



1.2.5 Space Weather 

Space weather is the name attributed to phenomena involving ambient plasma, 
magnetic field, radiation and other matter affecting the conditions in space as 



part of the vast field of heliophysics (Schrijver & Siscoe, 2010). It is predom- 



inantly due to the influences of the solar wind, flares and transients, and the 
interplanetary magnetic field (IMF). Earth's magnetosphere provides a natural 
shielding of the planet from the solar wind and space radiation, although large 
flares and coronal mass ejections (CMEs) that accelerate solar energetic particles 
(SEPs) can lead to geomagnetic storms at Earth which perturb the magneto- 
sphere and cause increased electric currents in the ionosphere. Auroral sightings 
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Figure 1.8: The Aurora Borealis, or Northern Lights, photographed above Bear 
Lake, Eielson Air Force Base, Alaska. Aurorae result from photon emissions of 
excited oxygen and nitrogen atoms in the Earth's upper atmosphere, as a result of 
geomagnetic storms and space weather. 
Image credit: Joshua Strang via Wikimedia Commons. 



throughout history are an indicator of geomagnetic storm occurrences (e.g., Fig- 



ure 



1.8), and by 1837 they were realised to be caused by electric currents in 



the upper atmosphere (Olmstead, 1837). Sabine (1852) went on to show that 



there was a detailed correlation between the sunspot cycle and the frequency 
of auroral displays, implicating solar activity as the ultimate cause of the auro- 
ral phenomenon. A key event in the realisation of how strongly solar activity 
can influence us on Earth, was the Carrington-Hodgson flare that occurred on 



1 September 1859 (Carrington, 1859), causing widespread sightings of aurorae 



down to latitudes as low as ~ 18° and the loss of a significant portion of the 



telegraph service for many hours (Green et al. , 2006). One model for the event 
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comprises the ejection of two CMEs from the Sun on 27 August and 1 September 
1859, whose interaction as the second CME ploughs through the first produced 
a shock responsible for the extreme nature of the 2-3 September auroral event 



(Green & Boardsen, 2006). More recently, a severe geomagnetic storm on 13 
March 1989 caused the collapse of the Hydro-Quebec power network due to a 
transformer failure from the geomagnetically induced currents (GIC). Six million 
people were left without power for nine hours, with a substantial economic loss. 
The cause was a CME ejected from the Sun on 9 March 1989 impacting the Earth 
several days later. Other CMEs have similarly knocked out communication satel- 
lites such as the Canadian Aniks El and E2 and the international Intelsat K on 
20 January 1994, and the AT&T Telstar 401 on 7 January 1997. Moreover, the 
increased radiation of solar and space weather storms poses a risk to astronauts 
and high-altitude flight passengers, particularly when travelling over the poles. In 
the modern era of technological advancement and increased dependency on satel- 
lites communications, GPS networks and power distribution grids, the influence 
of space weather is an increasing cause for concern. A recent report by the US 



National Research Council (Committee On The Societal & Economic Impacts Of 



Severe Space Weather Events^ 2008) indicates that the potential economic cost of 



a high-level geomagnetic storm could be up to $2 trillion. Thus the monitoring 
and forecasting of potentially hazardous events is of great importance to society 
at large, with particular emphasis on CMEs and the dynamics governing their 
propagation through space. 
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1.3 Coronal Mass Ejections 



"We define a coronal mass ejection to be an observable change in coronal structure 
that (1) occurs on a time scale between a few minutes and several hours and (2) 
involves the appearance [and outward motion] of a new, discrete, bright, white- 
light feature in the coronagraph field of view. " 



(Hundhausen et al, 1984) 



Through the association of early observed flares, notably the Carrington event 
of 1859, and the detections of geomagnetic storms at Earth, the theory was put 
forward that plasma transients may be ejected from the Sun and possibly impact 



the Earth's magnetic field several days later (e.g., Lindemann, 1919). Observa- 
tions of prominence disappearances provided evidence that the material may be 
rising through the corona with increasing velocity to eventually exceed the escape 



velocity of the Sun and erupt in to space ( Kiepenheuer 1953). Theories were de- 
veloped by scientists who postulated that the magnetic field should be affected by 
these ejections, as they might act to drag out field lines or sever completely from 
the Sun through magnetic reconnection, and potentially drive shock disturbances 



(Gold, 


1962 


Piddington 


1958 



advent of space-borne coronagraphs that these transients were realised to be a 
common occurrence on the Sun, and their potential geomagnetic effects on Earth 
led them to become a topic of great interest. An example of such observations is 



shown in Figure 1.9 



These coronal mass ejections (CMEs), as they became known, are the largest 
manifestation of the shedding of solar magnetic field during the Sun's 22 year 
cycle. Every 11 years the magnetic axis of the Sun flips, giving rise to periodic 
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Figure 1.9: Observation of a CME and prominence lift-off from the EUVI and 
COR1 instruments of the SECCHI suite on board the STEREO-A spacecraft. The 
field-of-view extends to ~ 4 Rq . The complexity of the magnetic field driving the 
eruption is clearly indicated by the twisted geometry of the bright ejecta. 
Image credit: http://stereo.gsfc.nasa.gov/ 
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patterns in the activity called solar minimum and maximum. At solar minimum 
a CME may occur up to once a week but during solar maximum they can be as 
frequent as three a day. CMEs that travel directly towards or away from the ob- 
server are seen to encircle the full disk of the Sun and are thus named halo CMEs. 
In the case of a halo CME coming toward the Earth, the particle densities and 
energies involved can cause geomagnetic storms, especially if the orientation of 
the CME's impacting magnetic field is oppositely directed to that of Earth's mag- 
netosphere since it can open up the field lines and penetrate deeper into Earth's 
atmosphere. This is referred to as space weather, and understanding this interac- 
tion is of considerable practical importance because technological systems, such 
as communications and navigation satellites, can suffer interruptions or damage. 
To this end, missions such as the Solar and Heliospheric Observatory (SOHO; 



Domingo et al, 1995) and Solar Terrestrial Relations Observatory (Kaiser et al. 



2008 ) have been launched to study CMEs. 



CMEs are observed as a typical three-part plasma structure of a bright leading 



front, dark cavity, and bright core (Illing & Hundhausen, 1985). This configura- 
tion is indicative of a magnetic loop system erupting off the Sun as a flux rope (a 
helically twisted bundle of field lines) or magnetic bubble with plasma embedded 
within and coronal material being swept up ahead of it. Thus the forces acting 
on CMEs are described within the context of magnetohydrodynamics (MHD) as 
outlined below. 
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1.3.1 Magnetohydrodynamic Theory 

The interplay between the plasma and magnetic fields of the Sun, notably in 
phenomena such as flares and CMEs, may be described through the coupling 
of the equations of electromagnetism with the theory of fluid motions. MHD 
attempts to combine Maxwell's equations with the fluid equations through the 
relative dependence on the electron motion in the currents set up in the plasma 
and the effects of the magnetic fields. Thus we obtain the induction equation for 
magnetised plasma and describe how the field may undergo non-ideal (resistive) 
MHD processes such as magnetic reconnection - an important basis of many CME 
models. 



1.3.1.1 Maxwell's Equations 

Maxwell's equations describe the interaction of magnetic field B and electric field 
E according to: 



1 <9E 

VxB = /io j + -— (1.18) 
c l at 

V-B = (1.19) 
<9B 

VxE = -— (1.20) 
V-E = -p (1.21) 



where j is the current density, p is the charge density, p is the magnetic per- 
meability of a vacuum, eo is the permittivity of free space, and c is the speed of 



light. The second term of Ampere's law (Equation 1.18) may be neglected if the 
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typical plasma velocities are much less than the speed of light: 

V x B = p j (1-22) 

1.3.1.2 Fluid Equations 

The mass continuity equation states that matter is neither created nor destroyed: 

^ + V-pv = (1.23) 

where p is the plasma density, and v the plasma velocity This can be expanded 
to give: 

|^ + (v-V)p + pV-v = (1.24) 

where for an incompressible fluid the convective time derivative is zero. (This is 
the derivative D/Dt = d/dt + v ■ V taken along a path moving with velocity v.) 
So the mass continuity equation reduces to: 

V-v = (1.25) 

The equation of motion (F = ma) for a CME may be written: 

P^ = -Vp + j x B + pg- l -pv 2 A cme C D (1.26) 

where p is the pressure, j x B is the Lorentz force, g is gravity, and the drag force 
depends on the cross-sectional area A cme and drag coefficient Co- We neglect 
viscous forces. 

In addition, Ohm's law couples the plasma velocity to the electromagnetic 
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fields by: 



j = a(E + v x B) 



1.27) 



where a is the electrical conductivity. 



1.3.1.3 The Induction Equation 



It is possible to eliminate the electric field E by combining Ampere's law (Equa- 



tion 1.18) and Ohm's law (Equation 1.27): 



E = -v x B + —V x B 



and substituting into Faraday's law (Equation 1.20) to obtain 



dB 

~dt 



V x (v x B) - V x (r?V x B) 

V x (v x B) - 77V x (V x B) 

V x (v x B) + rj [V 2 B - V (V • B)] 



1.28) 



(1.29) 
(1.30) 
(1.31) 



where 7] = l/^a is the magnetic diffusivity. Using the solenoidal constraint 
(Equation |1.19[ ) provides the induction equation: 



O 



dB 

~dt 
B 

1 



V x (v x B) + r]V 2 B 



(1.32) 



O 



vB 



O 



This equation forms the basis of any model that considers magnetised plasma 
motion on a variety of length scales, e.g., from magnetic confinement devices on 
Earth, to the dynamo action of the Sun's magnetic field. How a magnetic field 
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topology will respond to the forces of plasma motion, and vice versa, is governed 
by the ratio of the terms in the induction equation. 

1.3.1.4 The Magnetic Reynolds Number 

The magnetic Reynolds number is the ratio of the advection and diffusion terms 
in the induction equation: 



R m = VX ' V 2 ^ B) « 5* (1.33) 



for plasma speed v and length scale l . 

If R m ^> 1 the induction equation is approximated by: 



— = V x (v x B) (1.34) 



and the coupling of the magnetic field to the plasma motion is strong, so the 
topology of the field changes on the timescales of the plasma motion: 



B v l ( ^ 

=?■ Tmotion ~ \l.60) 



^motion ^0 ^0 

This is known as the 'frozen-in' condition, whereby field lines are carried with the 
plasma motion. For example, in the corona the length scales are Iq ~ 1,000 km, 
velocities are t>o ~ 1,000 m s" 1 , and the magnetic diffusivity r\ ~ 1 m 2 s _1 , resulting 
in a magnetic Reynold's number of Rm ~ 10 9 . 
If R m 1 the induction equation is approximated by: 



f = ^V 2 B (1.36) 
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and the magnetic field can diffuse through the plasma and change its topology: 

B Bg Zq 2 

~ ^7~~2 ^ 7~dif fusion ~ 
7~dif fusion ^0 V 

Magnetic diffusion is an important condition for magnetic reconnection to occur 
such as in a current sheet where the length scales are small enough (on the order 
of metres) to account for the observed restructuring of magnetic field (which is 
on the order of seconds in flaring active regions for example). 



1.3.1.5 Magnetic Reconnection 

The large energies released in flares and CMEs and observed restructuring of the 
coronal magnetic field are often attributed to the phenomenon of magnetic re- 
connection occurring on the Sun. Magnetic reconnection is generally defined as 
a change in connectivity of field lines in time. In an ideal plasma the 'frozen-in' 
condition is met and the magnetic field is coupled to the plasma motion. How- 
ever, when regions of opposite polarity flux come together, a boundary layer will 
form to separate the two regimes of magnetic field in a form of pressure balance. 
The high resistivity of such a system counteracts the currents within and allows 
the occurrence of non-ideal MHD processes and the formation of structures hav- 
ing small spatial scales, e.g., a thin current sheet. This leads to a low magnetic 
Reynolds number and allows diffusion to occur (illustrated by the shaded regions 



in Figure 1.10). The connectivity of field lines then changes to a more energeti- 
cally favourable configuration, and in doing so will eject plasma along resulting 
outflows as the reconnected field lines relax to a new equilibrium. The outflows 
create a low pressure in the diffusion region which in turn allows continued inflow 
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Sweet-Parker model 





Petschek model 




Figure 1.10: Geometry of the Sweet-Parker (top) and Petschek (bottom) recon- 



nection models, reproduced from Aschwanden (2005). 



of plasma and magnetic field forming a runaway process of reconnection until the 
system comes to rest in a new topology. 



A description of how magnetic reconnection occurs was put forward by Sweet 



(1958 


) and 


Parker 


(1957) 



tion. They estimated the rate of reconnection from a boundary layer analysis 



(top of Figure 1.10). For a reconnection layer of length A and thickness 5, the 



outflow must balance the inflow: 



v in A = v out S 



:i.38i 
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where vi n is the inflow reconnection velocity, and v out is the outflow velocity, 
which by conservation of energy (B^/2^ = pv^ ut /2) is equal to the Alfven velocity 



va = B /y/n p. From Ohm's law (Equation 1.27) the configuration of the straight 
field lines (V x B = 0) outside the layer is E z + vrB x = 0, and inside the 
layer (where there is a large current) is E z = r]J z . Integrating Ampere's law 



(Equation 1.22) around the layer gives B x = pj z 5. Combining the results of 



these two laws gives: 

E z 7]J Z 7] 



(1.39) 



B x /j,J z S fi5 

This can be written in terms of the outflow velocity (or interchangeably the Alfven 
velocity) by: 

The rate of reconnection is then written in terms of the Lunquist number S, 
which is the dimensionless ratio of an Alfven wave crossing timescale to a resistive 
diffusion timescale: 



1 where S = ^ (1.41) 



Vout VS V 

This is the Sweet-Parker result. The problem with the theory is that it predicts 
reconnection to take place on far too slow a timescale to reconcile with observa- 
tions. Consider solar flares, for example, with va ~ 1,000 km s _1 , and / ~ 10 4 km 
resulting in Sweet-Parker reconnection of tens of days when flare energy release 
is actually observed over minutes to hours. 



An extension to the Sweet-Parker model was put forward by Petschek (1964) 



in which the field lines do not have to reconnect along the entire length of the 
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boundary, but could merge over a shorter length A' < A (bottom of Figure 1.10) 



The remaining length of the boundary is occupied by slow shocks, where the mag- 
netic field tension accelerates the plasma to the Alfven velocity. The reconnection 
velocity in the Petschek model may be written: 



J out 



n 



v A /iA' 



Vout 



A 



Vs V A' 



1.42) 



a factor of y A/A' faster than the Sweet-Parker reconnection velocity. Issues 
with the description of the magnetic field and shock formation place limits on the 
plausibility of Petschek's formalism that A' is a free parameter which may be min- 



imised to obtain the maximum reconnection velocity (Kulsrud, 2001). Numerical 



simulations, such as that of Biskamp (1986), are in favour of the Sweet-Parker 



result, unless anomalous resistivity, for example, is considered. In either case, the 
process described by these theories forms the basis of many CME models in an 
effort to understand observations. 



1.3.2 Theoretical CME Models 



It is well known that CMEs are associated with filament eruptions and solar flares 



( Moon et al. , 2002 Zhang & Wang , 2002 ) but the driver mechanism remains elu- 



sive. Several theoretical models have been developed in order to describe the 
forces responsible for CME initiation and propagation, all of which are based on 
the idea that some form of instability must trigger the eruption. These models 
may be explained in terms of the following mechanical analogues, illustrated in 



Figure |1.11 
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Thermal Blast 





Dvnamo 





Mass Loading 



ft 





Telher Release 





Tether Straining 



Figure 1.11: Illustrations of the different mechanical analogues of CME eruptions, 



reproduced from Klimchuk (2001) 
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The Thermal Blast Model proposes that the increased thermal pressure produced 
from a flare overcomes the magnetic field tension and blows it open to cause a 
CME. Observations, however, have shown that not all CMEs are preceded by a 
flare, nor even necessarily associated with a flare at all. 

The Dynamo Model introduces the idea of magnetic flux injection or stressing 
of the field on a time-scale that is too fast for the system to dissipate the mag- 
netic energy before it builds to a critical point and erupts. 

The Mass Loading Model is concerned with the amount of material included in 
the eruption. Prominences, or regions of relatively higher electron density in the 
corona, overlaying a volume of lower density will erupt due to the Rayleigh- Taylor 
instability. 

The Tether Release Model is based on the restraining of the outward magnetic 
pressure by the magnetic tension of the overlying field. As 'tethers' are removed 
a loss-of-equilibrium occurs due to the magnetic pressure/tension imbalance and 
the system erupts. 

The Tether Straining Model is a variant on the tether release model whereby 
an increase in magnetic pressure due to flux injection or field shearing eventually 
overcomes the tension forces and the 'tethers' break and release the CME. 

The tether straining and release models are generally accepted as the most 
likely scenarios for CME initiation, being able to reproduce numerous observa- 
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tions of CMEs through the development of detailed 2D and 3D flux rope models, 

we can 



as discussed below. Within the context of MHD outlined in Section 1.3.1 



describe the solar plasma as a fluid with the assumption that there is negligible 



viscosity, so the motion of the flux rope is governed by the forces of Equation 1.26 



The Lorentz force is thought to be the dominant driver force in modelling CME 
eruptions, certainly during the early stages of propagation, before the drag force 
takes over and the CME propagates with the ambient solar wind through inter- 
planetary space. 



1.3.2.1 Catastrophe Model 



The 2D flux rope model is driven by a catastrophic loss of mechanical equilibrium 



as a result of footpoint motions in the photosphere (Forbes & Isenberg, 1991 



Forbes & Priest, 1995; Isenberg et al. 1993 Priest & Forbes, 1990a, 2000, 2002) 



The model is illustrated in Figure |1.12| as a coronal current filament channel 
and overlying magnetic field lines, in equilibrium due to the balance between the 
magnetic pressure and tension forces acting on the system. The description of the 
model's evolution in time may be split into a storage phase and an eruption phase. 
During the storage phase the footpoints of the system are slowly moved together 
such that the magnetic energy of the flux rope increases. The magnetic tension 
thus increases and causes the flux rope to move downwards, which builds up 
magnetic pressure until a critical footpoint distance is reached where equilibrium 
is lost. In the eruption phase the flux rope is accelerated upwards, stretching the 
magnetic field lines such that a current sheet forms behind it. If reconnection 
occurs in the current sheet then all the energy is released and the upward motion 
of the flux rope is unbounded. Otherwise it will come to equilibrium again at a 
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(d) X = 0.97 




Figure 1.12: Theoretical evolution of a 2D flux rope, reproduced from Forbes 



& Priest (1995). The flux rope footpoint separation A decreases, increasing the 
magnetic pressure until the flux rope becomes unstable and erupts away from the 
surface (b-d). The resulting height evolution of the flux rope is illustrated in (a). 
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greater height, or oscillate about this equilibrium height if it still has energy in 
excess of that required for the initial eruption. The configuration of this model 



outlined in Priest & Forbes (2000) predicts kinematics of the flux rope (provided 
it is 'thin' so that its radius is less than the scale-length Ao) prior to the formation 
of the current sheet (i.e. h/\ < 2) according to: 



h 



VAO 



In 



7T 

2 



2 tan 



■,1/2 



+ ho 



1.43) 



where h is the velocity of the flux rope, ho is the initial perturbation velocity, 
A is the source separation at the critical point, and vao is the Alfven speed at 
h = Xq. The kinematics may be further separated into an 'early' phase when the 
time scale is less than the Alfven time scale (t <C Xq/vaq), which gives: 



/ • \ 3/2 

V37T V A J 



(1.44) 



and a 'late' phase when h/\ Q ^> 1, but |ln/t| is still much smaller than |lna| 
which gives: 



h 



vao 



In 



1/2 



1.45) 



After the formation of the current sheet the system becomes too complicated for 
the kinematics of the CME's continued propagation to be analytically derived. 



1.3.2.2 Toroidal Instability 

An extension of the flux rope model to three-dimensions is illustrated in Fig- 
ure 



1.13 (Chen 1996; Chen & Krall 2003 Kliem & Torok 2006). The eruption 



of the flux rope is triggered by an increase in the poloidal magnetic flux of the 
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(x) B c 




Figure 1.13: A schematic of the 3D flux rope model, reproduced from Chen 



& Krall (2003), an extension of the 2D flux rope in Figure 1.12 when viewed 



end-on as indicated by the arrow from the right. The flux rope is rooted below 
the photosphere, and surrounded by the ambient coronal magnetic field B c and 
plasma density p c . Components of the current density J and magnetic field B are 
shown, where subscripts 't' and 'p' refer to the toroidal and poloidal directions 
respectively. The flux rope has a radius of curvature R, radius of cross-section a, 
apex height Z, footpoint separation sj, and the radial force outward is Fr. 



structure. The 3D flux rope consists of a current channel J and magnetic field B, 
and has major radius R and minor radius a such that for r < a the magnetic field 
lines are helical and can be described by their toroidal and poloidal components, 
but for r > a the field is purely poloidal (J t = 0). The major radius R is fixed, 
and the minor radius increases from a/ at the footpoints to a a at the apex. The 
footpoints are assumed to be immobile because of the high density photosphere 
(~ 10 23 m~ 3 ) relative to the corona (< 10 16 m~ 3 ). The poloidal field B p is also 
highly non- uniform in the photosphere since ^> 1. 
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The model may be directly compared to coronagraph observations as in |Krall 



et al. (2001), where the leading edge of the CME front is located at Z + 2a with a 
width of 4a when viewed end-on, or a width of 2R + 4a when viewed side-on. This 
definition arises from the fact that the poloidal field B p at r = 2a has decreased 
to about half the value of B pa at r = a and is then comparable to the ambient 
coronal field B c . This model sits well with observations, where the CME front 
corresponds to a plasma pileup ahead of the flux rope which appears as a darker 
cavity, and any erupting prominence material is suspended at the base of the flux 
rope and appears as the bright core of the CME. Background parameters such as 
coronal density and solar wind speeds are also specified in the model. The erup- 
tion is initiated by a poloidal flux injection that increases the toroidal current 
for a short period of time, increasing B pa while R does not change significantly, 
such that the radial force becomes more positive and exerts an upward net 
force on the structure. The eruption then proceeds through the corona as the 
external poloidal field decreases sufficiently rapidly in the direction of motion. 



The equation of motion (cf. Equation 1.26) may be written in terms of a radial 



force Fr, a gravitational force F g and drag force Fd, acting to cause the apex 
motion: 

= F R + F g + F d (1.46) 

where the radial force Fr results from the Lorentz magnetic force and pressure 
gradient of the system, and may be written: 
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where It is the toroidal current, c is the speed of light, R is the major radius 
as described above, and fn are the further collective pressure, magnetic and 



geometrical terms to be considered, detailed in Chen & Krall (2003). This for 



malism shows how the toroidal current increase will add to the upward force 
on the structure. The change in current affects the inductance of the flux rope 
(Fr oc I 2 oc L~ 2 ) and so, neglecting the gravity and drag terms, the acceleration 
may be expressed in terms of the geometrical size of the flux rope: 



d 2 Z 
dt 2 



$2 

V 



[Pin (8R/a f )] 



2/R 



1.48) 



Kliem & Torok (2006 ) show how the height of the flux rope during the very initial 



stages of the eruption may be approximated as a hyperbolic function: 



h{r) 



Pi 



sinh(Pir) 



h = H/H - 1< 1 



(1.49) 



where H is the height, and H the initial height, of the flux rope; r is the time 
normalised by the Alfven time; Pq comprises initial parameters on the flux rope 
dynamics; and Pi associates the external magnetic field profile. Their simulations 
show a fast rise and gradual decay phase of the CME accelerations due to the 



toroidal instability. However, Schrijver et al. (2008) demonstrate that tuning the 



initial parameters changes the acceleration profile from a fast initial rise to a more 
gradual rise phase. 
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t= 31902 s 



Figure 1.14: Schematic, and corresponding simulation snapshots, of the main 



stages of the axisymmetric 2.5D breakout model, reproduced from Lynch et al. 



(2008). (a) shows the initial multipolar topology of the system, (b) shows the 



shearing phase which distorts the X-line and causes breakout reconnection to begin, 
(c) shows the onset of flare reconnection behind the eruption that disconnects the 
flux rope, and (d) shows the system restoring itself following the eruption. 
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1.3.2.3 Breakout Model 



In the magnetic breakout model the CME eruption is triggered by reconnection 
between the overlying field and a neighbouring flux system through the shearing 



1.14 


(Antiochos et al. 


1999 


Lynch 



et al., 2004; MacNeice et al. 2004). It starts by shearing a potential field con- 
figuration consisting of a central arcade which will become the CME, two side 
arcades, and an overlying arcade, with a magnetic X-line separating the different 
topologies. This shearing adds magnetic pressure to the inner flux system and 
causes it to expand and distort the overlying field at the X-line, eventually form- 
ing a current sheet. As the current sheet grows, reconnection begins, transferring 
flux to the neighbouring arcades and creating a passage for the CME release 
as the central arcade erupts. A current sheet also forms beneath the erupting 
sheared field, creating a disconnected flux rope that escapes, and this is associ- 
ated with flare reconnection. An increase in the rate of outward expansion drives 
a faster rate of breakout reconnection, yielding the positive feedback required for 
an explosive eruption. At given distances from the Sun, the simulation, which is 



intrinsically 2.5D but has recently been extended to 3D by Lynch et al. (2008), 



produces a number of key observational properties to test against data (van der 



Hoist et al, 2007). Simulations run by Lynch et al. (2004) produced kinematics 



of the CME front which showed constant acceleration, and they fit quadratics to 
the height-time data of the form: 



h(t) = h + v (t - t ) + ^a(t - t ) 2 



(1.50) 
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However, the height-time profiles from the 3D simulations by Lynch et al. (2008) 
showed a more complex kinematic profile with separate rising and breakout phases 
of acceleration, proving a better match with observations. The velocity increases 
linearly from zero during the initial shearing and breakout phase, followed by 
an acceleration peak during the second stage of reconnection behind the CME. 



DeVore & Antiochos (2008) find a three-phase acceleration profile with a similar 
initial slow and fast acceleration profile during the breakout and flare recon- 
nections, followed by a short interval of fast deceleration as the magnetic field 
configuration is reformed after the eruption. 

1.3.3 CME Observations 

In order to test the validity of the theoretical CME models, and subsequently 
understand the forces governing their eruption through the solar atmosphere, 
comparisons must be made with observations. Generally the kinematics of events 
are determined from white-light coronagraph images, obtained through Thomson 
scattered emission of the plasma, by tracking the structure as it moves through 
the field-of-view. 

1.3.3.1 Thomson Scattering 

The low density, optically thin, plasma of the corona and solar wind is observable 
in white light through the process of Thomson scattering, whereby photons from 



the Sun are scattered by the free electrons of the coronal plasma (Billings, 1966 



Minnaert 1930 van de Hulst, 1950). Essentially, when light is incident on an 
electron, the electric field of the light waves will cause the electron to accelerate 
and re-radiate light in the plane perpendicular to the incident wave (illustrated in 



44 



1.3 Coronal Mass Ejections 





Figure 1.15: The Thomson scattering geometry for a single electron, reproduced 
from Howard Sz Tappin (2009). a) shows the electron and incident light with 
different observer positions indicated, b), c), d) show the resultant unpolarised, 
partially polarised, and polarised scattering of light seen by observers Oi, O2, O3 
respectively. 



Figure 1.15 ). Depending on the angle x to the observer, the scattered light may be 



unpolarised (Figure |i.l5[ j), partially polarised (Figure [1.15H ), or polarised (Fig- 
ure 



1.15f i). The tangential component of the scattered light intensity is isotropic, 



while the radial component varies as cos 2 x, resulting in the following expression 
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for the differential cross-section (Jackson, 1975): 



da If e 2 x ' 



(7-^—2) ( 1 + cos ^) (1.51) 



du 2 \AneQm e c 2 

where du is an element of solid angle at scattering angle \i an d eo is the permit- 
tivity of free space. Since the Sun is neither a point-source, nor is the intensity of 
light uniform across it, it is necessary to consider the three-dimensional geometry 
of the Thomson scattering process in detail to appreciate the effects it will have 



on CME observations, outlined in detail in Billings (1966) and Howard & Tappin 



(2009). In a full description of the corona the density profile must be considered, 
as well as the effect of limb darkening on the Sun whereby the light intensity from 
the photosphere decreases across the disk according to: 

I = 7 (1 -m + mcosV') (1-52) 



where I is the intensity observed at angle ip from the radial vector, 7 is the radial 
intensity, and u is the limb-darkening coefficient. 

Describing the Thomson sphere as the locus of all points that make an angle 
of x — 90° between the line-of-sight and the vector from the Sun to the scattering 



point P (Figure 1.16), the total intensity of scattered light is governed by three 
terms: 

1. The scattering efficiency which is minimised on the Thomson sphere. 

2. The incident intensity which is maximised on the Thomson sphere since 
that is where the line-of-sight is closest to the Sun. 

3. The electron density in the scattering region which is maximised on the 
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Thomson Scattering Geometry 



Thomson Surface 




Observer 



R = Sun-Observer Distance 

r = Heliocentric Distance of Scattering Center 

d = Impact Radius 

p. = Elongation 

I = Longitude (rel. to Solar Limb) 



Figure 1.16: Schematic of the Thomson surface, being the sphere of all points 
which are located at an angle of 90° between the Sun and the observer, reproduced 



from Vourlidas & Howard ( 2006 ) . An example line-of-sight is shown for an electron 
at point P, with radial distance r from the Sun, at longitude 4> relative to the solar 
limb. 
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Thomson sphere since the solar wind density drops off with radial distance 
from the Sun. 

The combination of these effects makes the Thomson sphere an important con- 
sideration when interpreting CME observations, especially out through the wide- 



angle fields-of-view of the heliospheric imagers on the STEREO spacecraft. VourL 



idas & Howard (2006) note the following: 



• CMEs that propagate along the solar limb and appear bright in near-Sun 
coronagraphs are unlikely to be detectable further out in the heliosphere. 

• Frontside events are always brighter than their backside counterparts, and 
ones at intermediate angles will exhibit approximately constant levels of 
brightness over a wide range of heliocentric distances. 

• The sky-plane assumption holds well for brightness observations out to at 
least ~ 70 R . 

If an expanding CME front moves along or crosses the Thomson sphere during 
its propagation, its observed brightness can change with regard to the changing 
location of its intersection with the sphere, as a result of the discussion above. 
This has implications for how the observed CME kinematics and morphology may 
be affected, certainly at large elongations from the Sun, and how CMEs may look 
quite distinct from different observers' points of view. 

An example of how the CME brightness can change due to the projection 
effects of an image and the consideration of Thomson scattering geometry is 



shown in Figure 1.17 It is based upon a hollow, unkinked, flux rope model with 



an electron density profile that is peaked toward the outer surface of the flux rope, 
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Figure 1.17: (a) Synthetic coronagraph image of the model flux rope, deduced 
from the integrated line-of-sight electron density, including a consideration of the 
Thomson scattering geometry, (b) The flux rope model represented by a mesh of 
helical field lines. Point A marks the true apex of the model. The axes are in units 
of R . Reproduced from Chen et al. ( 2000 ) . 



described in Chen et al. (2000). To obtain the synthetic coronagraph image, the 
electron density is integrated in the direction perpendicular to the plane-of-sky 
weighted by the angular dependence of Thomson scattering. The line-of-sight 
integrated electron density is greatest where the helical lines bend around the 
toroidal surface or where they bunch up. Thus the most prominent feature of the 
synthetic image is the bright outer rim, and the projection of the interior of the 
flux rope conveys the relatively dark cavity (no prominence material is included 
hence there is no bright core). In reality, non-linear densities and kinked magnetic 
field topologies add to the complexities of these observations and increase the 
difficulties of tracking CMEs as they propagate and evolve through sequences of 
images. 
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Figure 1.18: Three example CME observations (off the limb, partial halo, and 
full halo) and schematics showing how measurements of the plane-of-sky velocity 
Vps and expansion velocity V exp are skewed by projection effects, reproduced from 



Schwenn et al. 



(2005). 



1.3.3.2 CME Kinematics & Morphology 



A large number of CMEs have been studied since the advent of space-borne 
coronagraph observations, and they show speeds varying from tens up to a few 
thousand kilometres per second. Many of these exhibit a general multiphased 
kinematic evolution: CMEs tend to have an initial rise phase with possible high 
acceleration, and a subsequent constant-velocity cruise phase with another pos- 
sible low acceleration or deceleration in their continued propagation through the 



heliosphere. This initially led Sheeley et al. (1999) to distinguish CMEs as either 



'gradual' if their initial acceleration is low, or 'impulsive' if it is high. However, 
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statistics on a large sample of events do not show such a clear distinction but 
do indicate that slow CMEs tend to result from prominence lift-offs or streamer 
blowouts and speed up to the solar wind speed, while fast CMEs tend to result 



from flares and active regions and slow down to the solar wind speed (Gonzalez- 



2003 


Gopalswamy et al. 


2000 


Moon et al. 



analyses can provide a general indication of CME properties, for example Zhang 



& Dere (2006) study 50 CMEs and find an average acceleration of 330.9 m s 



with an average duration time of 180 minutes, and Gopalswamy et al. (2000) 
study 28 CMEs and derive a formula for their acceleration a related to their 
initial speed u by a = 1.41 — 0.0035m. However, plane-of-sky projection effects 
mean the measured kinematics are not representative of the true CME motion 



(Figure 1.18), with Wen et al. (2007) deducing that the error in CME leading- 
edge measurements grows roughly with the square of the distance from Sun centre 
within the first few solar radii and then varies approximately with the square root 



of the distance past ~5 R . In an effort to overcome plane-of-sky effects, dal 



Lago et al. (2003) use a sample of 57 limb CMEs to derive an empirical relation- 



ship between their radial and expansion speeds as V ra d = 0.88V^™, and Schwenn 



et al. (2005) similarly use 75 events to derive a formula for their transit time to 



Earth T tr = 203 — 20.77 In (V exp ). However, Vrsnak et al. (2007) show the inher- 
ent difficulties in performing and trusting such corrections for CME projection 
effects. This means individual CMEs must be studied with rigour in order to sat- 
isfactorily derive the kinematics and morphology to be compared with theoretical 
models. Examples of such rigourous analyses are seen in the following kinematic 
studies. 



Zhang et al. (2001 ) study the temporal relationship between CMEs and flares, 
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Figure 1.19: The CME velocity profile, and associated soft X-ray flare profile, for 



the event on 11 June 1998, reproduced from Zhang et al. (2001). The profiles in 



dicate a three-phase scenario of CME evolution: initiation, impulsive acceleration, 
and propagation. The datapoints are from LASCO/C1 (asterisks), C2 (triangles) 
and C3 (squares). 

with an emphasis on the three-phased scenario of CME evolution: initiation, im- 
pulsive acceleration, and propagation (Figure 1.19). The height range of ~1- 



3 R was not possible to observe previous to the launch of the SOHO/LASCO 
suite. They obtained height-time measurements from running difference images, 
and compare the derived CME velocity profile with the soft X-ray flux of its as- 
sociated flare, revealing a strong connection between the two phenomena. 



Following the loss of the innermost LASCO/C1 coronagraph, Gallagher et al 
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Figure 1.20: Height-time evolution of a CME observed with TRACE, UVCS 



and LASCO/C2 on 21 April 2002, reproduced from |Gallagher et aL| ( |2003| ). The 
datapoints have ± 5 pixel errorbars, and the best fits for an exponential varying 
and constant acceleration are plotted for comparison. 



(2003) include the use of low corona EUV observations for a CME on the 21 April 



2002 in order to track a CME from the very low corona out to almost 30 R . 
They consider the separate sets of observations collectively in their fitting, which 
is not strictly appropriate since the observations are due to different mechanisms 
of emission from hot plasma in EUV images and Thomson scattering of light in 



coronagraph images. Following Alexander et al. (2002) they first consider the 



simple case of constant acceleration in the early CME evolution (Figure 1.20), of 
the form: 

h(t) = h + v t+ -at 2 (1.53) 

and find that the best-fit to the data points does not adequately represent the 
observed profile. They then consider an exponentially varying acceleration, of 
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Figure 1.21: (a) Height-time, (b) velocity and (c) acceleration profiles for the 
CME on 21 April 2002, and (d) the GOES-10 soft X-ray flux for the associated 



XI. 5 flare during the interval 00:47-03:20 UT, reproduced from Gallagher et al. 



(2003). A three-point difference scheme is used on the data points, and a first- 



difference scheme is plotted with filled circles. The solid line is the best fit of the 



exponentially increasing and decreasing acceleration of Equation 1.55 
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the form: 

h{t) = h + v t + a T 2 exp(t/T) (1.54) 

which does provide an acceptable fit to the datapoints. They track the CME 
front through running-difference images and measure its height-time profile from 



~20 Mm to over 10 4 Mm above the Sun's surface, almost 30 R (Figure 1.21), 
resulting in a rising and falling acceleration. They model this resulting accelera- 
tion profile with a combined function of exponentially increasing and decreasing 
terms: 



a(t) 



1 1 



1.55) 



_a r exp(t/r r ) a d exp {-t/r d )_ 

where a r and are the initial accelerations, and r r and the e-folding times 
for the rise and decay phases, and show that the function provides a close ap- 
proximation to the trend in the datapoints (for completeness they show both 
a three-point difference and first-difference scheme). They determine an early 
acceleration peak of ~ 1,500 km s _1 and show it to correspond to the duration 
of the soft X-ray rise phase of the associated XI. 5 flare, having implications for 
either the thermal blast model, or a magnetically-dominated process such as re- 
connection, causing the eruption. Either way they highlight the importance of low 
coronal observations, and precise measurements, for revealing different regimes 
of CME propagation. 



Maricic et al. (2004) present a study of the initiation and development of 



a limb CME on 15 May 2001, tracking its leading edge, cavity and associated 



prominence from 0.32 R out to ~30 R (Figure 1.22). They distinguish a pre- 
acceleration characterised by the slow rising motion of the prominence, suggestive 
of an evolution of the system through a series of quasi-equilibrium states. They 
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Figure 1.22: The kinematics of the CME leading edge (circles), cavity (pluses) 



and prominence (triangles) of the 15 May 2001 event, reproduced from Maricic 



et al. (2004). (a) shows the height-time profiles, (b) shows the distance between 
the leading edge and the prominence (thick line), the cavity and the prominence 
(thin line), and the leading edge and the cavity (dashed line), (c) shows the veloc- 
ities determined by forward- difference technique upon the smoothed height-times, 
(d) shows the onset of acceleration against height, with a straight line fit. The 
horizontal bar between (a) and (c) indicates the period of the soft X-ray burst 
(dotted - precursor, full ■ 



mam rise 
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offer numerous possible explanations for this: a slow shearing or merging motion 



of the arcade footpoints (e.g., Lin et al. , 2004 Priest & Forbes, 1990b), twisting of 



the embedded flux rope (Vrsnak, 1990), emerging azimuthal flux (Chen & Krall 



2003), mass loss from the prominence body (Vrsnak, 1990), etc. This is followed 



by a rapid acceleration onset of the prominence, 380 ±50 m s~ 2 , simultaneous 
with the CME leading edge acceleration, 600 ±150 m s~ 2 . This simultaneity 
rules out the scenarios in which the prominence motion is merely a consequence 
of the disruption of the overlying magnetic structure, or that the prominence 
eruption itself drives the upper parts of the system (at least for this event, and 
the authors suggest a careful inspection of a larger sample of events). There is 
also a simultaneous soft X-ray burst which they deem an 'acceleration precursor' 



as postulated in the catastrophic evolution of a flux rope/arcade system (Lin & 



Forbes, 2000). A possible gradual deceleration of the CME in its later stages is 



also noted (-23 ± 1 m s 2 for the leading edge) and attributed to the possible 



effects of aerodynamic drag in the solar wind discussed in Cargill et al. (1996), 



Vrsnak (2001) and Vrsnak et al. (2004). Maricic et al. (2004) also measure the 



relative height difference between the CME leading edge, cavity and prominence, 
and report a kind of self-similar expansion of the event. 



Temmer et al. (2008) study the CME kinematics and associated flare flux pro- 



files for two fast events of 17 January 2005 and 6 July 2006 (the latter is shown 



in Figure 1.23). Using running-difference techniques they track the rising loops 
of the flares and consider them as the early height-time profile of the subsequent 
CME eruption. This results in a clear early acceleration peak for each event, of 
~ 4,400 m s -2 and ~ 1,100 m s~ 2 respectively. They find that the peaks of the 
hard X-rays correspond with the peak accelerations of the CMEs, indicating a 
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Figure 1.23: Kinematics of a CME observed with TRACE, EIT and LASCO/C2 
and C3 on 6 July 2006, plotted over the full time range on the left, and zoomed in 



for the duration of the acceleration on the right, reproduced from Temmer et al. 



(2008). The bottom panels include the hard X-ray flux of the associated flare 



(shown in red), peaking simultaneously with the acceleration. 

strong feedback relationship indicative of magnetic reconnection occurring in a 
possible current sheet behind the CME to drive the eruption under the dominant 
Lorentz force. 



More recently Lin et al. (2010) performed the same type of analysis on two 



CMEs in order to specifically test the kinematics of each against the predictions of 
theoretical models (catastrophe, breakout and toroidal instability). Tracking the 
CME front through running-difference images, the events of 17 December 2006 
and 31 December 2007 were found to have peak accelerations of ~60 m s~ 2 and 
~ 1,500 m s -2 within the first ~3 R Q height range, simultaneous with the soft 
X-ray flux. The authors show the difficulty in distinguishing any single model as a 
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basis for the eruption, since the predicted fits all fare relatively well in producing 
a profile within the scatter and error of the measured kinematics. The authors 
also acknowledge that the data is a two-dimensional projection of actual motion, 
and three-dimensional kinematics of CMEs are necessary for improved accuracy. 

The conclusions of such studies rely heavily on the kinematic uncertainties, 
since they are vital for providing a distinction between different models that at- 
tempt to fit the profiles. The above discussion highlights the importance of the 
image resolution, cadence and projection effects when performing CME analy- 
ses. It is also apparent that user-specific biases and different numerical methods 
for determining the kinematics and morphology can affect the resulting profiles 
and hence their interpretation. With these considerations in mind, the current 
aims of the community have been to develop methods for overcoming both the 
biases of individual users, and single viewpoint observations, in order to robustly 
determine the kinematics and morphology of CMEs with the greatest possible 
precision. These aims motivated the launch of the STEREO mission to provide 
twin-viewpoint observations of CMEs, with low-cadence observations of the inner 
corona, and off-limb imaging suites with fields-of-view extending out to the orbit 
of Earth at 1 AU. 

1.4 Outline of Thesis 

The work presented in this thesis improves the understanding of the kinematics 
and morphology of CMEs as they propagate through the solar corona. To date 
the quantification of these properties has been subject to various sources of error, 
most notably as a result of the innate difficulties in tracking CMEs with tra- 
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ditional image processing techniques. Moreover the projected two-dimensional 
nature of coronagraph observations has been a persistent source of error that 
scientists have striven to overcome when studying CMEs. This thesis outlines 
new methods of multiscale image processing and ellipse characterisation of the 
CME front in coronagraph observations so as to reduce the error on the derived 
kinematic and morphological parameters. This is extended to stereoscopic image 
data whereby an elliptical tie-pointing methodology for reconstructing the CME 
front in three-dimensions leads to a study of its true kinematics and morphology 
as it propagates into the heliosphere. 

Chapter 2 details the space-borne instrumentation used in this work for ob- 
serving CMEs through the solar corona and heliosphere. Chapter 3 discusses the 
current CME catalogues in use and their limitations, leading to the implemen- 
tation of new multiscale techniques and an ellipse characterisation for studying 
CME propagation. The efforts to automate this process for the development of a 
new cataloguing database are also discussed. These methods are the basis of the 
studies undertaken in Chapter 4, where a selection of CME events are presented 
and their kinematics and morphologies discussed in light of theory. Chapter 5 
presents an important extension of these methods to perform a three-dimensional 
reconstruction of a CME front, which overcomes the issues of projection effects 
and so provides insight into the true kinematics and morphology of the eruption. 
A detailed discussion of the results and conclusions drawn from this work are 
also presented. Finally Chapter 6 presents the main conclusions of the thesis and 
details possible future work that could follow on from these new developments. 
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Instrumentation 



CMEs were initially observed in the early 1970s with the launch of the first space- 
borne coronagraph onboard the seventh Orbiting Solar Observatory (OSO 7) pro- 
viding daily white-light coronal images with a field-of-view 2.8-10 R , observing 



about a dozen CMEs from 1971-1974 (Koomen et al. 1975 Tousey & Koomen 



1972). Skylab, a U.S. space station launched in 1973, housed a coronagraph 



with field-of-view 1.5-6 R that imaged the corona every 6-8 hours and con- 
clusively established these transient ejections as a common occurrence, observing 



~ 100 in 1973/74 (MacQueen et al. 1974). Following Skylab, several more coro- 
nagraphs were flown: SOLWIND on satellite P78-1 observed over 1,500 CMEs 
from 1979-1985 ( |Sheeley et al. 1980); the High Altitude Observatory Coron- 
agraph/Polarimeter onboard Solar Maximum Mission (SMM) observed ~ 1,350 



CMEs from 1980-1989 (MacQueen et al, 1980); the Large Angle Spectromet- 



ry Coronagraph suite (LASCO; Brueckner et al., 1995) onboard the Solar and 
Heliospheric Observatory (SOHO; Domingo et al. |1995 ) has observed thousands 
of CMEs from 1995 to present; and the COR1/2 coronagraphs of the Sun-Earth 
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Connection Coronal and Heliospheric Imaging suite (SECCHI; |Howard et al. 



2008a) onboard the Solar Terrestrial Relations Observatory (STEREO; Kaiser 



et al, 2008) have been observing CMEs from 2006 to present. Data from the 



instruments of SOHO/LASCO and STEREO/SECCHI are used in the research 
presented throughout this thesis. 



2.1 SOHO/LASCO 

SOHO is a joint European Space Agency (ESA) and National Aeronautics and 
Space Administration (NASA) mission, launched on 2 December 1995 to under- 
take scientific investigations of (1) helioseismology, the study of the interior solar 
structure, and (2) the physical processes that account for the heating and accel- 
eration of the solar wind (or, more broadly, the nature of evolutionary change in 
the Sun's outer atmosphere). SOHO is situated in orbit about the Lagrangian LI 
point approximately 1.5 x 10 6 km sunward of the Earth for an uninterrupted view 
of the Sun. Onboard are twelve complementary science instruments: three helio- 
seismology experiments to probe the Sun's inner structure through measurements 
of solar oscillations; three solar wind experiments to measure the in-situ proper- 
ties of the ambient wind (densities, speeds, charge states, etc.); and six telescopes 
and spectrometers to study the solar disk and atmosphere. A schematic of SOHO 



and its instrument suites is shown in Figure 2.1 

The LASCO instrument suite is a set of three coronagraphs CI, C2 and C3 
that image the solar corona from 1.1-3, 1.5-6 and 3.7-30 R respectively (how- 
ever the CI coronagraph has not been in operation since 1998 when contact with 
the SOHO spacecraft was lost for several weeks). The coronagraph was invented 
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Figure 2.1: A schematic of the SOHO spacecraft and onboard instrument suites, 



reproduced from Domingo et al. (1995). 



by French astronomer Bernard Lyot in 1939 to artificially eclipse the Sun for 
observing the solar corona. It essentially blocks light rays from the centre of 
the telescope field-of-view by occulting the solar disk, in order to increase the 
relative intensity of the surrounding coronal light which is on the order of 10 6 
times fainter. The externally occulted Lyot coronagraph design of LASCO/C2 



and C3 is illustrated in Figure 2.2 The top diagram demonstrates how the op- 
tical assembly images the coronal light, while the bottom diagram demonstrates 
how stray light is suppressed. Light is incident through aperture AO where the 
external occulter Dl eclipses the solar disk. The light then enters aperture Al 
and is focused by the objective lens 01, through the field stop, onto the inner 
occulter D2 which anodises the bright fringe of the external occulter. Field lens 
02 then collimates the light onto the Lyot stop A3 that intercepts the light rays 
diffracted off the entrance aperture Al. A relay lens 03 is placed behind A3 to 
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Figure 2.2: The externally occulted Lyot coronagraph, reproduced from Brueck- 
ner et al. ( 19951), showing: front aperture AO and external occulter Dl; entrance 



aperture Al and objective lens 01; the field stop; inner occulter D2 and field lens 
02; Lyot stop A3 and relay lens 03 with Lyot spot; filter and polariser wheels 
F/P; and the focal plane F. 



focus the coronal image on to the plane F. 03 contains the Lyot spot for inter- 
cepting residual diffracted light from Dl and ghost images created by 01. In front 
of the focal plane F are the colour filters and linear polarising filters F/P. The 
colour filters distinguish specific bandpasses of the coronal light, in the ranges 
400-850 nm for C2 and 400-1050 nm for C3. The polariser wheel is used to 
obtain total brightness B or polarised brightness pB images through combina- 
tions of polariser positions I a = —60°, h = 0°, and I c = 60°, according to the 
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Figure 2.3: A LASCO/C3 image of the solar corona out to ~ 30 R . 
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equations (Billings, 1966): 



B = -{Ia + h + Q 

pB = ^[(I a + I b + I c ) 2 -3(I a I b + I a I c + I b I c )} 1/2 



(2-1) 
(2.2) 



A CCD is placed at the focal plane F and the final images are 1024 x 1024 pixels, 
subtending an angle of 11.4 arcseconds per pixel in C2, and 56 arcseconds per 



pixel in C3 (see Section 2.3 for CCD details). A sample LASCO/C3 image is 



shown in Figure 2.3 



2.2 STEREO/SECCHI 



STEREO is the third mission in NASA's Solar Terrestrial Probes program (Kaiser 



et al, 2008). It was launched on 25 October 2006, and employs two nearly 
identical space-based observatories; one ahead of Earth in its orbit, and the other 
behind, separating at ± 22° each year. This arrangement provides the first ever 
stereoscopic observations of the Sun and inner heliosphere. The main objectives 
of STEREO are to: 

• Understand the causes and mechanisms of CME initiation. 

• Characterise the propagation of CMEs through the heliosphere. 

• Discover the mechanisms and sites of energetic particle acceleration in the 
low corona and the interplanetary medium. 

• Improve the determination of the structure of the ambient solar wind. 
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Figure 2.4: Payload diagram of one of the STEREO spacecraft, indicating the 
positions of the four instrument suites onboard: Sun-Earth Connection Coronal 
and Heliospheric Imagers (SECCHI); In-situ Measurements of Particles and CME 
Transients (IMPACT); Plasma and SupraThermal Ion Composition (PLASTIC); 
STEREO/WAVES radio burst tracker (SWAVES). 
Image credit: stereo.gsfc.nasa.gov. 



STEREO hosts four instrument suites to achieve this, as illustrated in Fig- 



ure 



2.4 The SECCHI suite comprises five scientific telescopes: firstly the Sun 



Centred Imaging Package (SCIP) consisting of an Extreme Ultraviolet Imager 
(EUVI) of the solar disk out to 1.7 R and two coronagraphs (COR1/2) with 
fields-of-view 1.4-4 and 2-15 _R ; and secondly the Heliospheric Imagers (HI) 
consisting of two wide-angle visible light imagers positioned on the sides of the 
STEREO spacecraft for fields-of-view extending out to Earth at 1 A.U. (astro- 
nomical unit, based on the distance from the Earth to the Sun which is approxi- 
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Figure 2.5: Diagram of the EUVI in the STEREO/SECCHI suite, reproduced 



from Howard et al. (2000) 



mately 1.49 xlO 8 km). 



2.2.1 EUVI 

The Extreme Ultraviolet Imager is a normal-incidence Ritchey-Chretien telescope 
that images the solar disk out to 1.7 R at four wavelengths of emission that 



span a temperature range of 0.1 to 20 MK (Wuelser et al, 2004). Radiation 
from the Sun enters through a thin aluminium filter of 150 nm thickness that 
suppresses most of the ultraviolet, visible and infrared wavelengths of light. The 
radiation passes one of four quadrants that are each optimised for one of the EUV 



wavelength lines (listed in Table 2.1). The primary and secondary mirrors direct 
the light through a filter wheel that has a redundant thin-film aluminium filter 
to remove the remainder of the visible and IR radiation. A shutter in the path 
controls the exposure time, and 2048 x 2048 pixel images are produced by the 
CCD subtending an angle of 1.6 arcseconds per pixel. 
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Table 2.1: Summary of EUVI wavelengths. 



Principal emission lines Wavelength 

Fe IX 172 A 

Fe XII 194 A 

Fe XV 284 A 

He II 304 A 



2.2.2 COR1 



CORl is a classic Lyot internally occulting refractive coronagraph (Thompson 



et al, 2003). Light enters through the front aperture of the telescope and is 



focused by the objective lens onto the occulter, with a series of baffles in place 



to minimise scattering of light within the telescope (Figure 2.6). The occulter is 
cone shaped to reject light from the centre of the field-of-view into a surrounding 
light trap. The field lens focuses the rest of the light down the telescope to the 
Lyot stop which removes light diffracted by the edge of the front aperture. A 
Lyot spot is also glued to the doublet lens immediately behind the Lyot stop in 
order to remove any ghosting of the objective lens. Two doublet lenses focus the 
light onto the CCD detector, with a bandpass filter 10 nm wide (centred on the 
Ha line at 656 nm) and a linear polariser in between them. To extract both total 
brightness B and polarised brightness pB images, three sequential images are 
taken with polarisations of I a = —60°, If, = 0°, and J c = 60°, and combined using 



Equations 2.1 and 2.2 The cut-on frequency (at 350 nm) of light through CORl 
is set by the transmission of the BK7-G18 glass in the objective lens, and the 
cut-off frequency (at 1100 nm) is set by the band gap of the silicon CCD detector. 
The final images are 2x2 binned onboard to 1024 x 1024 pixels, subtending an 
angle of 7.5 arcseconds per pixel. The field is unvignetted except for a small 



69 



2. INSTRUMENTATION 




70 



2.2 STEREO/SECCHI 



area around the edge of the occulter and near the field stop in the corner of the 
images. The average radial brightness profile for both instruments is well below 
1CT 6 B/B though some discrete ring-shaped areas of increased brightness in the 
CORl-Behind instrument are caused by features on the front surface of the field 
lens. 



2.2.3 COR2 

COR2 is an externally occulted Lyot coronagraph, similar to the LASCO/C2 and 



C3 telescopes (Figure 2.7). An array of internal baffles sits behind the external 
occulter to reduce stray light entering the telescope, and an internal occulter and 
Lyot stop minimise diffraction effects. The final images are produced at three 
different polarisations as in COR1 for creating total brightness and polarised 
brightness images. The final images are 2048 x 2048 pixels, subtending an angle 
of 14.7 arcseconds per pixel. The image is vignetted throughout the field-of-view, 
at a level of 40-50% around the occulter pylon, and reaching a minimum of 20% 
at about 10 R before increasing again towards the image edge. 



2.2.4 Heliospheric Imagers 



The Heliospheric Imagers (HI1/2; Eyles et al, 2009) are two small, wide-angle, 
visible-light camera systems mounted to the side of each STEREO spacecraft to 



image along the Sun-Earth line from elongations of 4-88.7° (Figure 2.8). This 
has provided several new opportunities for CME research, notably the ability to 
track their evolution as they propagate through the inner heliosphere and po- 
tentially impact at Earth or one of either STEREO spacecraft which allows a 
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2.2 STEREO/SECCHI 
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Figure 2.8: Diagram of the heliospheric imagers HI-1/2 in the STEREO/SECCHI 



suite, reproduced from Howard et al. (2000). 



comparison of in-situ data and white-light imagery of CMEs. The basic design of 
the HI comprises a number of occulting baffles that achieve the required level of 
light rejection for imaging the low intensities required to observe CMEs at large 



elongations (Figure 2.9). In order to image the low intensity CME signal suffi- 
ciently above the stellar background and zodiacal light (F-corona) whilst avoiding 
saturation by cosmic rays, a series of short exposures is taken and individually 
scrubbed of cosmic ray hits before being summed together. This increases the 
CME signal above the noise level of the relatively static F-corona, which may 
then be subtracted to reveal the CME intensity. The combination of summing 
(30 images for HI1, and 99 images for HI2) and 2x2 binning increases the signal- 
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Figure 2.9: The intensity profile of a CME compared to the K & F coronae 
observed at elongations up to 90°, and the corresponding fields-of-view of the He- 



liospheric Imagers (HI1/2), reproduced from Howard et al. (2000). 



to-noise ratio by about 14 times. The final images are 2x2 binned onboard to 
1024 x 1024 pixels, subtending an angle of 70 arcseconds per pixel for HI1 and 
4 arcminutes per pixel for HI2. 



2.3 CCD Detectors 

A charge-coupled device (CCD) is used in the LASCO and SECCHI instruments 
for detecting the incident photons and converting them to a digital output to 
generate images. Essentially a CCD converts light into electrons which are read 
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Figure 2.10: Illustration of a thick front-side illuminated CCD. 

Image credit: www.ing.iac.es. 



and converted into numeric values used to display image intensities. The CCD 
is a small silicon chip divided into a grid of cells, or pixels. The electrons in 
the silicon atoms lie in discrete energy bands. In the ground state the outermost 
electrons lie in the valence band and can be excited to the conduction band by the 
absorption of a photon, via the photoelectric effect, leaving behind a 'hole'. In a 
CCD an electric field is introduced to prevent recombination of the electron-hole 
pair. Thus an electric charge is accumulated proportional to the light intensity at 
that location. The charge is read out pixel-by-pixel to a charge amplifier which 
converts it to a voltage, then this voltage is digitised and stored in memory. 



A thick front-side illuminated CCD (Figure 2.10) is cheap to produce, but 
because photons are incident at the surface electrodes they can be reflected or 
absorbed, which gives low quantum efficiency (a measure of the percentage of 
photons detected: QE = N e /N v ). The LASCO/C2 and C3 detectors are front- 
side illuminated Textronix CCDs that have a quantum efficiency of about 0.3-0.5 
in the 500 to 700 nm spectral range. They are 1024 x 1024 pixels in size, each 
pixel being a square measuring 21 /im on a side. 
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Figure 2.11: Illustration of a thinned back-side illuminated CCD. 
Image credit: www.ing.iac.es. 



To increase the quantum efficiency back-side illumination is used so the elec- 
trodes do not obstruct the photons. But the silicon in a back-side illuminated 
CCD must be chemically etched down (thinned) to a thickness of about 15 /im, 



which is an expensive process (Figure 2.11). Silicon also has a high refractive 



index leading to strong photon reflection. It must therefore be coated with an 
anti-reflective material with a refractive index less than that of silicon (3.6) and 
preferably with an optical thickness of 1/4 at a chosen wavelength of 550 nm 
(close to the middle of the optical spectrum). Hafnium dioxide is regularly used 
to significantly reduce the reflectivity of the CCD. Due to their high quantum 
efficiency, almost all current astronomical CCDs are thinned and back-side illumi- 
nated. Each of the SECCHI instruments uses a back-side illuminated E2V 42-40 
CCD detector, that has a quantum efficiency of roughly 0.8 at 500 nm, 0.88 at 
650 nm, 0.64 at 800 nm, and 0.34 at 900 nm. They are 2048 x 2048 pixels in size, 
each measuring 13.5 /im on a side. This CCD has an operational temperature 
range of 153-323 K (-120-50° C). 

Sources of noise in CCD imaging must be noted when performing image anal- 
ysis. Thermal noise, or dark current, is due to thermal excitations of electrons in 
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the CCD. A dark frame must be generated to correct for thermal noise by taking 
a closed shutter exposure of some known duration to study the effects on the re- 
sultant image, though this form of noise is minimal for space-borne instruments 
operating at temperatures of < 200 K. 

Hot pixels can result from energetic particles or cosmic rays causing ionisation 
in the silicon, because the resulting free electrons from these hits are indistinguish- 
able from ones that are photo- generated. 

CCD read-out noise can occur when charge is converted to voltage since elec- 
tronic amplifiers are not perfect. A high charge transfer efficiency is also impor- 
tant during shift operations in the read-out process to minimise count errors. 

Calibrations of CCD images must be performed to remove imperfections. 
CCDs are not always linear (measuring one count for one photon incident). A 
flat-field calibration removes variations in sensitivity across the surface of the 
CCD, due to silicon or manufacturing defects and vignetting effects. Flat-field 
images are normally generated in the lab by taking an exposure when the CCD 
is evenly illuminated by a light source, and dividing this into future images for 
linearity. 

Similar to dark frames, bias frames may also be generated. A bias frame is 
a zero duration exposure taken with no light incident on the CCD (the shutter 
remains closed). Thus structures which appear in bias frames £1X6 duS du result of 
defects in the CCD electronics and must be removed from future images. 

The charge capacity of a CCD pixel is limited and when full it can overflow, 
leading to blooming. While this is somewhat unavoidable when taking long ex- 
posures, especially if a bright star or comet comes into view for example, most 
CCD design ensures blooming only occurs in one direction. 
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2.4 Coordinate Systems 

The pixel coordinates from the CCD images must be transformed to the relevant 
coordinate system for studying and interpreting observations, especially when 
comparing images from multiple viewpoints (such as STEREO and SOHO as 



discussed in Section 5.3.1). First the pixel coordinates (pi, p 2 , P3, ■■■) must be 
transformed to intermediate world coordinates (x\, x 2 , £3, ...), meaning they are 
converted into the relevant units (e.g., metres or arcseconds) but are not nec- 
essarily corrected for the reference point of the observations nor geometric or 
projection effects: 

»>.,•(/',• - r ; ) (2.3) 



N 

Xi Si 



where Si is the scale function, N is the number of axes, rriij is the transformation 



matrix, and r,- is the reference pixel (Calabretta & Greisen, 2002). These can 



then be transformed into one of the Sun-centred coordinate systems described 
below. 



2.4.1 Heliographic Coordinates 

Features on the Sun are located by the coordinates of latitude, 9, and longitude, 
$, with respect to the solar equator and rotational axis. In the Stonyhurst ap- 
proach, the zero point of longitudinal measurements is set at the intersection of 



the solar equator and central meridian as seen from Earth (Figure 2.12 1). In the 
Carrington approach, the central meridian is fixed to its observation on 9 Novem- 
ber 1853, and the rotations since then are counted and labelled as the Carrington 



rotation number (Carrington, 1863). 
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Figure 2.12: Schematics of two Sun-centred coordinate systems, reproduced from 



Thompson (2006). (a) Stonyhurst heliographic coordinates commonly used to lo- 
cate features on-disk. (b) Heliocentric-cartesian coordinates commonly used for 
spatially localising features in the vicinity of the Sun. 



2.4.2 Heliocentric Coordinates 

Heliocentric coordinates specify the location of a feature in space with respect 
to the centre of the Sun. The Heliocentric-Cartesian coordinate system (x, y, z) 
has the z-axis from Sun-centre along the Sun-observer line, the y-axis is per- 
pendicular to this and lies in the plane containing solar north, and the x-axis is 



perpendicular to both y and z and increases towards solar west (Figure 2.12 d). 
The Heliocentric-Radial coordinate system shares the same z-axis but measures 
features in cylindrical coordinates with radial distance p from the z-axis, and 
position angle ip counter-clockwise from solar north. 

With observations from the STEREO/SECCHI instrument suite centred on 
the Sun and extending out along the full Sun-Earth line, the optimal coordinate 
systems are the heliocentric coordinates, which can be defined in three possible 
manners depending on the user's preference. 
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2.4.2.1 Heliocentric Earth Equatorial (HEEQ) 

The HEEQ system is closely related to the Stonyhurst heliographic coordinates 
with the ^-axis parallel to the solar rotational axis, and the x-axis towards the 
intersection of the solar equator and central meridian as seen from Earth, obtained 
by the following transformations: 

Xheeq = rcos6cos<I> 

Yheeq = rcos0sin$ (2.4) 
Zheeq = rsin© 

It is thus useful when considering features in space with respect to the Sun-Earth 
line and the Sun's rotational axis. 

2.4.2.2 Heliocentric Earth Ecliptic (HEE) 

The HEE system has the x-axis towards the Earth from Sun centre, and the 
2-axis is perpendicular to the plane of Earth's orbit around the Sun called the 
ecliptic, and the y-axis perpendicular to both x and z. It is thus useful when 
considering features in space with respect to the Sun-Earth line and ecliptic. 

2.4.2.3 Heliocentric Aries Ecliptic (HAE) 

The HAE system has the x-axis towards the First Point of Aries (the direction 
to the point of intersection between Earth's equatorial plane and the plane of the 
ecliptic), the z-axis perpendicular to the ecliptic, and the y-axis perpendicular to 
both x and z. It is thus useful when considering features in space with respect 
to a fixed point on the celestial sphere (outlined below). 
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Figure 2.13: Schematic of the celestial sphere and ecliptic plane. The celestial 
equator is a projection of the Earth's equator, and the ecliptic is a projection of 
the Earth's orbit about the Sun. 



2.4.3 Helioprojective Coordinates 



When considered on a large scale it is more intuitive to take the projection of 
the Heliocentric-cartesian coordinates onto the celestial sphere; an imaginary 
sphere of arbitrarily large radius, centred on the Earth such that all observa- 



tions may be considered as projections upon it (Figure 2.13). This results in the 
Helioprojective-Cartesian coordinates (8 X , 9 y , (): 



81 



2. INSTRUMENTATION 



where d is the distance between the observer and the feature, and D Q is the 
distance between the observer and Sun centre. This similarly results in the 
Helioprojective- Radial coordinates (5 p , ip, (): 

5 P = 6 p - 90° where 6 p « (^^) ~ d ( 2 - 6 ) 

These sets of coordinates are thus useful when considering features in space with 
respect to the whole sky as seen from Earth. 
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Chapter 3 

Detecting and Tracking CMEs 



In coronagraph images CMEs are observed as outwardly moving regions of stronger 
brightness intensities relative to the background corona. Different approaches to 
thresholding the intensity of CMEs in these images have been employed in or- 
der to detect their appearance and track their motion through the field-of-view, 
leading to a cataloguing of their kinematics and morphology. However, these 
techniques suffer several drawbacks and, as such, different catalogues can vary 
significantly in their description of events. We introduce a method of multiscale 
analysis to overcome certain drawbacks of previous detection and tracking meth- 
ods. In multiscale decompositions of images, noise and small-scale features are 
removed to leave only larger-scale features of interest such as CMEs. This al- 
lows them to be tracked through the image sequences in order to determine their 
changing kinematics and morphology (see Chapter [4]). Unfortunately coronal 
streamers tend to appear on similar size scales to CMEs, making their automatic 
detection difficult. Streamers do, however, tend to remain static on timescales 
comparable to CME propagation through the field-of-view, and contain much less 
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angular information than the typically curved structure of CMEs, so they may 
be removed through spatio-temporal filtering of multiscale CME images. This 
chapter discusses the previous CME detection catalogues, and outlines our use 
of new methods of multiscale filtering to detect the CME edges in single images. 
We discuss our efforts to extend this to an automated CME detection algorithm. 
We also outline an ellipse characterisation of the CME front for study. 



3.1 CME Detection Catalogues 

Current methods of CME detection have their limitations, mostly since these 
diffuse objects have been difficult to identify using traditional image processing 
techniques. These difficulties arise from the varying nature of the CME mor- 
phology, the scattering effects and non-linear intensity profile of the surrounding 
corona, the presence of coronal streamers, and the addition of noise due to cosmic 
rays and solar energetic particles (SEPs) that impact the coronagraph detectors. 
The images are also prone to numerous instrumental effects and possible data 
dropouts. The following standard preprocessing methods are usually applied to 
optimise the images for CME studies. The coronagraph images are normalised 
with regard to exposure time in order to correct for temporal variations in the 
image statistics. A filter may be applied to remove pixel noise, for example to 
replace hot pixels with a median value of the surrounding pixel intensities, or to 
reduce the effects of background stars in the image. A correction for vignetting 
effects and/or lens distortion may be applied to the images. A background sub- 
traction may also be applied, obtained from the minimum of the daily median 
pixels across a time span of a month. The occulting disc is normally masked, along 
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Figure 3.1: Raw (left) and pre-processed image (right) of a CME observed by 
LASC0/C2 on 1 April 2004. The pre-processing includes normalising the image 
statistics, subtracting the background, and masking the occulter disk. The white 
circle (right) indicates the relative size and position of the Sun behind the occulter. 



with any data drop-outs in the images. These steps lead to a clear improvement 



in the image quality for CME study (Figure 3.1). 



3.1.1 CDAW 

The CME catalogue hosted at the Coordinated Data Analysis Workshop (CDAYVQ 
Data Center grew out of a necessity to record a simple but effective descrip- 



tion and analysis of each event observed by SOHO/LASCO (Gopalswamy et al 



2009b). The catalogue is wholly manual in its operation, with a user tracking the 
CME through C2 and C3 running-difference images and producing a height-time 
plot of each event. A linear fit to the height-time profiles provides a lst-order 
estimate for the plane-of-sky velocity, and a quadratic fit then provides a 2nd- 

1 http://cdaw. gsfc.nasa.gov/ CMEJist 
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order velocity fit and an acceleration for the event. The central position angle 
and angular width of the CME are also deduced from the images, and the event 
flagged as a halo if it spans 360°, partial halo if it spans > 120°, and wide if it 
spans > 60°. The catalogue itself lists each CME's first appearance in C2, cen- 
tral position angle, angular width, linear speed, 2nd-order speed at final height, 
2nd-order speed at 20 R Q , acceleration, mass, kinetic energy, and measurement 
position angle (the angle along which the heights of the CME are determined). 
While the human eye is supremely effective at distinguishing CMEs in coron- 
agraph images, errors may be introduced to the manual cataloguing procedure 
through the biases of different operators; for example, in deciding how the im- 
ages are scaled, where along the CME the heights are measured, or whether a 
CME is even worth including in, or discarding from, the catalogue. In an effort 
to overcome such biases, different automated catalogues have been developed to 
perform robust CME detections over large data-sets. This is also of great benefit 
for future missions where the data rate is expected to be too high for manual 
cataloguing to remain feasible. 



3.1.2 CACTus 



The Computer Aided CME Tracking catalogue (CACTusj |Robbrecht fc Bergh 



mans 



2004 ) was the first automated CME detection algorithm, in operation since 
2004. It is based upon the detection of CMEs as bright ridges in time-height slices 
(t, r) at each angle 9 around a coronagraph image. The images are preprocessed 
as standard, then a running-difference technique is applied and each image trans- 
formed into Sun-centred polar coordinates (r, 6), rebinned, and the C2 and C3 
1 ht t p : / / side . oma.be/cactus/ 
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(a) The (t, r) slice for a given angle 9, with a mirrored illustration of the resulting 
CME intensity ridge detections. 




(b) Left: the (t, r) slice for a given angle 9, with an example ridge drawn from onset time 
to with duration At across the field-of-view from r min to r max . Right: the corresponding 
accumulator space (to, At) where the ridge will appear as a point with a magnitude cor- 
responding to the ridge intensity. This modified Hough transform is used to threshold the 
most significant ridges in the slice, automatically detecting the CME in the coronagraph 
image. 



Figure 3.2: The top image (a) shows the detection of ridges in the (t, r) stacks 
of the CACTus catalogue, through the use of the Hough transform detailed in the 
bottom image (b), reproduced from Robbrecht &; Berghmans (2004). 
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fields-of-view combined. These are then stacked in time, and for each angle the 
corresponding (t, r) slice undergoes a modified Hough transform for detecting 
intensity ridges across it. This works by parameterising the (t, r) slice by the 
variables to and At, corresponding to the coordinate intersection point with the 
time axis, and the distance along the time axis respectively (together called the 



accumulator space; see Figure 3.2). So the equation of a line corresponding to an 



intensity ridge in the slice is given by: 



r max r min / , , \ , /o i N 

r = ^- {t-t )+r min (3.1) 



Thresholding the most significant ridges in the resultant accumulator space filters 
out the progression of CMEs, with the variables for each ridge characterised by 
onset time tR, the velocity vr (~ for angle Or, to give a characteristic 
variable Ir = (vr, Or, A 3D scatter plot (v, 0, t) of all detected ridges Ir 
is then integrated along the f-direction to identify clusters in the resulting (0, t) 
map which illustrates the angular span and duration time of the detected CMEs 
in the coronagraph data. A median velocity across the angular span is quoted as 
the CME speed. 

The running-difference cadence, the ridge intensity threshold, and the imposed 
limit on how many frames a CME may exist (and indeed the definition of a CME) 
all affect how successful the detection can be. However, Robbrecht &: Berghmansj 



(2004) show the algorithm to be robust in reproducing well the detections of a 
human user by direct comparison with the CDAW catalogue. The main drawback 
of the CACTus catalogue for studying CMEs is the imposed zero acceleration 
of the detection algorithm, since the Hough transform thresholds the ridges as 
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straight lines whose slopes provide a constant velocity. The velocity itself may 
also be underestimated since it is a median across the span of the CME. The 
angular spans are possibly over-estimated since side outflows in the images are 
enhanced by the running- difference and may also include streamer deflections. 
It is also difficult to distinguish when one CME has fully progressed from the 
field-of-view and another CME has entered it, so in some cases trailing portions 
of a CME are detected as separate events. 



3.1.3 SEEDS 



The Solar Eruptive Event Detection System (SEEDS^J |01medo et al. 2008) is 



an automated CME detection algorithm for tracking an intensity thresholded 
CME front in running- difference images from LASCO/C2. The images are pre- 
processed as standard, unwrapped into Sun-centred polar coordinates (r, 9), and 
a normalised running-difference technique is applied using the following equation: 



U: 



Ui - Ui-i 



a 

IT- 3.2 
At y ' 



where Ui is the running-difference image, n is the mean of the pixels in the entire 
field-of-view of the image n, At is the time difference between images (in minutes), 
a is a constant set to approximately the smallest time difference {At) between 
any image pair, where i is the current image and i — 1 the preceding image. This 
normalised difference ensures that the mean of the new image (u^ will effectively 
be zero. 

The pixel intensities (positive values only) are then summed along angles and 
1 http : / / spaceweat her . gmu . edu /seeds/ 
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Figure 3.3: Example of the SEEDS CME detection and height determination, 

(a) shows the running- difference im- 




reproduced from Olmedo et al. (2008). 
age unwrapped into Sun-centre polar coordinates, showing a CME observed by 
LASCO/C2 on 12 September 2002. The black line distribution across the image 
represents the positive value intensity count along each angle, and the two verti- 
cal black lines mark the angular span at one standard deviation above the mean 
intensity, (b) shows the new angular span following the region growing technique, 
(c) shows the intensity within the angular span averaged across heights, and the 
'Half-Max-Lead' is taken as the CME height in the image. 
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thresholded at a certain number of standard deviations above the mean intensity: 



fi + Na (cf. Equation 3.12) as in Figure 3.3 1. This determines the 'core angles 



of the CME, and a region growing technique based on a secondary threshold 
of intensities in the rest of the image is applied to open the angular span to 



include the full CME (Figure 3.3 d). Issues arise when streamer deflections occur 
that will offset the region growing technique and overestimate the CME angular 
width. An intensity average across the angles within the span of the CME is 
then determined, and where the forward portion of this intensity profile equals 



half its maximum value is taken as the CME height (Figure 3.3:). The velocity 
and acceleration are determined from the heights through consecutive images and 
these results are output with the CME position angle and angular width in the 
SEEDS catalogue. 

Along with the issues of streamer deflections and the tracking being limited 
to C2 images, the choice of the 'Half-Max-Lead' as the CME height is dependant 
on the overall CME brightness, and thus any brightness changes as the CME 
propagates will affect this measurement. This would add to the error on the 
height-time profile which, along with the error in time as a result of the running- 
difference technique, makes it difficult to accurately determine the velocity and 
acceleration. 

3.1.4 ARTEMIS 

The Automatic Recognition of Transient Events and Marseille Inventory from 
Synoptic maps (ARTEMISQ |Boursier et olj |2009b[ ) is an automated CME detec- 



tion algorithm that works by identifying signatures of transients in synoptic maps. 
1 http : / / www . oamp . fr /lasco / 
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These maps are generated as (t, 9) slices for specific heights r in the coronagraph 
images of LASC0/C2. The images are prepared through the standard prepro- 
cessing steps. Then at a specific height (e.g., r = 3 R ) the intensity is plotted 
across all angles 9 for each image through time t with transient events appearing 



as vertical streaks through the more persistent streamer intensities (Figure 3.4). 
A method of image filtering and intensity thresholding is applied to distinguish 
the streaks in the synoptic map, and image segmentation then discards small 
features and closes off regions-of-interest (ROIs) to produce a binary map of the 
streaks. Specific parameters of these streaks are also computed, such as their 
total radiances, areas and centres of gravity. Merging with high-level knowledge 
helps to associate ROIs of the same CME if three criteria are met: the difference 
between the ^-coordinates of the centres of two ROIs differs by less than two pix- 
els; the difference between the ^-coordinates of their centres differs by less than 
60 pixels (corresponding to a 60° angular span); and the ratio of their radiances 
calculated at their centres (on the original synoptic map) ranges from 0.25 to 4. 
The result is a binary CME detection map in (t, 9) space for different heights in 
the corona: 3, 3.5, 4, 4.5, 5 and 5.5 R Q . 

With the CME detections in place, estimates of the velocity may be made. A 
first estimate is taken by testing a range of constant velocities 50-2,000 km s" 1 
to determine which best matches the shifting of the CME detection in synoptic 
maps at subsequent heights through the corona. The binary maps are shifted 
by an amount corresponding to velocity steps of 10 km s _1 , such that the one 
which provides the maximum pixel value (with a minimum limit of 3) indicates 
the best velocity estimate of the event. A second estimate is taken by cross- 
correlating the detected CME ROIs on the original synoptic maps at 3 R and 
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Figure 3.4: An example of how the synoptic maps are generated for the ARTEMIS 
catalogue, reproduced from Boursier et al.\ ( 2009b| ). At a chosen height in the 
coronagraph image an annulus is unwrapped (indicated with the dashed line and 
blue square, circle and triangle) and these are then stacked together to illustrate 
how the intensity at that height changes through time. Vertical streaks represent 
transient events occurring on smaller time-scales than the more persistent streamers 
in the images. 
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5.5 R Q and inspecting the intensity shift in time (pixel shift in x-direction) to ob- 
tain the velocity estimate. A third estimate is taken by similar cross-correlation 
but specifically on each individual line of the ROIs to obtain a distribution of ve- 
locities across the angular span of the CME, the median of which is taken as the 



actual velocity. Boursier et al. (2009b) compare histograms of the three different 
velocity estimates for the ARTEMIS CME detections over a twelve year interval 
and find that, globally, the three estimates are highly consistent with each other. 

ARTEMIS is limited to the C2 field-of-view and it provides kinematics only 
in the 3-5.5 R range. The velocity determinations themselves are not specific 
to either the CME front nor any other identifiable feature, and carry all the in- 
accuracies resulting from the image rebinning, intensity averaging, filtering and 
segmentation techniques in generating the final detection masks. 



Due to the drawbacks of each of the catalogues above, the motivation exists 
to study the kinematics and morphology of CMEs with as great an accuracy as 
possible in order to better compare with theory. To this end we outline below 
our application of multiscale analysis to remove small scale noise/features and 
enhance the larger scale CME in single coronagraph frames, allowing the CME 
front edges to be detected and a geometrical characterisation applied to study its 
propagation with increased accuracy for deriving the kinematics and morphology. 



3.2 Multiscale Filtering 

In this section a new multiscale method of analysing CMEs is described. The 
use of multiscale methods in astrophysics have proven effective at denoising spec- 
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tra and images (Fligge & Solanki, 1997 Murtagh et al, 1995), analysing solar 



active region evolution (Hewett et al., 2008), and enhancing solar coronal im- 



ages (Stenborg & Cobelli, 2003 Stenborg et al, 2008). A particular application 



of multiscale decompositions uses high and low pass filters convolved with the 



image data to exploit the multiscale nature of the CME (Young & Gallagher 



2008). This highlights its intensity against the background corona as it propa- 



gates through the field-of-view, while neglecting small scale features (essentially 
denoising the data). It also leads to the use of non-maxima suppression to trace 



the edges in the CME images, and Young & Gallagher (2008) show the power of 



multiscale methods over previous edge detectors such as Roberts (Davis, 1975) 



and Sobel (Duda & Hart, 1973). With these methods for defining the front of 



the CME we can characterise its kinematics (position, velocity, acceleration) and 
morphology (width, orientation) in coronagraph images. Multiscale analysis also 
has the benefit of working on independent images without any need for differenc- 
ing, so the temporal errors involved are on the order of the exposure time of the 
instrument (~ a few seconds). 

The fundamental idea behind wavelet analysis is to highlight details apparent 
on different scales within the data. An example of this is the suppression of noise 
in images, which tends to occur only on the smallest scales. Wavelets have bene- 
fits over previous methods (e.g., Fourier transforms) because they are localised in 
space and are easily dilated and translated in order to operate on multiple scales, 
the basic equation being: 



1 , { t-a 



(3.3) 
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where a and b represent the shifting (translation) and scaling (dilation) of the 
mother wavelet i/j which can take several forms depending on the required use. 

We explore a method of multiscale decomposition in 2D through the use of 
low and high pass filters; using a discrete approximation of a Gaussian, 9, and its 



derivative, ip, respectively (Stenborg & Cobelli, 2003). Since 9(x,y) is separable 



i.e. 9(x,y) = 9(x)9(y), we can write the wavelets as the first derivative of the 
smoothing function: 



_ 2 d9(s l x) _j 
dx 



-9{s- l y) 



l 9(s- l x 



d9{s~ 1 y) 
dy 



(3.4) 
(3.5) 



where s is the dyadic scale factor such that s = 2 J where j = 1, 2, 3, J. Succes- 
sive convolutions of an image with the filters produces the scales of decomposition, 
with the high-pass filtering providing the wavelet transform of image I(x, y) in 
each direction: 



W S X I = W:i(x,y) = P x {x,y)*I(x,y) 



W s y l ee W s y I{x,y) = P v (x,y)*I{x,y) 



(3.6) 
(3.7) 



Akin to a Canny edge detector (Young & Gallagher, 2008), these horizontal and 



vertical wavelet coefficients are combined to form the gradient space, V s (x, y), for 
each scale: 

T s (x,y) = [W°I, Wyl] (3.8) 
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Figure 3.5: Top left, the horizontal detail, and top right, the vertical detail 
from the high-pass filtering at one scale of the multiscale decomposition (called 
the rows and columns respectively). Bottom left, the corresponding magnitude 
(edge strength) and bottom right, the angle information (0 - 360°) taken from the 



gradient space, for a CME observed in LASCO/C2 on 1 April 2004 (Byrne et al. 
20091). 
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Figure 3.6: The vectors plotted represent the magnitude and angle determined 
from the gradient space of the high-pass filtering at a particular scale. The CME 



of 2004 April 1 shown here is highlighted very effectively by this method (Byrne 



et al, 2009) 



The gradient information has an angular component a and a magnitude (edge 
strength) M: 



a s (x,y) = tan- 1 (W°I / W s x l) 



M s (x, y ) = \imiy + {w°iy 



(3.9) 
(3.10) 



The resultant horizontal and vertical detail coefficients, and the magnitude and 



angular information are illustrated in Figure 3.5 



At a particular scale the signal-to-noise ratio of the CME is highest and this 
is the optimum scale for determining the edges in the image. The angular com- 
ponent a of the gradient specifies a direction which points across the greatest 
intensity change in the data (an edge). A threshold is specified with regard to 
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this gradient direction in order to chain pixels along maxima to highlight the 
edges. The changes in magnitude and angular information may then be imple- 
mented in a spatio-temporal filter for distinguishing those edges corresponding 
to the CME only. Overlaying a mesh of vector arrows on the data shows how 
the combined magnitude and angular information illustrate the progression of the 
CME. Each vector is rooted on a pixel in the gradient space, and has a length 



corresponding to the magnitude M with an angle from the normal a (Figure 3.6). 
Using this information, it becomes possible to create a specific detection mask 
which is used to identify the edges along the CME front to study its propagation 
(as done for a sample of events in Chapter [4]). However, for the cases of faint 
CMEs or strong streamer deflections, the filter is presently limited by exploiting 
the information from only one scale and ignoring all other scales, meaning it cur- 
rently often requires the user to remove/include certain edges that the algorithm 
has mistakenly retained/discarded. Extending the algorithm to work on more 
than one scale may help alleviate this issue in order to develop a fully automated 
CME detection and characterisation routine, as outlined in the following section. 



3.3 Automated Multiscale CME Detection 



For the most part, CMEs exist on size scales larger than noise and any small 
scale features in coronagraph images, such as stars or planets, that are redundant 
for studying CME propagation. This fact has led to the development and imple- 
mentation of multiscale decompositions that highlight the CME in images from 



SOHO/LASCO and STEREO/SECCHI (Byrne et al. 2009 Stenborg fe Cobelli 



2003 Young & Gallagher, 2008). However, coronal streamers (plasma outflows 
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from open magnetic field regions on the Sun) can persist through coronagraph 
images with significant brightness intensities and tend to appear on similar scales 
as the CME in multiscale image analysis. If a CME propagates through an im- 
age with a strong streamer present, it becomes difficult to distinguish the two 
features by intensity thresholding alone, and this is one reason why differencing 
techniques have been widely used in CME analysis. In an effort to move away 
from differencing and the large errors involved in the subtraction of images from 
each other, since the goal is to obtain kinematics with the greatest precision, we 
endeavour to separate CME and streamer features from one another using mul- 
tiscale methods alone. These efforts involve exploiting the angular distribution 
that exists across a curved CME front compared to the more linear streamers 
in single independent images. To do this, the coronagraph images must first be 
normalised for their radial gradient in intensity, since the drop-off across the field- 
of-view is too steep to effectively segment a single entire streamer from the inner 
to outer edge of an image. This can serve to enhance the noise at the edge of the 
images, but this is again suppressed by the multiscale analysis. Occasions when 
the CME propagates directly in front of, or behind, a streamer remain problem- 
atic, as do strong streamer deflections that can occur when a CME propagates 
into or expands alongside a streamer. 



3.3.1 Normalising Radial Graded Filter (NRGF) 



Since the brightness drop-off of the corona is large, falling from approximately 
10~ 6 -10~ 9 B across heights of 1-6 R (Kimura & Mann, 1998), a method 



for radially normalising coronagraph images to enhance features across this steep 
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Figure 3.7: A normalised, background subtracted, LASC0/C2 image (left) of a 
CME on 1 April 2004, and the resulting NRGF image (right). The image radial 
intensity is scaled such that structure along streamers and the CME becomes visible 
across the field-of-view 



intensity gradient was developed by Morgan et al. (2006). It works by normalising 



the intensity in radial coordinates of the image according to the equation 



I'(r, 



I(r, 



Ur, 



<4» 



(3-11) 



where I'(r, 0) is the processed and I(r, <fi) is the original intensity at height r and 
position angle 0, and J(r)<^> and cr(r) < ^ )> are the mean and standard deviation 



of intensities calculated over all position angles at height r. Figure 3.7 shows 
the result when the NRGF is applied to a LASCO/C2 image of the 1 April 2004 
CME. 



The multiscale decomposition introduced in Young & Gallagher (2008) pro- 



vides magnitude and angular information of the edges in the image. This in- 
formation is combined to chain the strongest edges within the image on a scale 
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that provides an optimal signal-to-noise ratio for studying the CME. |Byrne et al. 



(2009) obtain the CME front edges in this manner, which are then used to fit 
an ellipse to characterise the CME propagation in the image sequence. In order 
to automate the algorithm, thresholds on the magnitude information (e.g., CME 
edges appear on larger scales than noisy features) and angular information (e.g., 
CME edges appear more curved than streamer edges) were investigated. The 
thresholding is strengthened by the inclusion of more than one scale in localising 
the CME and distinguishing it from the streamers, detailed below. 



3.3.2 Thresholding 

The magnitude information corresponds to the strength of the edges in the image, 
and so can be thresholded to discard the small scale noise. For the NRGF image 



(Figure 3.7), a hard threshold T is set at one standard deviation a of the mean \i 
of the image intensity to contour regions of interest that may be a CME according 
to the equations: 



1 n 

T = u + Na = - V Xi +N 



\ 



1 U 

(3.12) 



n 
i=i 



mean / ~, 7^ . ,. 

standard deviation 



where Xj are the pixel intensity values of the image, and we choose the number of 



standard deviations N = 1. The left image of Figure 3.8 illustrates this thresh- 
olding with a sample of contoured regions (outlined in red) on the multiscale 
decomposition of a CME observed in LASCO/C2 on 1 April 2004. 

It is apparent from the right of Figure |3.8| that the CME will contain edges 
whose normals are widely distributed across - 360° compared to the more linear, 
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Figure 3.8: A chosen scale of the decomposed NRGF image provides a magnitude 
image of the edge strengths displayed on the left, which is thresholded at one 
standard deviation from the mean intensity to obtain contoured regions of interest 
that could contain a CME (sample contours indicated in red). As shown, the 
streamers have edges which appear on the same scale as the CME edges in this 
image. The angular information from the decomposition is displayed on the right, 
and the contoured regions of interest overlaid for comparison. The grey scale 
indicates angles from 0-360° and it is clear that streamers tend to have a linear 
grey scale while the CME has a gradient of greys across the scale. 



radially directed, streamer edges. The angular distribution of each region is de- 
termined and then normalised and folded into 0-180° range, centred on 90°, due 
to the symmetry of the edge normals. This is illustrated for four selected contour 



regions in Figure |3.9[ The resulting angular distributions are then thresholded 
with regard to their median value, since the distribution of angles across the CME 
will be wider and have a higher median value than for a distribution of angles 
along a streamer. 

This thresholding is repeated across four scales of the multiscale decompo- 
sition, neglecting smaller scales dominated by noise and larger scales that are 
overly smoothed. Assigning a score to the regions that may contain a CME at 
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Medion 




0.01 1 




i 





Figure 3.9: Left: four contoured regions (at one standard deviation of the mean 
image intensity) highlighted on the magnitude information from the multiscale 
decomposition of the 1 April 2004 CME. Right: the corresponding angular distri- 
bution of each region, normalised and folded into the 0-180° range (centred on 
90°). The angular distribution may be thresholded with respect to its median value 
to distinguish regions corresponding to CMEs from those along streamers. 



each scale, it becomes possible to build a detection mask as in Figure [S\10 

The scoring system is chosen arbitrarily to work best with the chosen thresh- 
olds, and these may be changed and refined as an analysis of more CMEs is done. 
For example, the current thresholds from working on a sample of ~ 10 CMEs are 
as follows: 



1. The magnitude information is thresholded at one standard deviation (la) 
of the mean intensity across the image. 

2. The 15 largest contoured regions across the image are investigated (there 
are rarely more than ~ 5 streamers of similar intensity to a CME, and we 
allow for disjointed contours along structures). 



3. If the median angular value is > 20% of the angular distribution peak then 
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Figure 3.10: The resulting CME detection mask from combining the thresholded 
regions of strongest magnitude and angular distribution at four scales of the de- 
composition. The 3D representation on the right illustrates the pixel values of the 
mask. 



the region is deemed a CME and assigned a score of 3 (the pixels in that 
region of the mask are given the value 3). If it is > 10% the score is 2 
(potential CME structure), or > 5% the score is 1 (weak CME structure or 
portion thereof). 

4. The final CME mask through the combination of scores at each scale results 
in a dominant region that localises the CME front in the image and can 
be used to characterise the front, or input into a spatio-temporal filter if 
subsequent CME images are available in order to refine the masked region 
whenever streamers are still present. 

The resultant set of detection masks for three of the frames of the CME observed 



on 1 April 2004 are shown in Figure 3.11 
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Figure 3.11: The NRGF (left) and resulting detection masks (middle and right) 
for different frames of a CME observed by LASC0/C2 at times of 00:00 UT (top), 
00:40 UT (middle) and 01:20 UT (bottom) on 2 April 2004. The location of the 
CME front is highlighted very efficiently by this method, although the detection 
masks may contain artefacts of the chosen thresholds which must be discarded 
when characterising the CME front. 
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3.3.3 Faint CMEs and Streamer Interactions/Deflections 



Due to the nature of the hard thresholds in place on the magnitude and angu- 
lar information, there are problems which arise when the algorithm mistakenly 
disregards a CME or includes a streamer, or portions thereof. Firstly, if a CME 
is faint enough that the intensity falls below the la magnitude threshold, it will 
not be detected as a region of interest in the image. Secondly, if the CME inter- 
acts with a streamer, the two features may be contoured together and this will 
skew the angular distribution and affect the detection mask. And thirdly, if the 
CME causes a significant streamer deflection, it will lead to a wider distribution 
of angles along the streamer and the algorithm may thus detect it as part of the 
CME. This is why the above scoring system was introduced in an effort to min- 



imise these effects, which are highlighted in Figure 3.12 for a CME observed on 
23 April 2001. The event is too faint compared with the streamers across it for it 
to be easily distinguished in the image, and parts of the streamers are then mis- 
takenly included in the final detection mask. This is where the current algorithm 
requires a user to specify which parts of the edges correspond to the CME for 
characterisation. Such limitations in current wavelet analysis of CMEs may be 
overcome by extending these algorithms to work with ridgelets or curvelets that 



better suit the curved form of a typical CME front as discussed in Gallagher et al. 



(2010). Furthermore, it should be noted that a relatively small number of CMEs 
exhibit a narrow, spiky or strongly kinked structure which may not be readily 
distinguished from streamers in an automated fashion, nor have an adequately 
high angular distribution for automatic detection in the images. These CMEs 
may also not be satisfactorily characterised with an ellipse, discussed in the fol- 
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Figure 3.12: Example of the difficulty in detecting the faint CME observed by 
LASCO on 23 April 2001. If the CME is far-sided and/or of low intensity it becomes 
difficult to threshold its edges in the image compared to the coronal streamers, (a) 
is the pre-processed CME image, (b) is the NRGF image, (c) and (d) show 
the magnitude and angular information from the multiscale decomposition with 
intensity contours overlaid in red. (e) and (f) show the resulting CME detection 
mask in 2D and 3D respectively. 
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lowing section, however they are not typical of the commonly observed curved 
nature of CMEs and thus represent a small class of events outside the limits of 
these methods as they currently stand. 



3.4 Characterising the CME Front 



Using a model such as an ellipse to characterise the CME front across a sequence 
of images, has the benefit of providing the kinematics and morphology of a mov- 
ing and/or expanding structure. The ellipse's multiple parameters, namely its 
changeable axis lengths and tilt angle, is adequate for approximating the vary- 



ing curved structures of CMEs. Chen et al. (1997) suggest an ellipse to be the 



two-dimensional projection of a flux rope, and Krall & St. Cyr (2006) use ellipses 
to parameterise CMEs and explore their geometrical properties. It also serves as 
the observed projection of the base of the cone model applied to CME images 



Xie et al. 


2004 


Xue et al. 


2005 


Zhao et al. 


2002 



determined to be along the CME front by means of a Levenberg-Marquardt least 
squares algorithm. A kinematic analysis then provides height, velocity and accel- 
eration profiles; while the ellipse's changing morphology provides the tilt angle 



and angular width (see the example in Figure 3.14). Measuring these properties 
in the observed data is vitally important for accurate comparison with theoretical 
models. 



Following Schrank (1961), we may determine the polar equation of a tilted 



ellipse by starting with the standard equation for an ellipse with centre point 
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Figure 3.13: Ellipse tilted at angle 5 in the plane, with semimajor axis a, semimi- 
nor axis b, and radial line p inclined at angle oj to the semimajor axis. 



(xo,yo), semimajor axis a, and semiminor axis b: 



(g - ^o) 2 , (y - yo) 2 = 1 

a 2 b 2 



(3.13) 



This is written in polar coordinates by x = pcosu, y = psinu and centred on 
the origin (x = 0, y = 0) to give: 



p 2 cos 2 co p 2 sin 2 oj 
5 + = 1 



b 2 



(3.14) 



where p is a radial line from the centre to any point on the ellipse, at an angle u 



to the semimajor axis a (Figure 3.13). Allowing for a tilt angle S on the ellipse, 



we may define oj' = u + 5 to obtain: 



21,2 



a 2 b 



(^f) - cos(2w' - 26) 



(3.15) 
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This gives a first approximation which can then be used to iteratively solve the 
ellipse parameters until a best fit to the points along the CME front is obtained 
(Figure [314]. 
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Chapter 4 

The Kinematics and Morphology 
of CMEs using Multiscale 
Methods 

The diffuse morphology and transient nature of CMEs make them difficult to 
identify and track using traditional image processing techniques. We apply mul- 
tiscale methods to enhance the visibility of the faint CME front. This enables an 
ellipse characterisation to objectively study the changing morphology and kine- 
matics of a sample of events imaged by SOHO/LASCO and STEREO/SECCHI. 
The accuracy of these methods allows us to test the CMEs for non-constant ac- 
celeration and expansion. This chapter is founded on work published in Byrne et 
al, Astronomy & Astrophysics (2009). 
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4.1 Introduction 



To date, most CME kinematics are derived from difference images; a technique 
based either on the subtraction of a single pre-event image (fixed-difference) or the 
subtraction of each image from the next in an event sequence (running-difference). 
These techniques are applied in order to highlight regions of changing intensity, 
increasing the relative brightness of the CME against the background coronal fea- 
tures. However, drawbacks do exist. Numerical differencing can enhance noise to 
a level comparable to the signal. The noise can be suppressed to a certain degree 
by using a standard box-car or median filter, but this will also smooth out CME 
features such as structure along the CME front and its environs. An additional 
issue resulting from differencing is the introduction of spatio-temporal cross-talk 
in difference frames. Since it is used to highlight non-stationary features in both 
space and time, then the differencing of subsequent images of a moving feature 
will show a signature at the position where the feature was initially observed and 
a signature at the position that the feature has moved to when next observed. 
Since the signature of motion in the difference images is heavily dependent on 
the time between frames and how many pixels the feature has moved, it may be 
considered to blend spatial and temporal information in a non-trivial manner: 
an effect referred to as spatio-temporal cross-talk. This can serve to blur out 
CME features and introduce ambiguity in estimating positions and times, criti- 
cal to accurately deriving the kinematics of the event. Furthermore, user bias is 
introduced by the choice of intensity scaling and thresholding when determining 
the location of CME features by point-and-click methods or automated detection 



algorithms, as discussed in Section 3.1 



114 



4.1 Introduction 



In this work we apply multiscale methods for analysing CMEs as described in 
Chapter [3j which has the benefit of working on independent images without any 
need for differencing. Once the edges of the CME front are resolved, an ellipse 
characterisation is applied to determine the CME kinematics (position, velocity, 
acceleration) and morphology (width, orientation) in coronagraph images. CME 
height measurements are taken as the height of the furthest point on the ellipse 
from Sun centre. The angular width is taken as the opening angle of the ellipse 
from Sun centre, and the tilt of the ellipse is given by the calculated angle 5. 
(Note that in cases where the code produces an extremely large and oblate ellipse 
with one apex approximating the CME front, the width and tilt information is 
deemed redundant. Hence the resulting analysis of some events can have less 
data points included in the width and tilt plots than in the height-time plots.) 



Following previous concerns on the errors in CME heights (e.g., Wen et al. , 2007), 
multiscale methods allow us to determine the kinematics to a high degree of ac- 
curacy in order to improve confidence in their interpretation and comparison to 
theory. These methods also show potential for future automation. In this chap- 
ter a sample of CMEs observed by SOHO/LASCO and STEREO/SECCHI are 
studied, namely the gradual events of 2 January 2000, 18 April 2000, 23 April 
2001, 1 April 2004, 8 October 2007 and 16 November 2007, and the impulsive 



events of 23 April 2000 and 21 April 2002 (see Figures 4.1 and 4.2). We compare 
our results with the catalogues of CDAW, CACTus and SEEDS (note ARTEMIS 
was not included due to difficulties interpreting the entries in its database that 
correspond to the chosen CMEs). 
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4.1 Introduction 




4. THE KINEMATICS AND MORPHOLOGY OF CMES USING 
MULTISCALE METHODS 



4.2 Error Analysis 

The front of the CME is determined through the multiscale decomposition and 
consequent rendering of a gradient magnitude space. At scale 3 of the decompo- 
sition the smoothing filter is 2 3 pixels wide, which we use as our 3<j error estimate 
in edge position. This error is input to the ellipse fitting algorithm for weight- 
ing the ellipse parameters, and a final error output is produced for each ellipse 
fit. In the case of a fading leading edge the reduced amount of points along the 
front will increase the error on our analysis. The final errors are displayed in 
the height-time plots of the CMEs, and are used in the velocity and acceleration 
calculations. The derivative is a 3-point Lagrangian interpolation, so there is an 
enhancement of error at the edges of the data sets as explained below. It should 
also be noted that in the 3-point Lagrangian counter-intuitively results in larger 
errorbars on the velocity and acceleration profiles in cases where the number of 
height-time measurements are increased (i.e., for smaller cadences such as the 
inner coronagraphs) . This is due to the algorithm's inverse dependence on the 
spacing between points. 

4.2.1 3-Point Lagrangian Interpolation 

3-point Lagrangian interpolation is used on the discrete set of given data points in 
order to determine the first and second derivatives corresponding to the velocity 
and acceleration of the CME height-time measurements in a more robust manner 
than simple forward, reverse or centre difference techniques. Considering three 
data points, (xo,yo), (xi,yi), (^2,2/2), the Lagrangian interpolation polynomial is 
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given by: 



L ( x ) = ^vM^) where l i( x ) = n 

3=0 i=0,ift 1 



L(x) = VqIq{x) +y 1 l 1 (x) + y 2 l 2 (x) 

fx — Xi x — x 2 \ fx — Xo X — X2 \ (X — Xo X — X\ 

= 2/o +2/2 



^Xq - Xi Xq - X 2 J \X 1 - XqXx - X 2 J \X 2 - X Q X 2 - Xi 

So the derivative is determined to be: 



v - °gL (4,) 

2a; — rri — x 2 2x — x — x 2 2x — x Q — x\ 

= 2/0 7 w x+2/17 w v + 2/2- 



(x - xi) (x - x 2 ) (xi - x ) (xi - x 2 ) (x 2 - x ) (x 2 - Xi) 
And the edge point x = xo (and similarly for x = x n ) is weighted as follows: 



2x -x 1 -x 2 x -x 2 Xo-Xi 

* = »(x„ - x l} (x, - „) +y \ Xl - x„) {Xl - x 2) + *(x„ - „) (ll - X2) (4 - 3) 



In the case where the points are equally spaced this is simply: 



d o = I [-3i/o + 4yi - y 2 \ (4.4) 
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The error propagation equation is used to determine the errors on the resulting 
derivative points in L' = f(L(x),x): 

al , = al (M)\ al (?E)\... ( , 5) 



dx 2 dx 2 \dx 
Or more appropriately written in this context as 



dL J x \dx J 

(4. 6) 



ad = dx~ 2 + dx 2 ldx- ] (47) 



So the errors on the end points become: 

o 2 9 < + 16 <+< <+< f 3y - 4 Vl + y 2 \ 2 
d ° (x 2 -x ) 2 (x 2 -x ) 2 \ x 2 -x J 

a d n = ; r 2 + t 72 " — ; ( 4 - 9 ) 

[X n — X n ^ 2 ) \Xn-2 ~ x n) V x n-2 A n / 

This effect is reflected in the larger errorbars on the end points of the derived 
kinematics of Section 14.31 

The errors on the heights are used to constrain the best fit to a constant 
acceleration model of the form: 



h(t) = at 2 + v t + h (4.10) 

where t is time and a, vq and ho are the acceleration, initial velocity and initial 
height respectively. This provides a linear fit to the derived velocity points and a 
constant fit to the acceleration. An important point to note is the small time error 
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(taken to be the image exposure time of the coronagraph data) since the analysis 
is performed upon the observed data frames individually. Previous methods of 
temporal-differencing would increase this time error. With these more accurate 
measurements we are better able to determine the velocity and acceleration errors, 
leading to improved constraints upon the data and providing greater confidence 
in comparing to theoretical models. 



4.3 Results 

This section outlines events which have been analysed using our multiscale meth- 
ods. We use data from the LASCO/C2 and C3, and SECCHI/COR1 and COR2 



instruments, and preprocess the images as discussed in Section 3.1 Events were 
selected on the basis that they were strong candidates to test for non-constant 
acceleration and expansion due to their high signal-to-noise ratio out through 
the full coronagraph fields-of-view. The ellipse fitting algorithm applied to each 
event gives consistent heights of the CME front measured from Sun centre to 
the maximum height on the ellipse, and these lead to velocity and acceleration 
profiles of our events. The ellipse fitting also provides the angular widths and 
orientations, as shown below. The velocity (where the ranges are from initial to 
final speeds), acceleration and angular width results of each method are high- 



lighted in Tables 4.1, 4.2 and 4.3 In each instance we include the values from 
CACTus, CDAW and SEEDS. Note that CACTus lists a median speed of the 
CME; CDAW provide the speed at the final height and from the velocity profile 
we infer the speed at the initial height; and the SEEDS detection applies only 
to the LASCO/C2 field-of-view but doesn't currently provide a velocity range or 
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Table 4.1: Summary of CME velocities as measured by CACTus, CDAW, SEEDS 
and multiscale methods. Asterisks indicate analysis of SECCHI data rather than 
LASCO. 



Date 


CACTus 


CDAW 


SEEDS 


Multiscale 




km s _1 


km 


s- 1 


km s _1 


km 


s- 1 


02 Jan 2000 


512 


370 - 


- 794 


396 


396 - 


- 725 


18 Apr 2000 


463 


410 - 


- 923 


339 


324 - 


1049 


23 Apr 2000 


1041 


1490 


- 898 


595 


1131 - 


- 1083 


23 Apr 2001 


459 


540 - 


- 519 


501 


581 - 


-466 


21 Apr 2002 


1103 


2400 - 


- 2388 


702 


2195 - 


- 2412 


01 Apr 2004 


487 


300 - 


- 613 


319 


415 - 


- 570 


08 Oct 2007 


235* 


85 - 


331 


103 


71 - 


330* 


16 Nov 2007 


337* 


210 - 


- 437 


154 


131 - 


483* 



profile. Note also that the CMEs of 8 October 2007 and 16 November 2007 are 
analysed in SECCHI images by CACTus and our multiscale methods (marked by 
asterisks in the Tables), while CDAW and SEEDS currently only provide LASCO 
analysis, and so may be observing a different part of the CME structure so it is 
difficult to compare directly. Overall, it is clear that many of the CACTus, CDAW 
and SEEDS results lie outside the results of our analysis. Notably, since CDAW 
involves labour intensive analysis by qualified scientists, this is a prime indication 
that robust methods of determining the kinematics and morphology of CMEs are 
needed, which our methods of multiscale analysis and characterisation provide, 
especially given our rigourous error analysis. 



4.3.1 Arcade Eruption: 2 January 2000 

This CME was first observed in the south-west at 06:06 UT on 2 January 2000 
and appears to be a far-side event associated with an arcade eruption consisting 
of one or more bright loops. 
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Table 4.2: Summary of CME accelerations as measured by CACTus, CDAW, 
SEEDS and multiscale methods. Asterisks indicate analysis of SECCHI data rather 
than LASCO. 



Date 


CACTus 


CDAW 


SEEDS 


Multiscale 




m s -2 


m s -2 


m s -2 


m s~ 2 


02 Jan 2000 





21.3 


-5.8 


14.7 ± 3.6 


18 Apr 2000 





23.1 


17.5 


32.3 ± 3.5 


23 Apr 2000 





-48.5 


-8.9 


-4.8 ± 20.6 


23 Apr 2001 





-0.7 


-1.4 


-4.8 ± 4.1 


21 Apr 2002 





-1.4 


33.5 


32.5 ± 26.6 


01 Apr 2004 





7.1 


12.9 


4.4 ± 2.0 


08 Oct 2007 


0* 


3.4 


2.4 


5.7 ± 0.9* 


16 Nov 2007 


0* 


4.9 


11.0 


13.7 ± 1.7* 



Table 4.3: Summary of CME angular widths as measured by CACTus, CDAW, 
SEEDS and multiscale methods. Asterisks indicate analysis of SECCHI data rather 
than LASCO. 



Date 


CACTus 
degrees 


CDAW 
degrees 


SEEDS 

degrees 


Multiscale 

degrees 


02 Jan 2000 


160 


107 


96 


50 - 


- 95 


18 Apr 2000 


106 


105 


108 


68 - 


110 


23 Apr 2000 


352 


360 


130 


96 - 


130 


23 Apr 2001 


124 


91 


74 


55 - 


- 60 


21 Apr 2002 


352 


360 


186 


53 


- 65 


01 Apr 2004 


66 


79 


58 


44- 


- 38 


08 Oct 2007 


52* 


82 


59 


23 - 


60* 


16 Nov 2007 


68* 


78 


54 


40 - 


55* 



The height-time plot has a trend not unlike that of CDAW (overplotted in 



top Figure 4.3 with a dashed line). However, the offset of the CDAW heights - 
which puts them outside our error bounds - may be due to how the difference 
images are scaled for display. This is a problem multiscale methods avoid. From 



Figure 4J3 the velocity fit was found to be increasing from 396 to 725 km s 1 , 
giving an acceleration of 14.7 ±3.6 m s -2 . The ellipse fit spans approximately 
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50-70° of the field-of-view in the inner portion of C2, and expands to over 95° 
in C3. This expansion may simply be attributed to the inclusion of one or more 
loops in the ellipse fit as the arcade traverses the LASCO/C2 and C3 fields-of- 
view. The orientation of the ellipse as a function of time is shown in the bottom 



of Figure 4.3 It can be seen that the orientation angle of the CME increases to 
approximately 100° before decreasing toward 60°. 

The constant acceleration model is not a sufficient fit to the data in this event. 
The kinematics produced from the multiscale edge detection would be better fit 
with a non-linear velocity and a non-constant acceleration. This would show the 
CME to have a period of decreasing acceleration in the C2 field-of-view, levelling 
off to zero in C3 (if not decelerating further). 



4.3.2 Gradual/Expanding CME: 18 April 2000 

This CME was first observed off the south limb at 16:06 UT on 18 April 2000 
and exhibits a flux rope type structure. 

The height-time plot for this event has a trend similar to that of CDAW 



(overplotted in top Figure 4.4 with a dashed line). The velocity fit was found to 
be linearly increasing from 324 to over 1,000 km s _1 , giving an acceleration of 
32.3 ±3.5 m s~ 2 . The ellipse fit spans from 68° of the field-of-view in the inner 
portion of C2, to approximately 110° in C3. The orientation of the ellipse as a 
function of time is shown to increase from just above 0° to over 60° in Figure |4~4"1 
This CME is fitted well with the constant acceleration model but shows an 
increasing angular width, implying a type of over-expansion across the field-of- 
view. 
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Figure 4.3: Kinematic and morphological profiles for the ellipse fit to the multi- 
scale edge detection of the 2 January 2000 CME observed by LASCO/C2 and C3. 
The plots from top to bottom are height, velocity, angular width, and ellipse tilt. 
The CDAW heights are over-plotted with a dashed line. The height and velocity 
fits are based upon the constant acceleration model. 



125 



4. THE KINEMATICS AND MORPHOLOGY OF CMES USING 
MULTISCALE METHODS 




Figure 4.4: Kinematic and morphological profiles for the ellipse fit to the multi- 
scale edge detection of the 18 April 2000 CME observed by LASCO/C2 and C3. 
The plots from top to bottom are height, velocity, angular width, and ellipse tilt. 
The CDAW heights are over-plotted with a dashed line. The height and velocity 
fits are based upon the constant acceleration model. 
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4.3.3 Impulsive CME: 23 April 2000 

This impulsive CME was first observed in the west at 12:54 UT on 23 April 2000 
and exhibits strong streamer deflection. 

The height-time plot derived using our methods has a trend which diverges 



from that of CDAW (overplotted in top Figure 4.5 with a dashed line). The 
velocity fit was found to be linearly decreasing from 1,131 to 1,083 km s -1 , giving 
a constant deceleration of —4.8 ±20.6 m s -2 . The CME is present for one frame 
in C2 with an ellipse fit spanning 96°, increasing to approximately 120-130° in 
the C3 field-of-view, and the orientation of the ellipse as a function of time is 



shown to rise from 71° to 95° then fall to 64° (see bottom Figure 4.5). 

This event is modelled satisfactorily with a constant deceleration. However, 
due to the impulsive nature of the CME there are only a few frames available for 
analysis, making it difficult to constrain the kinematics. 



4.3.4 Faint CME: 23 April 2001 

This CME was first observed in the south-west at 12:39 UT on 23 April 2001 and 
exhibits some degree of streamer deflection. 

The height-time plot has a similar trend to CDAW (overplotted in top Fig- 
ure 4.6 with a dashed line). The velocity fit was found to be linearly decreasing 
from 581 to 466 km s _1 , giving a deceleration of —4.8 ±4.1 m s™ 2 . The ellipse fit 
spans approximately 55-60° of the field-of-view throughout the event, and the 
orientation of the ellipse as a function of time is shown to decrease from approx- 



imately 50° to almost 0° (see Figure 4.6) 



This CME is fitted well with the constant acceleration model. 
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Figure 4.5: Kinematic and morphological profiles for the ellipse fit to the multi- 
scale edge detection of the 23 April 2000 CME observed by LASCO/C2 and C3. 
The plots from top to bottom are height, velocity, angular width, and ellipse tilt. 
The CDAW heights are over-plotted with a dashed line. The height and velocity 
fits are based upon the constant acceleration model. 
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Figure 4.6: Kinematic and morphological profiles for the ellipse fit to the multi- 
scale edge detection of the 23 April 2001 CME observed by LASCO/C2 and C3. 
The plots from top to bottom are height, velocity, angular width, and ellipse tilt. 
The CDAW heights are over-plotted with a dashed line. The height and velocity 
fits are based upon the constant acceleration model. 
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4.3.5 Fast CME: 21 April 2002 

This CME was first observed in the west from 01:27 UT on 21 April 2002. 

The height-time plot follows a similar trend to that of CDAW (overplotted 



in top Figure 4.7 with a dashed line). The velocity fit was found to be lin- 
early increasing from 2,195 to 2,412 km s -1 , giving a constant acceleration fit of 
32.5 ±26.6 m s" 2 . The ellipse fit spans 53° in C2, and shows an increasing trend 
to 65° in C3. The orientation of the ellipse as a function of time is shown to 
scatter about 115° though it drops to approximately 81° in the final C3 image. 

The kinematics of this event are not modelled satisfactorily by the constant 
acceleration model, since the fits do not lie within all error bars. The argument 
for a non-linear velocity profile, with a possible early decreasing acceleration, is 
justified for this event, although the instrument cadence limits the data set avail- 



able for interpretation. The previous analysis of Gallagher et al. (2003) resulted 



in a velocity of ~ 2,500 km s 1 past ~ 3.4 R which is consistent with our results 



past ~6 R in Figure 4.7 



4.3.6 Flux-Rope/Slow CME: 1 April 2004 

This CME was first observed in the north-east from approximately 23:00 UT on 
1 April 2004, is in the field-of-view for over 9 hours, and exhibits a bright loop 
front, cavity and twisted core. 

The height-time plot follows a similar trend to that of CDAW (overplotted 



in top Figure 4.8 with a dashed line). The velocity fit was found to be linearly 
increasing from 415 to 570 km s _1 , giving an acceleration of 4.4 ± 2.0 m s~ 2 . Note 
also that the kinematics of this event exhibit non-linear structure clearly seen in 
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Figure 4.7: Kinematic and morphological profiles for the ellipse fit to the multi- 
scale edge detection of the 21 April 2002 CME observed by LASCO/C2 and C3. 
The plots from top to bottom are height, velocity, angular width, and ellipse tilt. 
The CDAW heights are over-plotted with a dashed line. The height and velocity 
fits are based upon the constant acceleration model. 
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the velocity and acceleration profiles. The ellipse fit spans approximately 44° in 
C2, stepping down to approximately 38° in C3. The orientation of the ellipse 
as a function of time is shown to jump down from approximately 130° in C2 to 
approximately 70-80° in C3. 

This event shows unexpected structure in the velocity and acceleration profiles 
which indicates a complex eruption not satisfactorily modelled with constant 
acceleration. 



4.3.7 STEREO-B Event: 8 October 2007 



This CME was first observed in the west from approximately 12:00 UT on 8 
October 2007, and is best viewed from the STEREO-B spacecraft. It is noted 
that the kinematics as measured by SOHO and STEREO will be different due 



to projection effects (Howard et al., 2008b Vrsnak et al, 2007). On this date 
STEREO-B was at an angular separation of 16.5° from Earth. 



The height-time plot for this event is shown in Figure |479| The velocity fit was 
found to be linearly increasing from 71 to 330 km s _1 , giving an acceleration of 
5.7 ±0.9 m s~ 2 . The ellipse fit in COR1 spans approximately 23° stepping up to 
a scatter about 40-50° which rises slightly to 50-60° in COR2. The orientation 
of the ellipse as a function of time is shown to increase from 55 - 110° then jumps 
to an approximately steady scatter about 180-190°. 

This CME is fitted well with the constant acceleration model. 
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Figure 4.8: Kinematic and morphological profiles for the ellipse fit to the mul- 
tiscale edge detection of the 1 April 2004 CME observed by LASCO/C2 and C3. 
The plots from top to bottom are height, velocity, angular width, and ellipse tilt. 
The CDAW heights are over-plotted with a dashed line. The height and velocity 
fits are based upon the constant acceleration model. 
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Figure 4.9: Kinematic and morphological profiles for the ellipse fit to the multi- 
scale edge detection of the 8 October 2007 CME observed by SECCHI/COR1 and 
COR2 on STEREO-B. The plots from top to bottom are height, velocity, angular 
width, and ellipse tilt. The height and velocity fits are based upon the constant 
acceleration model. 
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4.3.8 STEREO-A Event: 16 November 2007 

This CME was first observed in the west from approximately 08:26 UT on 16 
November 2007, and is best viewed from the STEREO-A spacecraft. On this 
date STEREO-A was at an angular separation of 20.3° from Earth. 



The height-time plot for this event is shown in Figure 4.10 The velocity fit 
was found to be linearly increasing from 131 to 483 km s -1 , giving an acceleration 
of 13.7 ± 1.7 m s -2 . The ellipse fit in COR1 spans approximately 40-50° stepping 
up slightly to a scatter about 45-55° in COR2. The orientation of the ellipse as 
a function of time is shown to start at 153° and end at 120° with the mid points 
scattered about 170°. 

This CME is fitted well with the constant acceleration model. 



4.4 Discussion & Conclusions 

A variety of theoretical models have been proposed to describe CMEs, especially 
their early propagation phase. Observational studies, such as those outlined 
above, are necessary to determine CME characteristics. We argue that the results 
of previous methods are limited in this regard due mainly to large kinematic errors 
which fail to constrain a model, an artefact of CME detection based upon either 
running- (or fixed-) difference techniques or other operations. Current methods 
fit either a linear model to the height-time curve, implying constant velocity and 
zero acceleration (e.g., CACTus) or a second order polynomial, producing a linear 
velocity and constant acceleration (e.g., CDAW, SEEDS). The implementation 
of a multiscale decomposition provides a time error on the scale of seconds (the 
exposure time of the instrument) and a resulting height error on the order of a 
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Figure 4.10: Kinematic and morphological profiles for the ellipse fit to the mul- 
tiscale edge detection of the 16 November 2007 CME observed by SECCHI/COR1 
and COR2 on STEREO-A. The plots from top to bottom are height, velocity, 
angular width, and ellipse tilt. The height and velocity fits are based upon the 
constant acceleration model. 
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few pixels. The height-time error is used to determine the errors of the velocity 
and acceleration profiles of the CMEs. It was shown that for certain events the 
results of CACTus, CDAW and SEEDS can differ significantly from our methods, 
as illustrated in the Tables of Section 14.31 

Our results clearly confirm that the constant acceleration model may not 
always be appropriate. The 2 January 2000 and 21 April 2002 CMEs are good 
examples of the possible non-linear velocity profile and consequent non-constant 



acceleration profile (see Figure 4.3 and Figure 4.7). Indeed these events are 
shown to have a decreasing acceleration, possibly to zero or below, as the CMEs 



traverse the field-of-view. Simulations of the breakout model outlined in Lynch 



et al. (2004) resulted in constant acceleration fits which do not agree with these 



observations. It may be further noted that the events of 23 April 2001 and 1 April 
2004 show a possible decreasing acceleration phase early on, though within errors 



this cannot be certain (see Figure 4.6 and Figure 4.8). Furthermore, the structure 
seen in some events would indicate that the CME does not progress smoothly. 



The velocities of the 1 April 2004 CME in Figure 4.8 and the 16 November 2007 



CME in Figure |4.10| show non-smooth profiles and may imply a form of bursty 
reconnection or other staggered energy release driving the CME. Other profiles 



such as Figure |4.3| and to a lesser extent Figures |4.6| and |4.7| may show a stepwise 
pattern, indicative of separate regimes of CME progression. None of the current 
CME models indicate a form of non-smooth progression, although the flux-rope 
model does describe an early acceleration regime giving a non-linear velocity to 



the eruption (see Figure 11.5 in Priest & Forbes, 2000). 

It may be concluded that the angular widths of the events are indicative of 
whether the CME expands radially or otherwise in the plane-of-sky. For the 
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CMEs studied above, the observations of 18 April 2000, 23 April 2000, and 21 
April 2002 show an angular width expansion (see Figure 



4.4 



Figure 4.5, and 



Figure 4.7). These events also show high velocities, obtaining top speeds of up 
to 1,000 km s -1 , over 1,100 km s -1 and 2,500 km s _1 respectively, and may 
therefore indicate a link between the CME expansion and speed. Furthermore, it 



is suggested by Krall h St. Cyr (2006) that the flux-rope model can account for 
different observed expansion rates due to the axial versus broadside view of the 
erupting flux system. 

The observed morphology of the ellipse fits may be further interpreted through 



the tilt angles plotted in Section 4.3 In knowing the ellipse tilt and the direction 
of propagation of the CME it is possible to describe the curvature of the front. For 
the events above, the changing tilt and hence curvature is possibly significant for 



the 18 April 2000, 1 April 2004, and 8 October 2007 events (see bottom Figure 4.4 



Figure 4.8, and Figure 4.9). The elliptical flux rope model of Krall et al. (2006) 
was shown to have a changing orientation of the magnetic axis which results in a 
dynamic radius of curvature of the CME, possibly accounting for these observed 
ellipse tilts. 

The work outlined here is an initial indication that the zero and constant 
acceleration models in CME analysis are not an accurate representation of all 
events, and the overestimated angular widths are not indicative of the true CME 
expansion. The ellipse characterisation has provided additional information on 
the system through its changing width and orientation. This work will be further 
explored and developed with STEREO data whereby the combined view-points 
can give additional kinematic constraints and lead to a correction for projection 
effects through 3D reconstructions (discussed in Chapter [5]). 



138 



Chapter 5 

Propagation of an Earth-Directed 
CME in Three-Dimensions 

CMEs are long known to be significant drivers of adverse space weather at Earth, 
but the physics governing their propagation is not fully understood. The launch 
of the STEREO mission in 2006 has provided new insight into their motion in the 
heliosphere, although the mechanisms governing their evolution remain unclear 
due to difficulties in reconstructing their true 3D structure. Here we use a new 
elliptical tie-pointing technique to reconstruct a full CME front in 3D, enabling 
us to quantify an early acceleration profile, non-radial motion, increasing angular 
width and 'pancaking' of the CME front as it propagates from 2-46 R . Be- 
yond 7 R , we show that its motion is determined by aerodynamic drag in the 
solar wind and, using our reconstruction as input for a 3D MHD simulation, we 
determine an accurate arrival time at the LI point near Earth. This chapter is 
founded on research published in Byrne et ah, Nature Communications (2010). 
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5.1 Introduction 



It is predominantly believed that magnetic reconnection is responsible for the 
destabilisation of magnetic flux ropes on the Sun, which then erupt through the 



corona into the solar wind to form CMEs (Moore & Sterling, 2006). There is 
much debate as to the specific processes which trigger the eruption of CMEs, 



and different models exist to explain these (Antiochos et al. 1999; Chen, 1996 



Forbes & Priest, 1995; Kliem & Torok 2006 van der Hoist et al. 2007). In the 



low solar atmosphere, it is postulated that high latitude CMEs undergo deflec- 
tion since they are often observed at different position angles with respect to their 



associated source region locations (Xie et al, 2009). It has been suggested that 
field lines from polar coronal holes may guide high-latitude CMEs towards the 



equator (Kilpua et al., 2009), or that the initial magnetic polarity of a flux rope 



relative to the background magnetic field influences its trajectory (Chane et al. 



2005 Filippov et al, 2001). During this early phase, CMEs are observed to ex- 



pand outwards from their launch site, though plane-of-sky measurements of their 



increasing sizes and angular widths are ambiguous in this regard (Gopalswamy 



et al, 2009a). This expansion has been modelled as a pressure gradient between 



the flux rope and the background solar wind (Odstrcil & Pizzo, 1999 Riley & 



Crooker, 2004). At larger distances in their propagation, CMEs are predicted to 



interact with the solar wind and the interplanetary magnetic field. Studies that 
compare in-situ CME velocity measurements with initial eruption speeds through 
the corona show that slow CMEs must be accelerated toward the speed of the 



solar wind, and fast CMEs decelerated ( Gonzalez-Esparza et al. , 2003 ; Maloney 



et al. , 2009 ). It has been suggested that this is due to the effects of drag acting on 
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the CME in the solar wind (Cargill 2004 Tappin, 2006). However, the quantifi- 
cation of drag, along with that of both CME expansion and non-radial motion, 
is currently lacking, due primarily to the limits of observations from single fixed 
viewpoints with restricted fields-of-view. 

Efforts to reconcile 2D plane-of-sky images with the true 3D morphology of 
CMEs have been underway since they were first observed in the 1970s. The inher- 
ent difficulties in this are predominantly due to the single, fixed-position imagers 
with restricted fields-of-view, as well as the difficulties in observing the optically 
thin coronal plasma of these dynamic events. Before the launch of STEREO, 
there was limited ability to infer the 3D CME morphology from the available ob- 
servations such as SOHO/LASCO. Coronagraphs mainly measure the Thomson 
scattered light of the free electrons in the coronal plasma, providing white-light 
images of CMEs against the plane-of-sky that are not trivial to deconvolve, and 
the projected 2D nature of these images introduces uncertainties in kinemati- 



cal and morphological analyses (Vrsnak et al., 2007). Some efforts were based 



upon a pre-assumed geometry of the CME, such as the cylindrical model (Cre- 



mades & Bothmer, 2004) or the cone model (Xue et al. 2005 Zhao et al. 2002), 



whose shapes were simply oriented to best match the 2D observations. Oth- 
ers used either a comparison of multiple events to infer a statistical relationship 



between plane-of-sky measurements and true CME motion (Howard & Tappin 



2005 Schwenn et al, 2005), or a comparison of observations with in-situ data 



and/or signatures on-disk (Demoulin et al. 2008 Howard et al. 2008b). One 
prominent method was the use of 3D polarisation analysis of LASCO images 



(Moran & Davila, 2004), whereby the line-of-sight averaged distance from the 
plane-of-sky is determined from the brightness ratio of polarised to unpolarised 
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electron scattered emissivity (K-corona). However, this lacks in details such as 
whether the feature is truly unique along the line-of-sight, and, if so, is it towards 
or away from the observer with respect to the plane-of-sky. Polarisation analy- 
sis itself is only acceptable up to heights of ~5 R , since beyond this distance 



the dust-scattered F-corona may no longer be considered unpolarised (Billings 



1966). These issues motivated the launch of the STEREO mission to further our 



understanding of CMEs. 



5.2 The STEREO Era 

The two near-identical spacecraft of the STEREO mission provide simultaneous 
observations of CMEs from independent viewpoints to better observe their true 
morphology. Unfortunately there are limitations on how much 3D information 
can be extracted from the combined plane-of-sky observations, especially when 
the object is optically thin and its boundaries ill-defined. In order to determine 
the morphology of an object in 3D from only two viewpoints, techniques must be 



applied within the context of an epipolar geometry, as illustrated in Figure 5.1 



(Inhester, 2006). This epipolar coordinate system for considering the 3D space 
observed from two independent viewpoints is built up as follows. A line is drawn 
to connect the two observers, called the stereo base line. The two observer loca- 
tions and any third object point or location in the observing space then define 
a plane. Numerous object points will define numerous planes that share an in- 



tersection with the stereo base line. These are the epipolar planes of Figure |5.1 



The plane-of-sky as seen by each observer then intersects the epipolar planes 
such that they appear as epipolar lines across the image, and will converge on a 
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epipolar; 
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Figure 5.1: Schematic of the epipolar geometry used to relate the observations 
from the two STEREO spacecraft, reproduced from Inhester (2006 ). This geometry 
enables us to localise features in 3D space by the triangulating lines-of-sight across 
epipolar planes. 



point along the stereo base line referred to as the epipole of that image. So if a 
line-of-sight from observer 1 is drawn across an epipolar plane, it will appear as a 
single point on image 1, but as a complete line across the corresponding epipolar 
plane in image 2 as seen by observer 2, who is then able to triangulate upon an 
object in 3D space by the intersecting lines-of-sight. This technique is known as 
tie-pointing. 

The technique of tie-pointing lines-of-sight across epipolar planes is best for re- 



solving a single feature, such as a coronal loop on-disk (Aschwanden et al. , 2008). 
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or 

observer 2 



observer 1 



Figure 5.2: Schematic of how tie-pointing a curved surface within an epipolar 
geometry is limited in its ability to resolve the true feature, since lines-of-sight will 
be tangent to different edges of the surface and not necessarily intersect upon it. 



Reproduced from Inhester (2006). 



Under the assumption that the same feature may be tracked in coronagraph im- 



ages many CME studies have also employed tie-pointing techniques (Liewer et al. 



2009 Mierla et al. , 2008 Srivastava et al. , 2009 Temmer et al. , 2009 Wood et al. 



2009). However, when measuring the kinematics of the CME front this technique 



alone doesn't hold true, since it is inevitable that the same part of the curved front 
cannot be confidently resolved from both viewpoints once the CME has traversed 
a certain distance in space, nor similarly once the spacecraft have moved beyond 



a certain angular separation during the mission (Figure 5.2). Furthermore, tri- 
angulating CME observations using only the COR images confines the kinematic 
and morphological analyses to within the 20 R field-of-view. The additional use 
of the heliospheric imagers allows a study of CMEs out to distances of ~ 1 AU, 
however (instrumental effects aside) a 3D analysis can only be carried out if the 
CME propagates along a trajectory between the two spacecraft so that it is ob- 
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served by both HI instruments. Otherwise, assumptions of its trajectory have 



to be inferred from either its association with a source region on-disk (Howard 



& Tappin, 2008) or its trajectory through the COR data (Maloney et al., 2009), 



or derived by assuming a constant velocity through the HI fields-of-view (Davis 



et al. , 2009). Triangulation of CME features using time-stacked intensity slices 



at fixed latitude, named 'J-maps' due to the characteristic propagation signature 



of a CME, has also been developed (Davis et al, 2010 Liu et al, 2010). This 
technique is hindered by the same limitation of standard tie-pointing techniques; 
namely that the curvature of the feature is not considered, and the intersection 
of sight-lines may not occur upon the surface of the observed feature. 

An alternative to tie-pointing is a method called forward modelling which pre- 
sumes a given shape of the CME and seeks to match it with observational data. 



Thernisien et al. (2006) employ a graduated cylindrical shell which is warped to 



form a flux rope model overlaid on CME images (Figure 5.3). The parameters 
governing the model's shape and orientation may be changed by the user to fit the 
model to STEREO- Ahead and Behind data simultaneously and obtain a 3D flux 
rope characterisation of the CME as it propagates, though this may not always 



be appropriate (Jacobs et al, 2009). Boursier et al. (2009a) outline a similar 
forward model which assumes one of three pre-assigned shapes: a hemispherical 
cap, a flux rope, or a cloud-like model. However, in each of these methods the 
predetermined shape of the CME model has a spherical cross-section and must 
adhere to some quasi-similarity (self-invariance) over the sequence of images. So 
while forward modelling better accounts for the curved nature of the CME being 
observed, the inherent restrictions of the imposed model still limit the analysis 
of the true 3D structure and dynamics of the CME as it propagates. 
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5.3 Elliptical Tie-Pointing 



In the epipolar geometry outlined above, 3D information may be gleaned from 
two independent viewpoints of a feature using tie-pointing techniques to trian- 
gulate lines-of- sight in space. However, when the object is known to be a curved 
surface, sight-lines will be tangent to it and not necessarily intersect upon it (Fig- 
ure 



5.2). Consequently CMEs cannot be reconstructed by tie-pointing alone, but 



rather their localisation may be constrained by intersecting sight-lines tangent to 



the leading edges of a CME (de Koning et al, 2009 Pizzo & Biesecker, 2004). 
Following the multiscale edge detection and ellipse characterisation outlined in 
Chapter |3j it is possible to extract the intersection of a given epipolar plane 
through the ellipse fits of both the STEREO-Ahead and Behind images. This 
defines a quadrilateral in 3D space which localises the ellipse characterisation of 
the CME front in that plane. 

Inscribing an ellipse within the quadrilateral such that it is tangent to all four 
sides (detailed below) provides a slice through the CME that matches the obser- 
vations from each spacecraft (Figure 5.4k .). A full reconstruction is achieved by 
stacking ellipses from numerous epipolar slices (Figure 5.4b ). Since the positions 
and curvatures of these inscribed ellipses are constrained by the characterised 
curvature of the CME front in the stereoscopic image pair, the modelled CME 
front is considered an optimum reconstruction of the true CME front. This is re- 
peated for every frame of the eruption to build the reconstruction as a function of 



time and view the changes to the CME front as it propagates in 3D (Figure 5.4:). 
Following Horwitz ( 2002| 2005), we inscribe an ellipse within a quadrilateral 



using the following steps (see Figure 5.5): 
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Figure 5.4: The elliptical tie-pointing technique developed to reconstruct the 
3D CME front, shown here for the 12 December 2008 event. One of any number 
of epipolar planes will intersect the ellipse characterisation of the CME at two 
points in each image from STEREO-A and B. (a) illustrates how the resulting 
four sight-lines intersect in 3D space to define a quadrilateral that constrains the 
CME front in that plane. Inscribing an ellipse within the quadrilateral such that 
it is tangent to each sight-line provides a slice through the CME that matches 
the observations from each spacecraft, (b) illustrates how a full reconstruction 
is achieved by stacking multiple ellipses from the epipolar slices to create a model 
CME front that is an optimum reconstruction of the true CME front, (c) illustrates 
how this is repeated for every frame of the eruption to build the reconstruction as 
a function of time and view the changes to the CME front as it propagates in 
3D. While the ellipse characterisation applies to both the leading edges and, when 
observable, the flanks of the CME, only the outermost part of the reconstructed 
front is shown here for clarity. 



148 



5.3 Elliptical Tie-Pointing 




Figure 5.5: An ellipse inscribed within a convex quadrilateral. An isometry of 
the plane is applied such that the quadrilateral has vertices (0,0), (A,B), (0,C), 
(s,t). The ellipse has centre (h,k), semimajor axis a, semiminor axis b, tilt angle 
5, and is tangent to each side of the quadrilateral. 

1. Apply an isometry to the plane such that the quadrilateral has vertices 
(0,0), (A, B), (0,C), (s,t), where A > 0, C > 0, s > and t > B. (Note, 
in the case of an affine transformation we set A — 1, B — and C — 1, 
with s and t variable.) 

2. Set the ellipse centre point (h, k) by fixing h somewhere along the open line 
segment connecting the midpoints of the diagonals of the quadrilateral and 
hence determine k from the equation of a line, for example: 
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3. To solve for the ellipse tangent to the four sides of the quadrilateral, we can 
solve for the ellipse tangent to the three sides of a triangle whose vertices 
are the complex points 



At — Bs 

Zl = o, z 2 = A + Bi, z 3 = —i (5.2) 

s — J\ 



and the two ellipse foci are then the zeroes of the equation 

p h ( z ) = ( s -A)z 2 -2(s-A)(h-ik)z-(B-iA)(s-2h)C (5.3) 

whose discriminant can be denoted by r(h) = ri(h) + ir 2 (h) where 

n = 4((s-Af-(t-B-Cf) (*^) 2 

+ A(s-A) (A(s-A) + B{B-t) + C(C-t)) 

+ 4 (s - A) (At + Cs + Bs- 2AB) 
+ 2A{s-Af(B-C) (5.5) 

Thus we need to determine the quartic polynomial u{h) = \r(h)\ 2 = ri(h) 2 + 
r2(h) 2 and we can then solve for the ellipse semimajor axis, a, and semiminor 
axis, b, from the equations 



a2 - 62 = '/(i6(^) MW M) 



+ (s- A) 2 (A 2 - (C - Bf) (5.4) 

h- A" 



r 2 =8(f- B -C){s- A) 



150 



5.3 Elliptical Tie-Pointing 



a 2 b 2 = - ( < =^ ) (2 (Bs -A(t-C))h- ACs) (2h - A) (2h - sl5.7) 



4 \( s - AY 

by parameterising R = a 2 — b 2 and W = a 2 b 2 to obtain 



l - (VR 2 + AW + i?) , b = (VR 2 + AW - R^j (5 



4. Knowing the axes we can generate the ellipse and float its tilt angle S until 
it sits tangent to each side of the quadrilateral, using the inclined ellipse 



equation (3.15) introduced in Section 3.4 



5.3.1 SOHO as a Third Perspective 

The elliptical tie-pointing technique was used to reconstruct the front of a CME 
observed by STEREO on 26 April 2008 in order to test its efficacy by comparing 
it with observations from SOHO - a third perspective on the event. The CME 
appears as a halo from STEREO-B, so a running-difference technique is used to 
highlight the faint CME front in the images. The front is defined in the images by 
a point-and-click methodology and characterised with an ellipse fit (outlined in 



Section 3.4). From STEREO-A the event appears off the east limb and shows a 
strong streamer deflection to the south-east (in fact the CME would probably be 
considered only as the northern portion of the erupting material in STEREO-A 
if it were not shown by STEREO-B to expand further south). 

With the ellipse characterisations determined for the CME front in the si- 
multaneous images from COR1 and COR2 onboard STEREO-A and B, the el- 
liptical tie-pointing technique is performed and the CME front reconstructed in 
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Figure 5.6: The back-projection of a STEREO 3D CME front reconstruction 
onto the SOHO/LASCO plane-of-sky, from the observations by STEREO-A (red) 
and STEREO-B (green) at 16:22 UT on 26 April 2008. 



3D. This reconstruction is then back-projected onto the LASCO plane-of-sky in 
order to compare it with observations of the CME from SOHO's vantage point 



at LI (Figure 5.6). This back-projection is performed by standard geometry of 
lines-of-sight from the observer position 0(xq, yo, zo) through the 3D reconstruc- 
tion points P(xi, yi, Zi) and determining where they intersect the plane-of-sky 
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Q(xj = 0, i/j, Zj) as follows: 



tana = = — = — =^> yj = Xi + y { (5.9) 



z i ~ z z j ~ Z _ z j ~ z i , { z i ~~ z 



tan/3 = = = => Zj = x t + z> (5.10) 

i_t. 2 «Xy Q J. ^ iAj Q iX> j J.. ^ \^ ■X' 2, **•' Q 



Due to the different instrument cadences of the SECCHI and LASCO coron- 



agraphs, frames which lie closest in time were chosen for comparison. Figure 5.7 



shows the COR2 frames from STEREO-A and B at 16:22 UT on 26 April 2008 
with the ellipse characterisations of the CME front (left and right panels), and the 
back-projected 3D front reconstruction as compared with the LASCO/C2 frame 
from SOHO at 16:30 UT (middle panel). The reconstruction from the STEREO 
observations adequately fits with the SOHO observations given the time offset, 
and so gives credence to the elliptical tie-pointing technique. 



5.4 Earth-Directed CME 

On 12 December 2008 an erupting prominence was observed by STEREO while 



the spacecraft were in near quadrature at 86.7° separation (Figure 5.8 1). The 
eruption is visible at 50-55° north from 03:00 UT in SECCHI/EUVI images, ob- 
tained in the 304 A passband, in the northeast from the perspective of STEREO-A 
and off the northwest limb from STEREO-B. The prominence is considered to be 
the inner material of the CME which was first observed in COR1-B at 05:35 UT 



(Figure 5.8 d). For our analysis, we use the two coronagraphs (COR1/2) and the 



inner Heliospheric Imagers (HI1) (Figure 5.8:). In each image the front of the 
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5.4 Earth-Directed CME 




Figure 5.8: Composite of STEREO- Ahead and Behind images from EUVI, C0R1, 



COR2, and HI1 ( |Byrne et al||2010[ ). (a) indicates the STEREO spacecraft loca- 
tions, separated by an angle of 86.7° at the time of the event, (b) shows the 
prominence eruption observed in EUVI-B off the NW limb from approximately 
03:00 UT which is considered to be the inner material of the CME. The multi- 
scale edge detection and corresponding ellipse characterisation are overplotted in 
COR1. (c) shows that the CME is Earth-directed, being observed off the east limb 
in STEREO- A and the west limb in STEREO-B. 
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Figure 5.9: STEREO-Ahead and Behind COR1 images at 07:35 UT on the 12 
December 2008. Overplotted are: the multiscale edge detections of the CME front 
(red); the ellipse characterisations (blue); and the resulting 3D reconstructions 
back-projected onto the plane-of-sky (white). 
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Figure 5.10: STEREO-Ahead and Behind COR2 images at 14:52 UT on the 12 
December 2008. Overplotted are: the multiscale edge detections of the CME front 
(red); the ellipse characterisations (blue); and the resulting 3D reconstructions 
back-projected onto the plane-of-sky (white). 
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Figure 5.11: STEREO-Ahead and Behind HI1 images at 01:29 UT on the 13 
December 2008. Overplotted are: the running difference edge detections of the 
CME front (red); the ellipse characterisations (blue); and the resulting 3D recon- 
structions back-projected onto the plane-of-sky (white). 
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5.5 Results 



CME is fitted with an ellipse that characterises its propagation across the plane- 



of-sky (Byrne et al. 2009). This ellipse fitting is sensitive predominantly to the 
leading edges of the CME but equal weight is given to the CME flank edges as 
they enter the field-of-view of each instrument. The 3D reconstruction is then 
performed using a method of curvature-constrained tie-pointing within epipolar 



planes containing the two STEREO spacecraft (detailed in Section 5.3). An ex- 
ample of the ellipse characterisation to the CME front and the corresponding 



back-projected 3D reconstruction is shown in Figure 5.9 for COR1, Figure 5.10 



for COR2, and Figure 5.11 for HI1. Corresponding frames from a 3D visualisation 



of the event are shown in Figure |5.12| for COR1, Figure [5. 13| for COR2, and Fig- 
ure 



5.14 for HI1, showing the relative locations of the Sun, Earth and STEREO 



spacecrafts in the inner heliosphere. Figure |5.15| illustrates the graphical user 
interface developed for the 3D visualisation of the CME front reconstruction. 



5.5 Results 

The resulting kinematics and morphology of the CME are measured along an an- 
gular span through the reconstructed CME front in the out-of-ecliptic plane along 



the Sun-Earth line (Figure 5.16). These were taken by first closing tangents to the 
CME front ('Northern/Southern Flanks'), and then measuring the height along 
an angle midway between these ('Midpoint of Front'), and then similarly along 
the two angles midway between the midpoint and the flanks ('Midtop/Midbottom 
of Front'). Although these measurements are fixed along the Sun-Earth line, in- 
vestigating how the CME height profile would change if taken along a trajectory 
slightly off the Sun-Earth line shows no significant deviation within the associ- 
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Figure 5.16: The 3D CME front reconstruction from the COR2 Ahead and Be- 
hind frames at 14:52 UT on 12 December 2008. The lines drawn from Sun-centre 
indicate the 'Midpoint of Front' (solid blue), the 'Northern / Southern Flanks' (solid 
red /brown), and the 'Midtop /Midbottom of Front' at angles in between (dashed 
red /brown). By taking these measurements across all frames we may determine 



the kinematics and morphology of the CME as plotted in Figures 5.18 and 5.19 



ated errors and thus has negligible effect on the kinematics. The same is true if 
the overall maximum height (of varying location) on each individual CME front 
is instead taken and the kinematics reanalysed. 
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Figure 5.17: The error trapezoid on the tie-pointing of two lines-of-sight in 3D 
space, reproduced from Inhester (2006). The error w in localising a point on each 
plane-of-sky results in a trapezoid with a diagonal measuring w/ sin(a/2) as shown, 
where a is the spacecraft separation angle. 



5.5.1 3D Error Propagation 

When considering the errors that propagate from the 2D plane-of-sky of each 
image onto the 3D quadrilateral localising the CME, we may assume that the 
lines-of-sight within the error range are essentially parallel. This means the error 
interval on the coordinate being tie-pointed in 3D is given by a trapezoid sur- 



rounding the intersection of the lines-of-sight, illustrated in Figure 5.17 This is 
done for each corner of the quadrilateral within a given epipolar plane. For a 
spacecraft separation of angle a and errorbar of magnitude w on the 2D image 
we can define the error trapezoid as having diagonals of length wj cos(a/2) and 
w/ sin(a/2). 
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So in the case of C0R1/2, the optimum filter size in the multiscale decompo- 
sition was 2 3 pixels wide, giving an error of ± 8 pixels, so w = 16. Over the course 
of the 12 December 2008 CME the average STEREO spacecraft separation was 
86.75°, so we calculate the error trapezoid as having diagonals of size: 



w 



w 



cos(f)' sin(f) 



[11.0, 11.6] 



(5.11) 



This provides a 3a height error of 11.6 pixels, so the corresponding la height 
error is given by 68(11.6)/99.7 = 7.9 pixels. The time error for the multiscale 
edge detections is given by the exposure time of the individual frames: 1.69984 
seconds for COR1 and 2.00090 seconds for COR2. 



In the case of HI1 a la plane-of-sky error of 3 pixels was determined (Maloney 



et al, 2009), so w = 6 and the error trapezoid is deduced to be [4.1, 4.4]. So the 
height error for HI1 is taken as 4.4 pixels and the time error is given by the thirty 



summed 60 second images to result in 1800 seconds (Eyles et a/. 2009) 



These errors are transformed first into arcseconds by multiplying by the plate 
scale of the instruments (7.5043001 arcsec/pixel for COR1, 14.7 arcsec/pixel for 
COR2, and 71.927554 arcsec/pixel for HI1), and then into metres, knowing the 
respective size of the Sun in arcseconds on the plane-of-sky observed by each 
instrument (given 1 R = 6.95508 x 10 8 m). The resultant errors are then prop- 



agated with the 3-point Lagrangian interpolation (detailed in Section 4.2.1) from 



the height-time curves into the velocity and acceleration profiles of Figure 5.19 



Due to the potentially large deviation of the end points from the general trend in 
3-point Lagrangian interpolation, the endpoints of the EUVI prominence data, 
the COR1/2 coronagraph data, and the HI data are each removed as outliers 
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from the velocity and acceleration plots. 

5.5.2 Prominence Sz CME Acceleration 

In determining the CME kinematics for the 'Midpoint of Front' we find a steep 
increase in the velocity corresponding to an early impulsive acceleration phase 



of 94 ±58 m s (bottom of Figure 5.18). The prominence rises with a velocity 



of ~ 50 km s" 1 before the system fully erupts and the prominence undergoes 
an acceleration of 40 ±5 m s~ 2 behind the CME front. This is indicative of the 
onset of explosive reconnection or other loss-of-equilibrium in the system whereby 
the internal magnetic pressure increases sufficiently and/or the external magnetic 
pressure decreases sufficiently to allow the eruption to proceed in the height range 



1.5-3 R (top panel of Figure 5.19). The velocity profile is synonymous with 



those produced by the 2D flux rope model as in Figure 11.5 of Priest fc Forbes 



(2000). The acceleration then reduces to scatter about zero as the explosive 



nature of the eruption due to the Lorentz force diminishes and the drag force due 



to the ambient solar wind pressure begins to dominate (see Equation 1.26). 



5.5.3 Non-radial Prominence Sz CME Motion. 



It is immediately evident from the reconstruction (illustrated in Figure 5.4:) that 
the CME propagates non-radially away from the Sun. The CME flanks change 
from an initial latitude span of 16-46° to finally span approximately ±30° of the 



ecliptic (middle panel of Figure 5.19 ). The mean declination, 9, of the CME is well 
fitted by a power-law of the form 9{r) = 6*or~ a92 (2 R Q < r < 46 R ) as a result 
of this non-radial propagation. Tie-pointing the prominence apex and fitting a 
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Figure 5.18: The kinematics of the prominence and 3D reconstructed CME front 
of 12 December 2008. The prominence is observed as the inner material of the CME, 
with both undergoing acceleration from ~ 06:00-07:00 UT, peaking at ~ 40 m s~ 2 
and ~ 94 m s -2 respectively, before reducing to scatter about zero. Measurement 
uncertainties are indicated by one standard deviation error bars. 
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Figure 5.19: Kinematic and morphological properties of the 3D reconstruction of 
the 12 December 2008 CME front. The top panel shows the velocity of the middle 
of the CME front with corresponding drag model and, inset, the early acceleration 
peak. Measurement uncertainties are indicated by one standard deviation error- 
bars. The middle panel shows the declinations from the ecliptic (0°) of an angular 
spread across the front between the CME flanks with a power-law fit indicative of 
non-radial propagation. The bottom panel shows the angular width of the CME 
with a power-law expansion. 
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power-law to its declination angle results in prom (r) = $P rom r ~°- 82 (\ R < r < 
3 R©), implying a source latitude of 6^ om {l R ) « 54° N in agreement with 
EUVI observations. Previous statistics on CME position angles have shown that, 
during solar minimum, they tend to be offset closer to the equator as compared 



to those of the associated prominence eruption (Gopalswamy et al, 2003). The 
non-radial motion we quantify here may be evidence of the drawn-out magnetic 
dipole field of the Sun, an effect predicted at solar minimum due to the influence 



of the solar wind pressure (e.g., Figure 8 in Pneuman & Kopp ( 1971 ) and Figure 2 



in 



Banaszkiewicz et au (1998)). Other possible influences include changes to the 



internal current of the magnetic flux rope (Filippov et al, 2001), or the orientation 



of the magnetic flux rope with respect to the background field (Chane et al. 



2005), whereby magnetic pressure can act asymmetrically to deflect the flux rope 



pole-ward or equator-ward depending on the field configurations. 



5.5.4 CME Angular Width Expansion 

Over the height range 2-46 R the CME angular width (A9 = 9 max — 8 min ) in- 
creases from ~ 30° to ~ 60° with a power- law of the form A6(r) = A# r a22 (2 R < 



r < 46 R ) (bottom panel of Figure 5.19). This angular expansion is evidence 



for an initial overpressure of the CME relative to the surrounding corona (coinci- 



dent with its early acceleration inset in top panel of Figure 5.19). The expansion 
then tends to a constant during the later drag phase of CME propagation, as it 
expands to maintain pressure balance with heliocentric distance. It is theorised 
that the expansion may be attributed to two types of kinematic evolution, namely 
spherical expansion due to simple convection with the ambient solar wind in a 
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diverging geometry, and expansion due to a pressure gradient between the flux 



rope and solar wind (Tappin, 2006). It is also noted that the southern portions of 
the CME manifest the bulk of this expansion below the ecliptic (best observed by 
comparing the relatively constant 'Midtop of Front' measurements with the more 
consistently decreasing 'Midbottom of Front' measurements in middle panel of 



Figure 5.19). Inspection of a Wang-Sheeley-Arge (WSA) solar wind model run 



( Arge & Pizzo , 2000 ) reveals higher speed solar wind flows (~ 650 km s _1 ) emanat- 
ing from open-field regions at high/low latitudes (approximately 30° north/south 
of the solar equator). Once the initial prominence/CME eruption occurs and is 
deflected into a non-radial trajectory, it undergoes asymmetric expansion in the 
solar wind. It is prevented from expanding upwards into the open-field high-speed 
stream at higher latitudes, and the high internal pressure of the CME relative to 
the slower solar wind near the ecliptic accounts for its expansion predominantly 
to the south. In addition, the northern portions of the CME attain greater dis- 
tances from the Sun than the southern portions as a result of this propagation in 
varying solar wind speeds, an effect predicted to occur in previous hydrodynamic 



models (Odstrcil & Pizzo, 1999). 



5.5.5 CME Drag in the Inner Heliosphere 



Investigating the midpoint kinematics of the CME front, we find the velocity 



profile increases from approximately 100-300 km s over the first 2-5 R(7 



before rising more gradually to a scatter between 400 - 550 km s 1 as it propagates 



outward (top panel of Figure 5.19). The acceleration peaks at approximately 
100 m s -2 at a height of ~ 3 R , then decreases to scatter about zero. This early 
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phase is generally attributed to the Lorentz force whereby the dominant outward 
magnetic pressure overcomes the internal and/or external magnetic field tension, 
while the subsequent increase in velocity, at heights above ~ 7 R , is predicted by 



theory to result from the effects of drag (Tappin, 2006). At large distances from 



the Sun, during this postulated drag-dominated epoch of CME propagation, the 
equation of motion can be cast in the form: 



■M cme ^Pswi^cme ^sui)|^cme ^ sw\A. cme C £) (5.12) 



This describes a CME of velocity v cme , mass M cme , and cross-sectional area A cme 
propagating through a solar wind flow of velocity v sw and density p sw . The drag 
coefficient, Cp, is found to be of the order of unity for typical CME geometries 



(Cargill, 2004), while the density and area are expected to vary as power-law 



functions of distance R. Thus, we parameterise the density and geometric varia- 



tion of the CME and solar wind using a power-law (Vrsnak & Gopalswamy, 2002) 
to obtain: 

= -aR~ p -^— {v sw - v cme y (5.13) 

where 7 describes the drag regime, which can be either viscous (7 = 1) or aero- 
dynamic (7 = 2), and a and /3 are constants primarily related to the cross- 
sectional area of the CME and the density ratio of the solar wind flow to the 
CME (psw/Pcme)- We determine a theoretical estimate of the CME velocity as a 
function of distance by numerically integrating this equation using a 4th order 
Runge-Kutta scheme and fitting the result to the observed velocities from ~7- 
46 R Q . The initial CME height, CME velocity, asymptotic solar wind speed, 
and a, (3, and 7 are obtained from a bootstrapping procedure which provides 
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a final best-fit to the observations and confidence intervals for the parameters. 
Best-fit values for a and (3 were found to be (A.hht 2 ^) x 10 ~ 5 and -^-^^-om which 
agree with values found in previous modelling work ( |Vrsnak 2001). The best-fit 
value for the exponent of the velocity difference between the CME and the solar 
wind, 7, was found to be 2.271q3q, which is clear evidence that aerodynamic drag 
(7 = 2), and not viscous drag (7 = 1) acts during the propagation of the CME 
in interplanetary space. 



5.5.6 CME Arrival Time 

The drag model provides an asymptotic CME velocity of 555iy> 4 km s _1 when 
extrapolated to 1 AU, which predicts the CME to arrive one day before the Ad- 
vanced Composition Explorer (ACE) or WIND spacecrafts detect it at the LI 
point. We investigate this discrepancy by using our 3D reconstruction to simu- 
late the continued propagation of the CME from the Alfven radius (~21.5 R Q ) 



to Earth using the ENLIL with Cone Model flXie et al\ |2004[ ) at NASA's Com- 
munity Coordinated Modeling Center. ENLIL is a time-dependent 3D MHD 
code that models CME propagation through interplanetary space. An ideal fluid 
approximation is used to describe the solar wind plasma, under time-dependent 
MHD processes (neglecting microscopic processes). The plasma is treated as a 
fully ionised hydrogen gas with equal electron and proton densities (n = n e = n p ) 
and temperatures (T = T e = T p ), and the basic equations of MHD theory applied 



(such as outlined in Section 1.3.1). 



We use the height, velocity, and width from our 3D reconstruction as ini- 
tial conditions for the simulation, and find that the CME is actually slowed to 
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Figure 5.20: The 3D CME front parameters are used as initial conditions for an 



ENLIL with Cone Model MHD simulation (Xie et al. , 2004) and the output density 
(top) and velocity (bottom) profiles of the inner heliosphere are illustrated here for 
the time-stamp of 06:00 UT on 14 December 2008. Beyond distances of ~ 50 R© 
the CME is slowed by its interaction with the upstream, slow-speed, solar wind 
flow along its trajectory towards Earth, and this accounts for its arrival time as 
detected in-situ by the ACE and WIND spacecraft at the LI point near Earth. 
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Figure 5.21: The in-situ solar wind plasma and magnetic field data observed by 
the WIND spacecraft. From top to bottom the panels show proton density, bulk 
flow speed, proton temperature, and magnetic field strength and components. The 
red dashed lines indicate the arrival time of the density enhancement predicted from 
our ENLIL with Cone Model run providing 08:09 UT on 16 December 2008, with 
a potential offset error between our reconstruction and the derived model height 
profiles accounting for an arrival time up to 13:20 UT. We observe a magnetic cloud 
signature behind the front, as highlighted with blue dash-dotted lines. 



175 



5. PROPAGATION OF AN EARTH-DIRECTED CME IN 
THREE-DIMENSIONS 



~342 km s" 1 at 1 AU. This is as a result of its interaction with an upstream, 
slow-speed, solar wind flow at distances beyond 50 R , as seen by inspection of 
the solar wind profile along the trajectory of the CME in the ENLIL simulation 



(Figure 5.20). This CME velocity is consistent with in-situ measurements of solar 
wind speed (~330 km s _1 ) from the ACE and WIND spacecraft at LI. Track- 
ing the peak density of the CME front from the simulation gives an arrival time 
at LI of ~ 08:09 UT on 16 December 2008. Accounting for the offset in CME 
front heights between our 3D reconstruction and ENLIL simulation at distances 
of 21.5 R Q < r < 46 R gives an arrival time in the range 08:09-13:20 UT on 16 
December 2008. This prediction interval agrees well with the earliest derived ar- 
rival times of the CME front plasma pileup ahead of the magnetic cloud flux rope 



from the in-situ data of both ACE and WIND (Figure 5.21 ) before its subsequent 



impact at Earth (Davis et al, 2009 Liu et all 2010). 



5.5.7 CME 'Pancaking' 

From the ENLIL simulation it is apparent that the CME undergoes an effect 
known as 'pancaking' whereby the middle portion of the CME may be slowed 
while the flanks of the CME maintain or increase speed such that the front distorts 
to become concave outwards in shape. This is illustrated in the density plot for 



the cross-section of the CME along the Sun-Earth line in top of Figure 5.20 and 
the effect increases with distance from the Sun. We investigate the curvature 
of the 3D CME front reconstruction along the Sun-Earth line in the distance 2 - 
46 R by fitting the front in this plane with an ellipse and inspecting the changing 
morphology with distance. A plot of this characterisation against height is shown 
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in Figure 5.22 where it can be seen that the curvature of the front decreases as the 
CME propagates, with the front initially optimised by high-curvature horizontal 
ellipses, then becoming better optimised by more spherical ellipses, before finally 
being optimised by smoother vertical ellipses. The observations through the 
latter half of the HI1 into the HI2 fields-of-view also show this pancaking effect 
on the CME (though geometrical and instrumental effects must be considered 
when interpreting these images). 



5.6 Discussion & Conclusions 

Since its launch, the dynamic twin-viewpoints of STEREO have enabled studies 
of the true propagation of CMEs in 3D space. Our new elliptical tie-pointing 
technique uses the curvature of the CME front as a necessary third constraint 
on the two viewpoints to build an optimum 3D reconstruction of the front. Here 
the technique is applied to an Earth-directed CME, to reveal numerous forces at 
play throughout its propagation. 

The early acceleration phase results from the rapid release of energy when the 
CME dynamics are dominated by outward magnetic and gas pressure forces. Dif- 
ferent models can reproduce the early acceleration profiles of CME observations 



though it is difficult to distinguish between them with absolute certainty (Lin 



et al, 2010 Schrijver et al, 2008). For this event the acceleration phase coin- 
cides with a strong angular expansion of the CME in the low corona, which tends 
toward a constant in the later observed propagation in the solar wind. While, 



statistically, expansion of CMEs is a common occurrence (Bothmer & Schwenn 



1994), it is difficult to accurately determine the magnitude and rate of expansion 
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Figure 5.22: Top: Ellipses characterising the 3D CME front reconstruction in 
the out-of-ecliptic plane along the Sun-Earth line. Bottom: The ellipse tilt angle 
is indicative of the initial effects of 'pancaking', as the curvature of the CME front 
decreases with increasing height due to its changing morphology in the solar wind. 
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across the 2D plane-of-sky images for individual events. Some studies of these 
single-viewpoint images of CMEs use characterisations such as the cone model 



(Xie et al., 2004 Xue et al, 2005; Zhao et al., 2002) but assume the angular 
width to be constant (rigid cone) which is not always true early in the events 



(Byrne et al, 2009 Gopalswamy et al, 2009a). Our 3D front reconstruction 
overcomes the difficulties in distinguishing expansion from image projection ef- 
fects, and we show that early in this event there is a non-constant, power-law, 
angular expansion of the CME. Theoretical models of CME expansion generally 
reproduce constant radial expansion, based on the suspected magnetic and gas 
pressure gradients between the erupting flux rope and the ambient corona and so- 



lar wind ( |Berdichevsky et al] [20031 |Cargill et al] |2000| |Odstrcil fc Pizzo] p99| ) . 
To account for the angular expansion of the CME, a combination of internal 
overpressure relative to external gas and magnetic pressure drop-offs, along with 
convective evolution of the CME in the diverging solar wind geometry, must be 



considered (Riley & Crooker, 2004). 

During this early phase evolution the CME is deflected from a high-latitude 
source region into a non-radial trajectory as indicated by the changing inclina- 



tion angle (middle panel of Figure 5.19). While projection effects again hinder 
interpretations of CME position angles in single images, statistical studies show 
that, relative to their source region locations, CMEs have a tendency to deflect 



toward lower latitudes during solar minimum ( Gopalswamy et al. , 2003 Yashiro 



et al., 2004). It has been suggested that this results from the guiding of CMEs 



towards the equator by either the magnetic fields emanating from polar coronal 



holes (Kilpua et al. , 2009 Xie et al. , 2009), or the flow pattern of the background 



coronal magnetic field and solar wind/streamer influences (Cremades & Bothmer 
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20041 iMacQueen et all 119861 |Xie et al.\ 120091). Other models show that the in- 



ternal configuration of the erupting flux rope can have an important effect on its 
propagation through the corona. The orientation of the flux rope, either normal 
or inverse polarity, will determine where magnetic reconnection is more likely to 
occur, and therefore change the magnetic configuration of the system to guide 



the CME either equator- or pole- ward (Chane et al., 2005). Alternatively, mod- 
elling the filament as a toroidal flux rope located above a mid-latitude polarity 
inversion line results in non-radial motion and acceleration of the filament, due to 



the guiding action of the coronal magnetic field on the current motion (Filippov 



et al. , 2001). Both of these models have a dependence on the chosen background 



magnetic field configuration, and so the suspected drawn-out magnetic dipole 



field of the Sun by the solar wind ( Banaszkiewicz et al. , 1998 Pneuman & Kopp 



1971) may be the dominant factor in deflecting the prominence/CME eruption 



into this observed non-radial trajectory. 

At larger distances from the Sun (> 7 R ) the effects of drag become impor- 
tant as the CME velocity approaches that of the solar wind. The interaction 
between the moving magnetic flux rope and the ambient solar wind has been 
suggested to play a key role in CME propagation at large distances where the 



Lorentz driving force and the effects of gravity become negligible (Chen, 1996). 
Comparisons of initial CME speeds and in-situ detections of arrival times have 



shown that velocities converge on the solar wind speed ( Gonzalez-Esparza et al. 



2003 Maloney et al., 2009). For this event we find that the drag force is in- 



deed sufficient to accelerate the CME to the solar wind speed, and quantify that 
the kinematics are consistent with the quadratic regime of aerodynamic drag 
(turbulent, as opposed to viscous, effects dominate). The importance of drag 
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becomes further apparent through the CME interaction with a slow-speed solar 
wind stream ahead of it, slowing it to a speed that accounts for the observed 
arrival time at LI near Earth. This agrees with the conjecture that Sun-Earth 
transit time is more closely related to the solar wind speed than the initial CME 



speed (Vrsnak et ai, 2009). Other kinematic studies of this CME through the HI 
fields-of-view quote velocities of 411 ±23 km s _1 (Ahead) and 417 ±15 km s _1 
(Behind) when assumed to have zero acceleration during this late phase of prop- 



agation (Davis et al., 2009), or an average of 363 ±43 km s when triangulated 



in time-elongation J-maps (Liu et al, 2010). These speeds through the HI fields- 
of-view, lower than those quantified through the COR1/2 fields-of-view, agree 
somewhat with the deceleration of the CME to match the slow-speed solar wind 
ahead of it in our MHD simulation. Ultimately we are able to predict a more 
accurate arrival time of the CME front at LI. 

A cohesive physical picture for how the CME erupts, propagates, and expands 
in the solar atmosphere remains to be fully developed and understood from a theo- 
retical perspective. Realistic MHD models of the Sun's global magnetic field and 
solar wind are required to explain all processes at play, along with a need for 
adequate models of the complex flux rope geometries within CMEs. Addition- 
ally, spectral observations of CME onset signatures using such instrumentation 



as the Coronal Diagnostic Spectrometer (CDS; Harrison et al, 1995) on SOHO 



and EUV Imaging Spectrometer (EIS; Culhane et al., 2007) on Hinode, along 



with future space exploration missions such as Solar Orbiter (ESA; |Hassler et al. 



2009) and Solar Probe± (NASA; McComas et al. , 2008), will be required to give 



us a better understanding of the fundamental plasma processes responsible for 
driving CMEs and determining their adverse effects at Earth. 
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Chapter 6 



Conclusions & Future Work 



This thesis has sought to increase our understanding of solar activity and its 
effects on Earth through a study of the phenomenon known as CMEs. This was 
undertaken by studying CME propagation with new data and techniques. This 
chapter presents the main results and conclusions, and outlines possible future 
directions for this work. 



6.1 Principal Results 

The primary objective of this study was to further our understanding of the kine- 
matic and morphological evolution of CMEs as they propagate from the Sun into 
the heliosphere. This was done by applying new methods of multiscale image 
processing to CME observations in order to identify and track the CME front 
and characterise its propagation through coronagraph data. Following this, the 
development of a new elliptical tie-pointing technique applied to STEREO ob- 
servations allowed a reconstruction of a CME front in 3D, overcoming projection 
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effects and revealing its true 3D motion. The principal results arising from these 
studies may be summarised as follows: 

1. The multiscale nature of CMEs was revealed through the application of a 



high- and low-pass image filtering technique (Young & Gallagher, 2008). 
The multiscale filtering was shown to effectively suppress noise and small- 
scale features in order to reveal CME structure on a scale that best identifies 
the CME front. The specific implementation of multiscale filtering intro- 



duced by Young & Gallagher (2008) also allowed the chaining of pixels 



along edges in the decomposed images to reveal the CME front for tracking 
through time. Such an algorithm provides a robustness in CME front detec- 
tion that alleviates issues of subjective user biases and unreproducibility of 
results. The technique is also more accurate than running-difference tech- 
niques which are widely used for determining CME heights, both spatially 
since it requires no arbitrary scaling and/or thresholding to find the edges, 
and temporally as it operates on individual images without a need for sub- 
tracting antecedent frames. The technique was also extended for use as a 
potentially automated CME front detection algorithm, discussed further in 
Section 16.2.11 below. 

2. An ellipse characterisation of the CME front in coronagraph data was im- 
plemented and shown to be an effective method for retrieving information 
on the changing CME morphology through sequences of observations. An 
ellipse was chosen for the innate freedom in its parameters, having the abil- 
ity to best fit the varying curvature of CME fronts across different image 
sequences, while still being constrained to close upon itself. The character- 
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isation provided a robust method for obtaining CME heights since it is not 
affected by deviations along the often ill-defined and/or kinked CME front. 
The opening cone angle to the ellipse also provided a measure of angular 
width for testing CME expansion, and the eccentricity of the ellipse pro- 
vided information on the changing CME front curvature as it was observed 
to propagate across the plane-of-sky. 

3. The degree of accuracy provided by the multiscale methods and charac- 
terisation of the CME front allowed a test of the constant acceleration 
model upon a variety of CMEs, and it was found not to be true of all 
cases. This has implications for how CMEs may be modelled theoretically 
since the forces acting must have different regimes of dominance to cause 
non-constant acceleration in certain events. This warrants further investi- 
gation, especially if projection effects can be corrected for using STEREO 
data, and the kinematics determined with the best available accuracy (see 



the discussion of Sections 6.2.3 and 6.2.4 below). 



4. The accuracy of these methods further revealed an early acceleration phase 
for some CMEs, and those with high speeds tended to reach greater angular 
widths, indicative of an initially dominant outward driving force ramping 
up the CME speed and expansion within the first few solar radii. Testing for 
this outward force that drives the CME acceleration and causes it to expand, 
is vitally important for comparing with theory in an effort to understand 
the interplay of forces acting on the CME as it propagates through the 
corona and heliosphere. 

5. The newly developed elliptical tie-pointing technique was shown to over- 
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come the plane-of-sky projection effects of previous single vantage-point 
observations in a manner better suited to studying the CME of 12 De- 
cember 2008 than previous stereoscopic efforts. By using the characterised 
curvature of the CME front in the observations from the twin-viewpoints 
of STEREO it was possible to constrain a reconstruction of the true CME 
front between the two planes-of-sky. This revealed the CME's true 3D mo- 
tion, thus greatly increasing the accuracy with which we can interpret its 
kinematics and morphology. 

6. The acceleration phase of the prominence eruption and CME was deter- 
mined in 3D to occur within the first 3 R , with the speed of the outer 
CME front being higher than the speed of the prominence that forms the 
inner core material of the CME. Their different speeds are indicative of the 
CME expansion, and indeed the early phase acceleration corresponds to 
the event's strongest characterised (out-of-ecliptic) angular width increase. 
This is again indicative of an initially dominant outward driving force that 
causes the CME speed and width to increase dramatically (generally consid- 
ered to be the Lorentz magnetic force of CME initiation and propagation). 

7. The deflected motion of the prominence eruption and CME was quantified 
in 3D; originating at high solar latitudes of ~ 55° but determined to move 
on a trajectory along the ecliptic at ~0°. This is indicative of the drawn- 
out background magnetic field of the quiet Sun. The pressure of the solar 
wind acts to drag out field lines in the (3 > 1 coronal regime, influencing 
the overall magnetic field configuration right down to the surface of the 
quiet Sun. At these low heights the magnetic pressure of the drawn-out 
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field guides the motion of the gradual eruption into a non-radial trajectory, 
placing it on a course towards Earth. This highlights the importance of the 
overall solar magnetic field with respect to the CME and its source region, 
and reveals how necessary it is to understand their interdependence, both 
in the context of CME physics and also space weather monitoring. 

8. The effects of drag on the 3D CME motion were investigated and it was 
found that, subsequent to the early acceleration phase of the CME, its speed 
continued to increase towards the speed of the ambient solar wind. A com- 
parison of the CME velocity with the in-situ detection of its arrival at LI 
implied the CME had to be further slowed along its trajectory, revealed to 
be true when investigated with the 3D MHD 'ENLIL with Cone' Model. 
This highlights the importance of drag on the CME throughout its prop- 
agation as it was found to undergo further deceleration by the slow-speed 
solar wind stream ahead of it as it propagated through interplanetary space. 
Understanding these drag effects is thus of great importance for predicting 
CME arrival times at Earth and improving space weather forecasts. 

6.2 Future Work 

The methods developed and implemented in this thesis are a first use of multiscale 
analysis and characterisation in obtaining CME kinematics and morphology with 
better constrained errors than previous efforts, of benefit in testing theoretical 
CME models and in forecasting space weather at Earth or elsewhere in the helio- 
sphere. The possibility for automation has been demonstrated and could provide 
a new catalogue of CMEs working in realtime with SolarMonitor.org for exam- 
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pie. Extending these multiscale techniques to curvelets or ridgelets may better 
suit the detection of the typically curved nature of CMEs in coronagraph data. 
Furthermore, the data from the twin viewpoints of the STEREO mission have 
allowed us to overcome plane-of-sky projection effects in studying CMEs with 
our newly developed elliptical tie-pointing technique, resulting in a very cohesive 
description of how the CME of 12 December 2008 propagates from the Sun to the 
Earth. There are numerous candidate events which should also be studied in this 
manner to test the conclusions drawn from the results of the CME studied here 
and gain further insight into CME dynamics. Other methods for deriving the 
velocity and acceleration profiles of CMEs also warrant investigation, since the 
3-point Lagrangian interpolation is sensitive to scatter in the data (though less 
so than the standard forward/reverse difference) and bootstrapping or inversion 
techniques, for example, may help overcome this. 

6.2.1 Automation 

As discussed in Chapter [3] there is great benefit in implementing an automated 
CME detection and tracking algorithm for cataloguing their kinematics and simi- 
lar properties of interest in large data sets. The algorithms (specifically CACTus, 
SEEDS and ARTEMIS) remove the need for a manual inspection of the images 
(as performed in the CDAW catalogue), which can be both laborious and suffer 
from user-specific biases. The automated techniques also produce a robust output 
of parameters for large statistical analyses of CME properties over long periods of 
solar activity, since the thresholds for detection and height /width measurements 
are hard-coded into the algorithm. However, the ultimate aim is to determine 
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Figure 6.1: An example interface of the automated multiscale detection algorithm 
based upon thresholding the magnitude and angular information, discussed in Sec- 
tion |3.3| Left: The algorithm can be made to cycle through different scales, and 
threshold any number of contoured regions of magnitude (edge strength). Right: 
Contours that contain a wide distribution of angles signal a CME detection. The 
angular information is normalised to 1 and folded into a 0-180° range due to 
symmetry of the edge normals. A threshold on the normalised angular distribu- 
tion is specified to flag regions as CMEs or otherwise (e.g., >20%). A detection 
mask is then built through multiple scales. The limitations to be overcome for a 
robust automation of this technique include developing dynamic thresholds such 
that multiple contours of CME edges are not fragmented, and increased angular 
distributions due to streamer deflections are not mistakenly labelled as CMEs (a 
non-trivial task). Moreover, an automated pixel chaining of the CME edges must 
be included in order to produce an ellipse characterisation of the CME front in the 
image. 
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Figure 6.2: An example interface which could be developed for the potentially 
automated multiscale filtering and ellipse tracking of CMEs observed by STEREO- 
A and B. Images from the respective instruments appear on the left, and the 
resulting parameters from the characterisation of the CME front appear on the 
right. Any manner of information may be chosen for display, and clocked through 
time as the CME progresses. 
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the kinematics and morphology of CMEs with as great an accuracy as possible 
with the available data, thus operations such as image rebinning, smoothing, and 
differencing are not ideal. Also, accurately measuring the CME height (and sub- 
sequently deriving the velocity and acceleration) can be difficult since the CME 
front is often diffuse and ill-defined. 

These issues motivated the application of multiscale filtering to enhance the 
CME front in coronagraph images without reducing the quality of the image since 
noise and small-scale features can be removed in the multiscale decomposition for 
optimum CME detection. The multiscale technique outlined in Chapter [3] also re- 
sults in an edge detection on the image akin to the Canny edge detector as a result 
of the horizontal and vertical directions in which the multiscale filter is applied. 
Regions with intensities above a specific threshold in the image are contoured 
as possible CME candidates, though this often includes streamers since they ap- 
pear on the same scales as CMEs. CME detection is then shown to be effective 
through the thresholding of angular spreads across the contoured regions of the 
image, since CMEs will have a large distribution of edge normals compared with 



the generally radial nature of coronal structures such as streamers (Figure 6.1). 
Ideally the thresholds should be made dynamic so as to better detect the varying 
CME intensities and minimise the fragmentation of intensity contoured regions 
of interest. Too strong a hard threshold will neglect faint CMEs or parts thereof. 
The algorithm also currently lacks a satisfactory technique for discarding the re- 
gions of an image that do not correspond to the detected CME and thus still 
requires a user to perform an inspection of the edges to be maintained/discarded 
in the images. Including multiple scales of the decomposition, rather than just 
the one with the best signal-to-noise ratio for the CME, allows a scoring system 
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Pixels 




Figure 6.3: The magnitude (edge strength) from the multiscale filtering of a CME 
observation, unwrapped into polar coordinates (r, 6) with axis units of pixels and 
degrees, respectively. A scanning-line runs over the image to produce the scatter 
of edge normals at the bottom of the image, where each normal is plotted with 
a vector magnitude and associated angular information from the multiscale filter 



(cf. the edge normals of Figure 3.6). An end-projection of the normals along 
the scanning-line is plotted in the bottom right of the image. The scanning-line is 
located at angle 328°, passing over a CME in the image. The resulting edge normals 
show a slice through the angular distribution of the CME, and as the scanning- line 
moves along the image, the end-projection shows a continuous rotation of angles 
due to the curvature of the CME front. Detecting this continuous rotation of angles 
may be used to distinguish CMEs from streamers which show an abrupt, stepwise 
change in angle as the scanning-line crosses over them. This may alleviate a need 
for the intensity thresholding and contouring which can fail to detect faint CMEs 
or portions thereof. 
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to refine the CME detection mask, although cases of streamer interaction or de- 
flection are prone to skew the resulting detection. The algorithm needs to be 
improved with a form of image segmentation included to make the CME detec- 
tion masks more robust, and minimise the effects of streamers. Then the ellipse 
characterisation can be automatically applied for studying the CME front kine- 
matics and morphology with greater accuracy than traditional image processing 
techniques, and produce a multiscale-methods based catalogue of events (e.g., 



Figure 6.2). 

In an effort to overcome the hard-thresholding biases of the intensity con- 
touring of CMEs, a detection based solely on the angular distribution may be 
considered in future work. First a multiscale filtered coronagraph image, like 



that of Figure 3.8, is unwrapped into polar coordinates (x, y) — > (r, 9). Then the 
image is scanned with regard to the combined magnitude and angular information 
of the multiscale filtering process, resulting in a dynamic change of vector arrows 



(edge normals), as shown in Figure 6.3 When scanning over the quiet corona or 
streamers, their radial nature means the angles along the scanning-line will all 
predominantly point in one direction and undergo abrupt, stepwise changes. As 
the scanning-line moves across a CME the angles will change in a more linear, 
continuous fashion. When such a continuous angular change is measured, the 
region under the scanning-line can be flagged as a CME detection. 



6.2.2 Ridgelets/Curvelets 

The implementation of multiscale analysis has been demonstrated for its efficacy 
in highlighting CMEs and coronal structure against noise and small-scale features 



193 



6. CONCLUSIONS & FUTURE WORK 




in coronagraph images (Byrne et al., 2009 Stenborg & Cobelli, 2003 Stenborg 



et al, 2008 Young h Gallagher, 2008). However, wavelets are better suited 



to identifying point-like features in images, but their extension to ridgelets and 
curvelets has been shown to better resolve the visibility of the curved form of a 



typical CME front (Gallagher et al. , 2010 ). Thus their development may provide a 



more reliable CME detection algorithm than the multiscale technique investigated 
in this work. 

The continuous wavelet transform of an image may be defined as: 



w (s, a, b) 



I (x, y) ip 87ajb (x,y) dxdy 



(6.1) 
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Figure 6.5: The curvelet filtered image of a CME observed by LASC0/C2 on 18 
April 2000, reproduced from Gallagher et al. (2010). The detail along the curved 



CME front is enhanced as a result of the curvelet technique following the removal 
of certain coefficients probably due to noise. 



where w(s,a,b) are the wavelet coefficients of the image I(x,y), ip s (x,y) is the 
mother wavelet, and s is the term describing scale at a position (a, b). The 
mother wavelet can take many forms, e.g., the Morlet wavelet or the Mexican hat 
wavelet. The ridgelet transform takes a similar mathematical form to the wavelet 
transform, i.e., it is a convolution of an image with a predefined basis function, 



but they are anisotropic and thus more directionally sensitive (Figure 6.4). The 
ridgelet uses a radon transform that transforms lines into points, upon which a 
wavelet transform may be applied to provide a sparse representation of the points. 
The basis function of the ridgelet takes the form: 



-1/2 , i'xcos6 + ysin9 -b\ 
i>8,b,o(x,y) = s i ip[ (6.2) 



The ridgelet is constant along lines x cos 9 + y sin 9 = const, and the ridgelet 
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coefficients are given by the convolution: 



Ri (s,b,9) 




I (x, y) ip s , bfi (x, y) dx dy 



(6.3) 



Curvelets generalise the idea of ridgelets to multiscale curves in images. The 
detailed maths on the application of curvelets to CME images is discussed in 



Gallagher et al. (2010) based on the developments by Candes & Donoho (1999). 



Figure 6J3 shows how well the curvelet filter performs on an image of a CME. It 
enhances the structure within the curved CME front while reducing the noise and 
small-scale intensity features elsewhere in the image. This may provide a more 
accurate method of determining the CME front location in individual images, 
rather than relying on differencing techniques which can be prone to scaling 



errors and spatio-temporal cross-talk (discussed in Section 4.1). 



6.2.3 3D CME Reconstruction 

The development of the elliptical tie-pointing technique has proven effective at 
reconstructing the CME front of the 12 December 2008 event, discussed in detail 
in Chapter [5} Such studies are necessary for obtaining the true kinematics and 
morphology of CMEs and minimising the uncertainties in their 3D quantifica- 



tion, of benefit in comparing observations with theoretical models (Figure 6.6). 
This event was an ideal case for study since the STEREO spacecraft were at an 
angular separation of almost 90°, so the lines-of-sight intersected to form opti- 
mum quadrilaterals localising the CME front. Very small, or very large, angles of 
separation would not be as effective since the quadrilaterals will be skewed and 
so too will the corresponding inscribed ellipses. So while a reconstruction would 
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Figure 6.6: A schematic of the possible 3D magnetic field topology of a CME, 
as inferred from typical observations of a leading edge, dark void and bright core. 
There may also be an associated on-disk post-flare arcade at the base of the pos- 
tulated current sheet. 3D reconstructions of CMEs are important for testing the 
validity of theoretical models. 
Image credit: http://www.nrl.navy.mil 

match the observations of both spacecraft it may not perfectly represent the true 
curvature of the CME front (although it will still offer a better approximation 
than tie-pointing alone). 

The fact that the 12 December 2008 CME propagated along the Sun-Earth 
line midway between the two spacecraft meant the reconstruction could be per- 
formed out to heights of almost 50 R Q before the 'pancaking' of the CME (along 
with concerns on the scattering geometry and instrumental effects) meant the 
ellipse fit was no longer appropriate. If we can correct for these effects, or in- 
clude them in the uncertainty, then it is possible to use more than one ellipse 
fit to characterise the different portions of the CME front and perform ellipti- 
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cal tie-pointing on these individually. This could provide an insight into the 3D 
propagation through the rest of the HI1 field-of-view. It would also be interesting 
to test how the method fares with the observations of HI2. 

A number of other CMEs will warrant studying with the elliptical tie-pointing 
technique throughout the lifetime of the STEREO mission. Many events have 
been observed through the COR1 and COR2 fields of view and been studied by 



several authors through a variety of stereoscopic methods (Boursier et al., 2009a 



de Koning et al. 


2009 


Liewer et al. 


2009 


Mierla et al., 


2008 


Srivastava et al. 



2009 Temmer et al., 2009 Wood et al., 2009). Direct comparisons may be made 



with the results of these studies and future stereoscopic CME analyses as the 
spacecraft move into orbits on the far side of the Sun, approaching quadrature 
again as the Sun becomes more active through the rise of solar cycle 24. In con- 
junction with further 3D CME reconstruction analyses, it will also be useful to 
perform a simulation of different model CME structures, with varying signal-to- 
noise ratios, cadences, image compressions, etc., and at different spacecraft sep- 
aration angles. Such simulations will provide a method for testing the inherent 
uncertainties associated with the different events being reconstructed at different 



stages of the STEREO mission. Furthermore, as discussed in Section 5.3.1, the 
observations from SOHO/LASCO can also be used as a third perspective on the 
events. 



In particular, the polarisation technique of Moran & Davila (2004) could be 
used in conjunction with the elliptical tie-pointing reconstruction, to potentially 
reveal detail of the structure of the CME behind the front. The technique relies 
upon the geometric dependence of the polarisation of Thomson-scattered light, 
whereby the polarisation fraction in CME emission provides a line-of-sight av- 
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Figure 6.7: 'Top' view of the polarimetric reconstruction method applied to both 
STEREO-A and B observations of a CME at 01:30 UT on the 31 December 2007 



(rotated to the Sun-Earth frame), reproduced from Moran et al. (2010). The red 
points are the reconstruction from STEREO-A, and the blue points from STEREO- 
B, with the green points showing the regions of overlap. 



eraged 'mean distance to the plane of the sky' for selected, or all, CME pixels. 
The validity of the method was tested and proven using STEREO observations of 
two CMEs and shown to be in good agreement with other triangulation methods 



(Moran et al. , 2010). A combined use of the technique with the geometric local- 



isation on COR2 beacon data has been explored by De Koning et al. (2009), al- 



though the polarisation technique becomes unreliable at heights > 5 R (Billings 



1966). An example of how the technique applies to STEREO data is shown in 



Figure 6/7 for a CME on 31 December 2007. It can be seen how the final polari- 
metric reconstruction could be combined with the 3D CME front reconstruction 
for a more cohesive rendering of the overall CME structure. 

Due to the nature of the STEREO mission, telemetry limitations are an im- 
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portant point to note as the spacecraft separate, specifically regarding image 
degradation and the consequences for CME studies. The overall mission is split 
into three phases: early operations, prime science phase, and extended mission 
phase. During the early operations phase, when the spacecraft were still close 
to the Earth, telemetry rates were restricted compared to the rates of the prime 
science mission, being 30 kbps in the first two weeks. Once the spacecraft were in 
heliocentric orbits the prime science phase began and nominal telemetry rates of 
720 kbps were employed. As the spacecraft separate, the data rate is diminished, 
going down to 160 kbps in early 2011, and 120 kbps by 2013. In the case of the 
SECCHI suite, progressive telemetry reductions are introduced through increased 
image binning and compression, reduced image cadences, and prioritisation of to- 
tal brightness coronagraph images, as per the STEREO Science Operations PlarQ 
With regards the elliptical tie-pointing technique, the uncertainties involved in 
identifying and tracking CMEs will increase due to the image binning and reduced 
cadences. When the spacecraft are in near-quadrature (the optimum positioning 
for the technique) on the far side of the Sun again in late 2012 to early 2013 the 
low-rate telemetry stream will limit the precision with which the true kinematics 
of future events may be derived. 



6.2.4 Deriving CME Kinematics 

The standard method for determining the kinematics of a CME is to obtain a 
sequence of height-time measurements and perform a 3-point Lagrangian inter- 
polation to derive the velocity and acceleration of the event (detailed in Sec- 
tion 4.2.1). This method alone is somewhat dated since the advent of strong 



^ttpi/Zstereo-ssc. nascom.nasa.gov/publications.shtml 
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Figure 6.8: A theoretical model for a CME with constant acceleration 2 m s 
and initial velocity 300 km s _1 , and two simulations of how the resulting profiles 
for a noisy sample of datapoints behave using 3-point Lagrangian interpolation. 



numerical computing power and development of bootstrapping and spline-fitting 
techniques, for example. As a first approximation the 3-point Lagrangian pro- 
vides a good estimate of how the velocity and acceleration corresponding to a 
height-time curve are likely to behave, by revealing the trends in the profiles that 
indicate increasing/decreasing velocity and/or acceleration. However, this is only 
true if the scatter and errors of the height measurements are not unreasonably 
large, since a large scatter would be enhanced by the derivatives, and a large 
error will increase the uncertainty on the derivative points, such that trends may 
become untestable. A quick simulation of how the 3-point Lagrangian fares for a 
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Figure 6.9: The resulting velocity profile for the 3D reconstructed CME front 
of the 12 December 2008 using the inversion technique of Kontar & MacKinnon 



(2005). 



theoretical model CME undergoing a constant acceleration of 2 m s~ 2 and initial 
velocity of 300 km s _1 reveals the unreliability of the resulting kinematic profiles. 
A scatter of height-time datapoints is chosen with varying levels of noise up to 
~20%. An errorbar on each datapoint is determined by its distance from the 
theoretical height-time profile. Various instances of datapoint scatters result in 
erroneous trends in the velocity and acceleration profiles - even with the proper 



error treatment. Figure [678] shows two examples of how different the derived kine- 
matics can be, in comparison with the theoretical model and each other. They 
show how different scatters of the datapoints can result in what appear to be 
completely opposing acceleration trends, meaning the nature of the scatter is not 
satisfactorily reflected in the derived errorbars of the resulting profiles. This war- 
rants further investigation for future CME, and other, kinematic studies, since 
the implications for how we interpret observations could be profound. 

It should also be noted that the 3-point Lagrangian counter-intuitively in- 
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creases the errorbars on the resulting velocity and acceleration profiles when the 
number of height-time measurements are increased (i.e., for smaller cadences 
which new missions continue to provide). This is due to the algorithm's inverse 
dependence on the spacing between points. It is therefore worth investigating 
how other techniques might be applied to more confidently derive the kinemat- 
ics of CMEs. Bootstrapping of a presumed model fit to the CME height-time 
profile would result in the best match of parameters to the data, but questions 
would remain on the appropriateness of the chosen model itself. Alternatively a 
simulation of data could be bootstrapped regarding the resulting derived kine- 
matics and an estimate of the errors involved could be deduced to apply to true 
observations, though this again may be model dependant. An inversion tech- 



nique could also be investigated for obtaining derivatives (Kontar & MacKinnon 



2005), an example of which is shown in deriving the velocity of the 12 December 



2008 CME (Figure 6.9). It works by essentially solving for the smoothest spline 
fit that minimises the distance between the end-points while still being bound 
by any constraints on the data. If a quantification of the height-time errors is 
provided, then inversion techniques may result in a more robust determination of 
the kinematic profiles with statistically sound uncertainties. 
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